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Fiir meinen Vater 



Granger stood looking back with Montag. "Everyone must leave something behind when he 
dies, my grandfather said. A child or a book or a painting or a house or a wall built or a 
pair of shoes made. Or a garden planted. Something your hand touched some way so your 
soul has somewhere to go when you die, and when people look at that tree or that flower you 
planted, you're there. It doesn't matter what you do, he said, so long as you change something 
from the way it was before you touched it into something that's like you after you take your 
hands away. The difference between the man who just cuts lawns and a real gardener is in the 
touching, he said. The lawn-cutter might just as well not have been there at all; the gardener 
will be there a life-time." 



Ray Bradbury, Fahrenheit 451 
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Preface 



This thesis has been divided into two parts both being apphcations of perturbative Quan- 
tum Chromo Dynamics (QCD). 

The first part 'Heavy Quark Production in CC and NC DIS' is a continuation of work 
presented in [1] and has been done in collaboration with S. Kretzer. Part I is based on 
the following pubhcations [2-4]: 

• S. Kretzer and I. Schienbein, Charged- current leptoproduction of D mesons in the 
variable flavor scheme, Phys. Rev. D56, 1804 (1997) [Chapter 2]. 

• S. Kretzer and I. Schienbein, Heavy quark initiated contributions to deep inelastic 
structure functions, Phys. Rev. D58, 094035 (1998) [Chapter 3]. 

• S. Kretzer and I. Schienbein, Heavy quark fragmentation in deep inelastic scattering, 
Phys. Rev. D59, 054004 (1999) [Chapter 4]. 

The second part 'The Structure of Real and Virtual Photons' has emerged from a 
collaboration with Prof. M. Gliick and Prof. E. Reya. Part II is based on the following 
pubhcations [5-8]: 

• M. Gliick, E. Reya, and 1. Schienbein, Pionic parton distributions revisited, Eur. 
Phys. J. CIO, 313 (1999) [Chapter 10]. 

• M. Gliick, E. Reya, and I. Schienbein, Radiatively generated parton distributions for 
real and virtual photons, Phys. Rev. D60, 054019 (1999); (E) D62 019902 (2000) 
[Chapter 11]. 

• M. Gliick, E. Reya, and I. Schienbein, The Photon Structure Function at small-x, 
Phys. Rev. D64, 017501 (2001) [Chapter 11]. 
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• M. Gliick, E. Reya, and I. Schienbein, Has the QCD RG-improved parton content 
of virtual photons been observed?, Phys. Rev. D63, 074008 (2001) [Chapter 12]. 

The results of Refs. [5,6] have also been summarized in [9, 10]. The work in Section 7.3 
is completely new and has been done in collaboration with C. Sieg. 



Part I 

Heavy Quark Production in CC and 

NC DIS 



Chapter 1 

Introduction and Survey 



Part I of this thesis is devoted to heavy quark production in charged current (CC) and 
neutral current (NC) deep inelastic scattering (DIS) where a quark h with mass nih 
is viewed as heavy if nih ^ Aqcd contrary to the light -u, rf, s quarks with rriu^dys <^ 
Aqcd- Thus, the heavy quark mass provides a hard scale allowing for a perturbative 
analysis. In deep inelastic heavy quark production a second hard scale is given by the 
virtuality of the probing boson (photon, Z-boson, W-boson) such that we have to deal 
with the theoretically interesting problem of describing a two-hard-scale process within 
perturbative QCD (pQCD). In the following, we will prominently deal with charm quarks 
being the lightest of the heavy quarks {h — c, b, t) in the standard model. 

In addition to these general theoretical arguments, there are several important phe- 
nomenological reasons for studying heavy quark production in DIS: 

• Charm contributes up to 30% to the total structure function at small Bjorken-x 
as measurements at the ep coUider HERA at DESY have shown [11-13]. For this 
reason a proper treatment of charm contributions in DIS is essential for a global 
analysis of structure function data and a precise extraction of the parton densities 
in the proton. 

• NC charm production offers the possibility to extract the gluon distribution in the 
proton from a measurement of F| [14, 15]. 

• CC charm production is sensitive to the nucleon's strange sea. The momentum 
(z) distributions of D-mesons from the fragmentation of charm quarks produced in 
neutrino deep inelastic scattering have been used recently to determine the strange 



quark distribution of the nucleon s{x, Q'^) at leading order (LO) [16] and next-to- 
leading order (NLO) [17]. 

In the past few years considerable effort has been devoted to including heavy quark 
effects in DIS. If one could sum the whole perturbation series (keeping the full mass 
dependence) one would arrive at unique perturbative QCD predictions. However, at any 
finite order in perturbation theory differences arise due to distinct prescriptions (schemes) 
[18-25] of how to order terms in the perturbative expansion. These schemes, in turn, enter 
global analyses of parton distributions [26-28] and hence feed back on the resulting light 
parton distributions [28]. It is therefore necessary to work out the various schemes as 
much as possible and to compare them since a most complete possible understanding or 
-even better- reduction of the theoretical uncertainties is required to make concise tests 
of pQCD predictions against heavy flavor tagged deep inelastic data. 

In this thesis we concentrate on the ACOT variable flavor number scheme (VFNS) 
[19,20] which wc work out to full order 0{al). Recently, this scheme has been proven 
by Collins [21] to work at all orders of QCD factorization theory. We will perform all 
required calculations with general couplings and masses in order to be able to describe 
both CC and NC processes in one framework. 

The outhne of Part I will be as follows: 

• In Chapter 2 we present formulae for the momentum (z) distributions of D-mesons 
produced in neutrino deep inelastic scattering off strange partons. The expressions 
are derived within the variable flavor scheme of Aivazis et al. (ACOT scheme) [20], 
which is extended from its fully inclusive formulation to one-hadron inclusive lep- 
toproduction. The dependence of the results on the assumed strange quark mass 
nis is investigated and the — > limit is compared to the corresponding MS re- 
sults. The importance of 0{as) quark-initiated corrections is demonstrated for the 
nis — case. 

• In Chapter 3 we perform an explicit calculation of the before missing 0{al) correc- 
tions to deep inelastic scattering amplitudes on massive quarks within the ACOT 
scheme using general masses and couplings thereby completing the ACOT formalism 
up to full O{os). After identifying the correct subtraction term the convergence of 
these contributions towards the analogous coefficient functions for massless quarks. 



obtained within the modified minimal subtraction scheme (MS), is demonstrated. 
Furthermore, the importance of these contributions to neutral current and charged 
current structure functions is investigated for several choices of the mass factoriza- 
tion scale fx as well as the relevance of mass corrections. 

In Chapter 4 we turn to an analysis of semi-inclusive production of charm (mo- 
mentum (z) distributions of D-mesons) in neutral current and charged current deep 
inelastic scattering at full 0{al). For this purpose we generalize the results of Chap- 
ter 3 from its fully inclusive formulation to one-hadron inclusive leptoproduction. 
We review the relevant massive formulae and subtraction terms and discuss their 
massless limits. We show how the charm fragmentation function can be measured in 
CC DIS and we investigate whether the charm production dynamics may be tested 
in NC DIS. Furthermore, we also discuss finite initial state quark mass effects in 
CC and NC DIS. 

Finally, we summarize our main results in Chapter 5. Some details of the calculation 
and lengthy formulae are relegated to the Appendices A and B. 



Chapter 2 



Charged Current Leptoproduction of 
D— Mesons in the Variable Flavor 
Scheme 



2.1 Introduction 



The momentum (z) distributions of D-mesons from the fragmentation of charm quarks 
produced in neutrino deep inelastic scattering (DIS) have been used recently to deter- 
mine the strange quark distribution of the nucleon s{x, Q^) at leading order (LO) [16] 
and next-to-leading order (NLO) [17]. A proper QCD calculation of this quantity re- 
quires the convolution of a perturbative hard scattering charm production cross section 
with a nonperturbative c — > D fragmentation function Dc{z) leading at 0{as) to the 
breaking of factorization in Bjorken-x and z as is well known for light quarks [29,30]. 
So far experimental analyses have assumed a factorized cross section even at NLO [17]. 
This shortcoming has been pointed out in [31] and the hard scattering convolution kernels 
needed for a correct NLO analysis have been calculated there in the MS scheme with three 
massless flavors {u, d, s) using dimensional regularization. In the experimental NLO anal- 
ysis in [17] the variable flavor scheme (VFS) of Aivazis, CoUins, Olness and Tung (ACOT) 
[20] for heavy flavor leptoproduction has been utilized. In this formalism one considers, 
in addition to the quark scattering (QS) process, e.g. W~^s c, the contribution from 
the gluon fusion (GF) process W'^g — > cs with its full ms-dependence. The coUinear 
logarithm which is already contained in the renormalized s{x, Q^) is subtracted off nu- 
merically. The quark-initiated contributions from the subprocess W~^s — > eg (together 
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with virtual corrections) which were included in the complete NLO (MS) analysis in [31] 
are usually neglected in the ACOT formalism. The ACOT formalism has been formu- 
lated explicitly only for fully inclusive leptoproduction [20]. It is the main purpose here 
to fill the gap and provide the expressions needed for a correct calculation of one-hadron 
(D-meson) inclusive leptoproduction also in this formalism. 



2.2 Gluon-Fusion 



In the following we will stick closely to the ACOT formalism as formulated in [20] ex- 
cept that we are not working in the helicity basis but prefer the standard tensor basis 
implying the usual structure functions Fi=i,2,3- We are not considering kinematical effects 
arising from an initial state quark mass in the W'^s — > c quark scattering contribution, 
i.e., s{x, Q"^) represents massless initial state strange quarks. This latter choice must be 
consistently kept in the subtraction term [20] to be identified below from the ^ limit 
of the W'^g — > cs gluon fusion contribution. The fully massive partonic matrix elements 
have been calculated for the general boson-gluon fusion process Bg — > Q1Q2 in [32] where 
B — W"^, Z. When they are convoluted with a nonperturbative gluon distribution 
g{x,ii^) and a fragmentation function Dq^{z), one obtains the GF part of the hadronic 
structure function Fi{x, z,Q^) describing the momentum (z) distribution of a hadron H 
containing the heavy quark Q2: 

'max[z,Cmin{x/x')] 

dx' r^rna.{x/x') 

'max[z,Cminix/x')] 

with the fractional momentum variables z = pn -Pn/i-Pn and C = VQi -Pn/q-Pn, Pn and 
q being the momentum of the nucleon and the virtual boson, respectively. The structure 
functions Fi{x^z^Q'^) generalize the usual fully inclusive structure functions Fi^x.Q"^), if 
one considers one-hadron (H) inclusive leptoproduction. The partonic structure functions 
Mx', C, Q'^) are given by 



J ax X Jmax[zXmin{x/x')] S •'^ S 

F,«^(a;,^,g^) = /^/ ^a;'5(a;',/.^)/2(4,C,g^)/^Q.(^) (2.1) 

J ax J max\z,Cmin(x/x')] S S 



/i=l,2,3 {x'X,Q^) 



TT 



A, Bi C, D 



(2.2) 
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with 



Ci(x',g2) 



E2 

C3 

-B- 1,2 
3 

D. 1,2 

3 



X 



/2 2 

X mj 



= -r 7^ 1 + 



4 Q 



J 



0+ 

4 



__^'(l_^')____(l_2x') + 



/ Am^ X 



2^/ \ 2 



^ ( -1 + 2x' - ) 
4 

g+ x 



i2 



X 



ml (if ^rri 



.2\ 2 



A/;r — mf 1 

^2 



2 

1 - 2x'{l -x') + ^(l + 8x' - 18x'^) 



+ ■^[l-4x' + 6x"j 2x'{l - 3x') + 4x'(l - 5x') 



+ 



) Q 

Am^Am^ ,2 <Z- 2mim2 



^ q+x 



.2t' - — 

l + 3x'{l-x') 
,2 Am^ + g^ 



R„ 



\-x\l-x') 



Am^ 



x'{l-2x') + 



Am^ 



„/2 



g^ 



a; 



\Q') = 

g^) = ±A (a:;', g^) [rni ^ ms] 
V,g2) = ±C3 (x', g^) [mi ^ m2] 



where Am" = m2 — m" , mi 2 being the mass of the heavy quark gi,2- The kinematical 
boundaries of phase space in the convolutions in Eq. (2.1) are 



ax — 



1 + 



(mi + m2)' 



with — 1 



(mi + m2)^ x' 



mm.max 



^ ' 2 



Am^ x' 
1 H — : ± vv 



g2 i-x' 



(2.3) 



v' = 1 



(mi - m2 



|2 x' 



1-x' - ^ g^ l-x^ ■ The vector (y) 

and axialvector (A) couphngs of the 7ju(^ — ^75) quark current enter via q± — V'^±A'^, 
Rq — VA. ^ If the partonic structure functions in Eq. (2.2) are integrated over ( the well 
iPor the 7*Z interference term the couphngs read: q± = VW^ ± A^A^, Rg = 1/2{V^A^ + V^A'^). 



2.3 CC Leptoproduction within ACOT 



8 



known inclusive structure functions [33-35] for heavy flavor production are recovered: 

/ ^ dC fi^.,2{x\C.Q^)^±U{x\Q^) , (2.4) 
where the fi{x',Q'^) can be found in [33]. 

2.3 CC Leptoproduction within ACOT 

In the following we will consider the special case of charged current charm production, 
i.e., mi = rris, m2 = rric {q± = 2,0; i?g = 1 assuming a vanishing Cabibbo angle). Of 
course, all formulae below can be trivially adjusted to the general case of Eqs. (2.1) and 
(2.2). The nis limit of the partonic structure functions in Eq. (2.2) is obtained by 
keeping terms up to O(m^) in the Ai, Ci and in Cmax due to the singularity of the phase 
space integration stemming from C ^ 1- One obtains 

hm^ ^/,(a;',C,g') = Q i/f (^, C, m^. A) 



*o a, A 



= Q 5(1-0 i^J)(y) In ^^ + 0(m°) (2.5) 
A mi 



where P^^^x') = -[x'^ + (1 - x'% A = gV(Q' + Ci = 1/2, = x'/\, c, = 1 
and the Hf are the dimensionally regularized MS (m^ = 0) gluonic coefficient functions 
obtained in [31]. The q arise from different normalizations of the fi and the and are 
such that the infrared-safe subtracted [see Sec. 2.4] convolutions in Eq. (2.1) converge 
towards the corresponding ones in [31] as — > if one realizes that x'/X = x/X ~ 
^ = x{l + m'^/Q'^). Taking also the limit m^ — > in Eq. (2.5) gives -besides the coUinear 
logarithm already present in Eq. (2.5)- finite expressions which agree with the massless 
results of [29,30]. 

In the ACOT formalism the GF convolutions in Eq. (2.1) coexist with the Born level 
quark scattering contributions F^^{x,z,Q'^) — ki s{^,fj,'^) Dc{z), A;i=i,2,3 = 1, 2^, 2. The 
overlap between the QS and the GF contributions is removed by introducing a subtraction 
term (SUB) [20] which is obtained from the massless hmit in Eq. (2.5) 



ZTT m'^ 



s 



J |^9(x'./'^X'(| 
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The complete 0{as) structure functions for the z distribution of charmed hadrons (i.e., 
dominantly D-mesons) produced in charged current DIS are then given in the ACOT 
formahsm [20] by 

pACOT _ pQS _ pSUB _^ pGF ^2 7) 

It is worthwhile noting that general results for polarized partonic structure functions 
analogous to the unpolarized ones in Eq. (2.2) have been obtained in [36] allowing for 
a formulation of polarized one-hadron inclusive heavy flavor leptoproduction within the 
ACOT scheme along the same hues. 



2.4 Numerical Results 

In Fig. 2.1 we show the structure function F^'^'-''^ at experimentally relevant [16] values 
of X and for several finite choices of together with the asymptotic — > limit. 
For Dc we use a Peterson fragmentation function [37] 

D,{z) ^n[z[1- Z-' - £e/(l - z)fY' (2.8) 

with Be = 0.06 [38-40] normalized to dzD^{z) = 1 and we employ the GRV94(HO) 
parton distributions [41] with rric = 1.5 GeV. Our choice of the factorization scale is 
fi^ = + which ensures that there is no large ln(Q^ + m^) / //^ present in the difference 
GF-SUB. As can be seen from Fig. 2.1 the effects of a finite strange mass are small and 
converge rapidly towards the massless MS limit provided < 200 MeV as is usually 
assumed [17]. 

In Fig. 2.2 we show the effects of adding the quark-initiated 0{as) correction from 
the process W'^s — > eg (together with virtual corrections) to the asymptotic (m^ — > 0) 
F^^O'^, employing the CTEQ4(MS) densities [42] with = 1.6 GeV. The 0(q;,) quark 
contribution is usually neglected in the ACOT formahsm since it is assumed to be effec- 
tively suppressed by one order of ctg with respect to the gluon fusion contribution due 
to s{x,iJ?)/g{x,ii'^) ~ 0{as). To check this assumption for the quantity F2{x,z,Q'^) we 
show, besides the full result, the contributions from the distinct processes (using again 

^ Q2 _^ j^2^ r|.]^g y^+g _^ contribution corresponds to GF-SUB in Eq. (2.7). The 
quark-initiated 0{as) contribution has been calculated in the MS scheme according to [31] 
which is consistent with the asymptotic gluon-initiated correction in the ACOT scheme 
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Figure 2.1: The structure function F^^^'^{x, z,Q'^) as defined in Eq. (2.7) using the 
GRV94(HO) parton densities [41] with rric — 1.5 GeV and a Peterson fragmentation function 
[37] with Ec — 0.06. Several finite choices for are shown as well as the asymptotic — > 
limit. 

due to Eq. (2.5). It can be seen that the quark-initiated correction is comparable in size 
to the gluon-initiated correction around the maximum of F2. Since most of the exper- 
imentally measured D-mesons originate from this region the 0{as) quark contributions 
should not be neglected in a complete NLO calculation. 

2.5 Summary 

To summarize we have given formulae which extend the ACOT scheme [20] for the lep- 
toproduction of heavy quarks from its fully inclusive formulation to one-hadron inclusive 
leptoproduction. We have applied this formulation to D-meson production in charged 
current DIS and studied finite corrections to the asymptotic — > limit. The cor- 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 ^ 0.9 1 

Figure 2.2: The structure function F2{x,z,Q'^) for charged current leptoproduction of D- 
mesons at 0{as) for nis — using the CTEQ4(MS) parton distributions [42] and a Peterson 
fragmentation function [37] with Ec = 0.06. The full 0{as) result is shown as well as the 
individual contributions from the distinct quark- and gluon-initiated processes. 

rections turned out to be small for reasonable choices of rris < 200 MeV and we have 
shown that the rris limit reproduces the dimensionally regularized MS (m^ = 0) 
gluonic coefficient functions [31]. Furthermore we have investigated the quark-initiated 
0{as) corrections for = using the relevant MS fermionic coefficient functions [31]. 
The latter corrections turned out to be numerically important at experimentally relevant 
values of x and [16] and should be included in a complete NLO calculation of charged 
current leptoproduction of D-mesons. In the next chapter, the quark-initiated diagrams 
will be calculated up to in the ACOT scheme allowing to study finite effects 

stemming from these contributions as has been done in this chapter for the 0{as) gluon 
contributions. 



Chapter 3 



Quark Masses in Deep Inelastic 
Structure Functions 



In this chapter we consider heavy quark contributions to inclusive deep inelastic struc- 
ture functions within the ACOT variable flavor number scheme to which we contribute 
the calculation of the before missing Bremsstrahlung corrections (inch virtual graphs) 
off initial state massive quarks. The calculation exemplifies factorization with massive 
quark-partons as has recently been proven to all orders in perturbation theory by Collins 
[21]. After identifying the correct subtraction term the convergence of these contribu- 
tions towards the analogous coefficient functions for massless quarks, obtained within 
the modified minimal subtraction scheme (MS), is demonstrated. Furthermore, the phe- 
nomenological relevance of the contributions to neutral current (NC) and charged current 
(CC) structure functions is investigated for several choices of the factorization scale. The 
results presented in this chapter are taken from Ref. [3]. 



3.1 Introduction 



Leptoproduction of heavy quarks has become a subject of major interest in QCD phe- 
nomenology both for experimental and theoretical reasons. Heavy quark contributions are 
an important component of measured neutral current (NC) [11-13] and charged current 
(CC) [43] deep inelastic (DI) structure functions at lower values of Bjorken-x, accessible to 
present experiments. Charm tagging in NC and CC deep inelastic scattering (DIS) offers 
the possibility to pin down the nucleon's gluon [44] and strange sea [16, 17, 45, 46] density, 
respectively, both of which are nearly unconstrained by global fits to inclusive DI data. 
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Theoretically it is challenging to understand the production mechanism of heavy quarks 
within perturbative QCD. The cleanest and most predictive method [18] of calculating 
heavy quark contributions to structure functions seems to be fixed order perturbation 
theory (FOPT) where heavy quarks are produced exclusively by operators built from 
light quarks (u,d,s) and gluons (g) and no initial state heavy quark lines show up in any 
Feynman diagram. Heavy quarks produced via FOPT are therefore also called 'extrin- 
sic' since no contractions of heavy quark operators with the nucleon wavefunction are 
considered (which in turn would be characteristic for 'intrinsic' heavy quarks). Besides 
FOPT much effort has been spent on formulating variable flavor number schemes (VFNS) 
[19,20,22-24] which aim at resumming the quasi-collinear logs ^n{Q'^ /m^); Q and m be- 
ing the virtuality of the mediated gauge boson and the heavy quark mass, respectively] 
arising at any order in FOPT. All these schemes have in common that extrinsic FOPT 
induces the boundary condition [23, 47] 

q{x,Ql^m^)^Q + 0{al) (3.1) 

for an intrinsic heavy quark density, which then undergoes massless renormalization group 
(RG) evolution. Apart from their theoretical formulation VFNS have to be well under- 
stood phenomenologically for a comparison with FOPT and with heavy quark tagged DI 
data. We will concentrate here on the scheme developed by Aivazis, CoUins, Olness and 
Tung (ACOT) [19-21]. In the ACOT scheme full dependence on the heavy quark mass 
is kept in graphs containing heavy quark lines. This gives rise to the above mentioned 
quasi-collinear logs as well as to power suppressed terms of 0[{m'^ / Q'^Y]. While the latter 
give mass corrections to the massless, dimensionally regularized, standard coefficient func- 
tions (e.g. in the MS scheme), the former are removed by numerical subtraction since the 
coUinear region of phase space is already contained in the RG evolution of the heavy quark 
density. Up to now explicit expressions in this scheme exist for DIS on a heavy quark at 
0(q;°) [19] as well as for the production of heavy quarks via virtual boson gluon fusion 
(GF) at C(q;]) [20]. In Section 3.2 we will give expressions which complete the scheme 
up to 0{a\) and calculate DIS on a heavy quark at first order in the strong coupling, 
i.e. B*Qi — > Q2g (inch virtual corrections to B*Qi — > Q2) with general couplings of the 
virtual boson B* to the heavy quarks, keeping all dependence on the masses m\^2 of the 
quarks (5i,2- It is unclear whether (heavy) quark scattering (QS) and GF at 0{a].) should 
be considered on the same level in the perturbation series. Due to its extrinsic prehis- 
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tory QS^^) (bracketed upper indices count powers^ of ag) includes a coUinear subgraph of 
GF(2) e.g. eg contains the part of Yd ^ cc^f, where the gluon sphts into an almost 

on-shell cc pair. Therefore QS at 0{a\) can be considered on the level of GF at 0{al). 
On the other hand the standard counting for light quarks is in powers of ag and heavy 
quarks should fit in. We therefore suggest that the contributions obtained in Section 3.2 
should be included in complete experimental and theoretical NLO-analyses which make 
use of the ACOT scheme. Theoretically the inclusion is required for a complete renor- 
malization of the heavy quark density at 0(aj). However, wc leave an ultimate decision 
on that point to numerical relevance and present numerical results in Section 3.3. Not 
surprisingly they will depend crucially on the exact process considered (e.g. NC or CC) 
and the choice of the factorization scale. Finally, in Section 3.4 we draw our conclusions. 
Some technical aspects and lengthy formulae are relegated to the Appendix A. 

3.2 Heavy Quark Contributions to ACOT Structure 
Functions 

In this section we will present all contributions to heavy quark structure functions up 
to 0{a\). They are presented analytically in their fully massive form together with the 
relevant numerical subtraction terms which are needed to remove the coUinear divergences 
in the high limit. Section 3.2.1 and 3.2.3 contain no new results and are only included 
for completeness. In Section 3.2.2 we present our results for the massive analogue of the 
massless MS coefficient functions Cf^^. 

3.2.1 DIS on a Massive Quark at 0{al) 

The C(q;°) results for B*Qi — > Q2, including mass effects, have been obtained in [19] 

within a helicity basis for the hadronic and partonic structure functions. For completeness 

and in order to define our normalization we repeat these results here within the standard 

tensor basis implying the usual structure functions Fj=i^2,3- The helicity basis seems to be 

advantageous since in the tensor basis partonic structure functions mix to give hadronic 

structure functions in the presence of masses [19]. However, the mixing matrix is diagonal 

[19] for the experimental relevant structure functions -Fi=i,2,3 and only mixes with F5 

^For the reasons given here we refrain in this chapter in most cases from using the standard terminology 
of 'leading' and 'next-to-leading' contributions and count explicit powers of ««. 
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Figure 3.1: Feynman diagrams for the QS^^^ (a) and QS'^^'' [(b), (c)] contributions to ACOT 
structure functions in Eqs. (3.3) and (3.9), respectively. 

which are both suppressed by two powers of the lepton mass. We neglect target (nucleon) 
mass corrections which are important at larger values of Bjorken-x [19] where heavy quark 
contributions are of minor importance. 

We consider DIS of the virtual Boson B* on the quark Qi with mass mi producing 
the quark Q2 with mass At order 0(«g) this proceeds through the parton model 
diagram in Fig. 3.1 (a). 

Finite mass corrections to the massless parton model expressions are taken into account 
by adopting the Ansatz given in Eq. (17) of [19] 

W^^ = I fQii^,f^')^n,t=iP^ ■ (3.2) 
W^'^ is the usual hadronic tensor and lj'^'^ is its partonic analogue. Here as in the following 
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a hat on partonic quantities refers to unsubtracted amplitudes, i.e. expressions which still 
contain mass singularities in the massless limit, and are the light-cone momentum 
components of the incident quark Qi and the nucleon, respectively. Generally the '+' 
light-cone component of a vector v is given by v+ = {v^ + v^)/\f2. 

Contracting the convolution in Eq. (3.2) with the projectors in Appendix A.l gives 
the individual hadronic structure functions Fi=i,2,3- In leading order (LO) the latter are 
given by [19] 



F^'''\x,Q^) = 2R,Q,ix,Q') (3.3) 



with 

E±± ^Q^ ±ml±ml . (3.4) 

In Eq. (3.3) we use the shorthand A = A[m^, m|, —Q^] , where the usual triangle function 
is defined by 



A[a, b, c] = ^/a'^ + b'^ + c'^- 2{ab + bc + ca) . (3.5) 

The vector (V) and axial vector (A) couplings of the Q2ln{V — ^75)Qi quark current 
enter via the following combinations: 

S± = VV'±AA' 

= (VA'±V'A)/2 (3.6) 

where V,A = V, A' in the case of pure B scattering and V,A^ V, A' in the case of 
B, B' interference (e.g. 7, interference in the standard model). The scaling variable x 
generalizes the usual Bjorken-x in the presence of parton masses and is given by [19]: 

X = ^(S,. +A) . (3.7) 

The mass dependent structure functions in Eq. (3.3) motivate the following definitions 

_ 2A p 

^2 = ^T.F2 \ ^ Qi{x,Q') + 0{al) (3.8) 
such that J^i — Tj, i,j — 1, 2, 3, will be finite of 0{as) in the limit mi,2 — > 0. 
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3.2.2 DIS on a Massive Quark at 0{al] 



At 0{a\) contributions from real gluon emission [Fig. 3.1 (b)] and virtual corrections 
[Fig. 3.1 (c)] have to be added to the 0{a^g) results of Section 3.2.1. The vertex correction 
with general masses and couplings [Fig. 3.1 (c)] does to our knowledge not exist in the 
literature and is presented in some detail in Appendix A. 2. The final result (virtual+real) 
can be cast into the following form: 



•^1^1,2,3 {^iQ ) A* ) 



(3.9) 



= Qi(x,At') + 



27r 



with ^' = I and 

Hl{i\m^,m2) = CF 



{S, + Vi)6{l-i') + 



Qi 



X 2 



Hf (^',1711, 1712) 



rrin 



(3.10) 



(1 - e)+ 8s 

where s = (pi + q)'^ and the Si, Vi, Ni and are given in Appendix A. 3. The factorization 
scale fi^ will be taken equal to the rcnormalization scale throughout. The '+' distribution 
in Eq. (3.10) is a remnant of the cancellation of the soft divergences from the real and 
virtual contributions. It is defined as usual: 



/' de m') [giai ^ f di' [m - /(i)] g{o 

Jo Jo 



(3.11) 



As indicated by the hat on H^, the full massive convolution in Eq. (3.9) still contains the 
mass singularity arising from quasi-collinear gluon emission from the initial state quark 
leg. The latter has to be removed by subtraction in such a way that in the asymptotic 
limit Q"^ ^ 00 the well known massless MS expressions are recovered. The MS limit is 
mandatory since all modern parton distributions -and therefore all available heavy quark 
densities- are defined in this particular scheme (or in the DIS scheme [48], which can be 
straightforwardly derived from MS). The correct subtraction term can be obtained from 
the following limit 



lim 

mi—*0 



di'. 



1 + e 



/2 



1-e 



ln^-l-21n(l-e'; 



+ 

\g2 



(3.12) 
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where A = Q"^ / {Q"^ + m\) , x/\ = x 



mi- 



-_o and the Hf'^^ can be found in [49, 50]. Obviously 
the MS subtraction term for a 'heavy quark inside a heavy quark' is given not only by the 
splitting function Pqq^ = Ci?[(l + — ^')]+ times the coUinear log ln{fi'^/m'l) but also 

comprises a constant term. Herein we agree with Eq. (3.15) in [51]^, where this was first 
pointed out in the framework of perturbative fragmentation functions for heavy quarks. 
We therefore define 

'/2 



+ 



e.(f..^) 

(3.13) 



such that 



lim 

mi— »0 



,MS, 



(3.14) 



where the superscript Q\ on '^^ refers to that part of the inclusive structure function 
Ti which is initiated by the heavy quark Qi, i.e. which is obtained from a convolution 
with the heavy quark parton density. Note that the limit in Eq. (3.12) guarantees that 
Eq. (3.14) is also fulfilled when m\ — m2 (e.g. NC leptoproduction of charm) since 



(3.15) 



where Cf are the standard massless coefficient functions in the MS scheme, e.g. in 
[30,48]. 



Compcirison to Existing NC and CC Results 

We have performed several cross checks of our results against well known calculations that 
exist in the literature [30, 48, 49, 52, 53]. The checks can be partly inferred from the above 
paragraph. Nevertheless we present here a systematic list for the reader's convenience 
and to point out discrepancies of our calculation with [53] . 

In the charged current case V — A — 1 our results in Eq. (3.9) reduce in the limit 

mi ^ to the corresponding expressions in [52] , or in [49, 50] if the scheme dependent 

term represented by Eq. (3.13) is taken into account. The latter agrees with Eq. (3.15) 

^We also agree with the quark initiated coefficient functions in [52] where quark masses have been 
used as regulators. 
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in [51]. For mi 2 ^ we reproduce the well known MS coefficient functions, e.g., in [30, 
48,54]. The vertex correction in Appendix A. 2 is implicitly tested because it contributes 
to any of the final results. However, as an independent cross check the well known QED 
textbook result can be reproduced for mi = m2, A = 0. 

Initial state parton mass effects in NC DIS at 0{a\) have been ffist considered in [53] 
within the scenario [55] of intrinsic nonperturbative cc pairs stemming from fluctuations of 
the nucleon Fock space wavefunction. Although we do not consider such a scenario here 
we note that our results could be easily transferred to corresponding applications [56]. 
The main difference would be an inclusion of kinematical target mass effects which are 
important at larger x [19] where a possible nonperturbative charm component is expected 
[55] to reside. We list a detailed comparison of our calculation to the one in [53] in 
Appendix A. 4. Since our calculation does not fully agree with [53] for reasons which we 
are completely able to trace back and given the amount of successful independent tests 
of our results we regard the disagreement with [53] as a clear evidence that the results in 
[53] should be updated by our calculation. 

3.2.3 Gluon Fusion Contributions at 0{al) 

The gluon fusion contributions to heavy quark structure functions, depicted in Fig. 3.2, 
{B*g — s> Q1Q2) are known for a long time [33-35] and have been reinterpreted in [20] within 
the helicity basis for structure functions. Here we only briefly recall the corresponding for- 




Figure 3.2: Feynman diagrams for the production of a massive quark-antiquark pair via 
boson-gluon fusion. 

mulae in the tensor basis for completeness. The GF component of DI structure functions 
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is given by 

F^'ix^Q') = £ |^eW,/^^)/2(|,Q^) (3.16) 

where ax = [1 + (mi + m-i)^ jQ^\x and the can be found for general masses and cou- 
phngs in [33]. The corresponding Tf^ are obtained from the Ff^ by using the same 
normalization factors as in Eq. (3.8). Along the lines of [20] the GF contributions coexist 
with the QS contributions which are calculated from the heavy quark density, which is 
evolved via the massless RG equations in the MS scheme. As already pointed out in 
Section 3.2.2 the quasi-coUinear log of the fully massive GF term has to be subtracted 
since the corresponding mass singularities are resummed to all orders in the massless RG 
evolution. The subtraction term for the GF contribution is given by [20] 

Tt--'i.,Q\fi) = E ^ £ f <'«') S (f .''^) . (3-17) 

where Pf^{^) = 1/2 [^'^ + (1 - ^')^]. Note that Eq. (3.17) as well as Eq. (3.13) are 
defined relative to the Ti in Eq. (3.8) and not with respect to the experimental structure 
functions Fj. The sum in Eq. (3.17) runs over the indices of the quarks for which 
the quasi-coUinear logs are resummed by massless evolution of a heavy quark density, i.e. 
A; = 1, A; = 2 or A; = 1,2. 

3.2.4 ACOT Structure Functions at C>(aJ) 

As already mentioned in the introduction to this chapter it is not quite clear how the per- 
turbation series should be arranged for massive quarks, i.e. whether the counting is simply 
in powers of ccs as for light quarks or whether an intrinsic heavy quark density carries an 
extra power of due to its prehistory as an extrinsic particle produced by pure GF. We 
are here interested in the QS^^ component of heavy quark structure functions. Usually 
the latter is neglected in the ACOT formalism since it is assumed to be suppressed by 
one order of as with respect to the GF contribution as just explained above. We have, 
however, demonstrated in Chapter 2 within MS that this naive expectation is quantita- 
tively not supported in the special case of semi-inclusive production of charm (dimuon 
events) in CC DIS. We therefore want to investigate the numerical relevance of the QS^^^ 
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contribution to general heavy quark structure functions. In this chapter we present results 
for the fully inclusive case, relevant for inclusive analyses and fits to inclusive data. We 
postpone experimentally more relevant semi-inclusive (2;-dependent) results to Chapter 
4. Our results at full 0(q;J) will be given by 

^(1) ^ _^ pGF _ pSUBg _ pSUB, /g 

with Ff^^'-°^'\ F^'^, Ff^^" and Ff^^' given in Eqs. (3.9), (3.16), (3.13), and (3.17), 
respectively. Furthermore, we will also consider a perturbative expression for Fi which is 
constructed along the expectations of the original formulation of the ACOT scheme, i.e. 
QS^^^ is neglected and therefore F^^^'' need not be introduced 

^(0)+GF-SUBg ^ pQS^O) ^ ^Qp _ pSUBg 



3.3 Results for NC and CC Structure Functions 

In this section we present results which clarify the numerical relevance of QS^^^ contribu- 
tions to inclusive heavy quark structure functions in the ACOT scheme. We will restrict 
ourselves to NC and CC production of charm since bottom contributions are insignificant 
to present DI data. Our canonical parton distributions for the NC case will be CTEQ4M 
[42] (Figs. 3.5 and 3.6 below), which include 'massless heavy partons' above the scale 
= ml. Figures 3.3 and 3.4, however, have been obtained from the older GRV92 [57] 
distributions. The newer GRV94 [41] parametrizations do not include a resummed charm 
density since they are constructed exclusively along FOPT. GRV94 is employed in the CC 
case. Section 3.3.2. The radiative strange sea of GRV94 seems to be closest to presently 
available CC charm production data [31]. Furthermore, the low input scale of GRV94 al- 
lows for a wide range of variation of the factorization scale around the presently relevant 
experimental scales, which are lower for CC DIS than for NC DIS. Qualitatively all our 
results do not depend on the specific set of parton distributions chosen. 



3.3.1 NC Structure Functions 

For our qualitative analysis we are only considering photon exchange and we neglect the 
Z^. The relevant formulae are all given in Section 3.2 with the following identifications: 



Qi,2 c 
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Figure 3.3: x-dependence of the subtracted QS^^^ contribution to the NC charm structure 
function (solid line). Q'^ = fj? = 10 GeV^ is fixed. For comparison the MS analogue in 
Eq. (3.20) is shown (dashed line). The GRV92 parton distributions have been used. 



V 



mi,2 ^ = 1.6 (1.5) GeV for CTEQ4 (GRV92) 
■ V',A^A ^,0 



and wc use = if not otherwise noted. We consider contributions from charmed 
quarks and anti-quarks which are inseparably mixed by the GF contribution. This means 
that in Eq. (3.17) the sum runs over k = 1,2 and the relevant expressions of Section 3.2.1 
and 3.2.2 have to be doubled [since c{x,ii^) = c{x,ii^)]. 

First we investigate the importance of finite mass corrections to the limit in Eq. (3.14). 



In Fig. 3.3 the difference 

^(c+c)(i),MS ^ 4 



-K 



SUB„ 



can be compared to its MS analogue which is 



X 



27r 



{x,Q') 



(3.20) 



where denotes the usual (massless) convolution. From Fig. 3.3 it is obvious that the 
relative difference between ACOT and MS depends crucially on x. It can be large and 
only slowly convergent to the asymptotic MS limit as can be inferred from Fig. 3.4. Note 
that the solid curves in Figs. 3.3, 3.4 are extremely sensitive to the precise definition of 
the subtraction term in Eq. (3.13), e.g. changing x ^ ^ -which also removes the collinear 
singularity in the high limit- can change the ACOT result by about a factor of 5 
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X 10 




10 lo' lo' Q'[GeV'] 

Figure 3.4: The same as Fig. 3.3 but varying Q^{— jj?') for fixed x. 

around ^ 5 GeV^^ ^his is an example of the ambiguities in defining a variable flavor 
number scheme which have been formulated in a systematic manner in [22] . 

The relative difference between the subtracted QS^^'' contribution calculated along 
ACOT and the corresponding MS contribution in Eq. (3.20) appears, however, phe- 
nomenologically irrelevant if one considers the significance of these contributions to the 
total charm structure function in Fig. 3.5. The complete 0{al) result (solid hne) is shown 
over a wide range of Q'^ together with its individual contributions from Eq. (3.18). It can 
be clearly seen that the full massive QS^^'' contribution is almost completely cancelled by 
the subtraction term SU Bg (Indeed the curves for QS^^^ and SU Bq are hardly distinguish- 
able on the scale of Fig. 3.5). The subtracted quark correction is numerically neghgible 
and turns out to be indeed suppressed compared to the gluon initiated contribution, which 
is also shown in Fig. 3.5. Note, however, that the quark initiated corrections are not unim- 
portant because they are intrinsically small. Rather the large massive contribution QS^^^ 
is perfectly cancelled by the subtraction term SUBq provided that jj? — is chosen. 
This is not necessarily the case for different choices of \j? as we will now demonstrate. 

In Fig. 3.6 wc show the dependence of the complete structure function and its com- 
ponents on the arbitrary factorization scale /i^. Apart from the canonical choice //^ = 
^The subtracted gluon contribution GF changes by about a factor of 2 under the same replacement. 
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Figure 3.5: The complete 0{a\) neutral current structure function F2 and all individual 
contributions over a wide range of Q"^, calculated fronn the CTEQ4M distributions. Details of 
the distinct contributions are given in the text. 

(which was used for all preceding figures) also different scales have been proposed [58, 59] 
like the maximum transverse momentum of the outgoing heavy quark which is approx- 
imately given by (p™"^)^ ~ {1/x — 1) (5^/4. For low values of 2;, where heavy quark 
structure functions are most important, the scale (p™"^)^ ^ . The effect of choosing 
a which differs much from can be easiest understood for the massless coefficient 
functions Cf^'^^ which contain an unresummed \n{Q^ / ji'^) . The latter is of course absent 
for /i^ = but becomes numerically increasingly important, the more /i^ deviates from 
Q^. This logarithmic contribution cannot be neglected since it is the unresummed part 
of the collinear divergence which is necessary to define the scale dependence of the charm 
density. This expectation is confirmed by Fig. 3.6. For larger values of the subtracted 
QS^^^ contribution is indeed still suppressed relative to the subtracted GF contribution. 
Nevertheless, its contribution to the total structure function becomes numerically signifi- 
cant and reaches the ~ 20 % level around (p^"^)^. Note that in this regime the involved 
formulae of Section 3.2.2 may be safely approximated by the much simpler convolution 
in Eq. (3.20) because they are completely dominated by the universal collinear logarithm 
and the finite differences ACOT — MS from Figs. 3.3 and 3.4 become immaterial. In prac- 
tice it is therefore always legitimate to approximate the ACOT results of Section 3.2.2 by 
their MS analogues because both are either numerically insubstantial or logarithmically 



3.3 Results for NC and CC Structure Functions 



25 




Figure 3.6: iJ? dependence of the complete 0{q\) NC structure function (CTEQ4M) in Eq. 
(3.18) (solid line) and of the structure function in Eq. (3.19) (dot-dashed line) where the 
subtracted QS'^^^ contribution is neglected. Also shown are the different subtracted 0{a\) 
contributions GF and QS^^\ 

dominated. Finally we confirm that the scale dependence of the full 0{a\) structure func- 
tion F^^ in Eq. (3.18) is larger than the scale dependence of p^^)+^^-^^^9 (3.19) 
which was already pointed out in [60]. Nevertheless, the subtracted QS^^^ contribution 
should be respected for theoretical reasons whenever ln((5^//i^) <^ 1. 

3.3.2 CC Structure Functions 

Charm production in CC DIS is induced by an s — > c transition at the 1^-Boson vertex. 
The strange quark is not really a heavy quark in the sense of the ACOT formalism, i.e., 

the production of strange quarks cannot be calculated reliably at any scale using FOPT 
because the strange quark mass is too small. It is nevertheless reasonable to take into 
account possible finite effects into perturbativc calculations using ACOT since the 
subtraction terms remove all long distance physics from the coefficient functions. Indeed 
the ACOT formalism has been used for an experimental analysis of CC charm production 
in order to extract the strange sea density of the nucleon [17]. Along the assumptions of 
ACOT (55**^^^ contributions have not been taken into account. This procedure is obviously 
questionable and has been shown not to be justified within the MS scheme in Chapter 2. 
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With our results in Section 3.2.2 we can investigate the importance of quark initiated 
0{a\) corrections within the ACOT scheme for inclusive CC DIS. As already mentioned 
above, results for the experimentally more important case of semi-inclusive (z-dependent) 
DIS will be presented in Chapter 4 of this thesis. In the following we only introduce 
subtraction terms for coUinear divergencies correlated with the strange mass and treat 
all logarithms of the charm mass along FOPT. We do so for two reasons, one theoretical 
and one experimental: First, at present experimental scales of CC charm production 
ln((5^/w2^) terms can be safely treated along FOPT and no introduction of an a priori 
unknown charm density is necessary. Second, the introduction of a subtraction term for 
the mass singularity of the charm quark would simultaneously require the inclusion of the 
c — s> s QS'-transition at the W^-vertex with no spectator-like c-quark as in GF. This 
contribution must, however, be absent when experiments tag on charm in the final state. 
CC DIS on massive charm quarks without final-state charm tagging has been studied in 
[61]. 

The numerics of this section can be obtained by the formulae of Section 3.2 with the 
following identifications: 

Qi — s , Q2 — c 
m2^m^ = 1.5 GeV (GRV94) 
V^V',A^A' 1,1 

and the strange mass mi = will be varied in order to show its effect on the structure 
function F2. 

In Fig. 3.7 we show the structure function F| and its individual contributions for two 
experimental values of x and [16] under variation of the factorization scale //^. Like in 
the NC case we show the complete 0{a\) result as well as p(^)+'^^-^^^s ^j^Q^e QS^^^ has 
been neglected. The thick curves in Fig. 3.7 (a) have been obtained with a regularizing 
strange mass of 10 MeV. They are numerically indistinguishable from the rUs — MS 
results along the lines of [49] . For the thin curves a larger strange mass of 500 MeV has 
been assumed as an upper limit. Finite mass effects can therefore be inferred from the 
difference between the thin and the thick curves. Obviously they are very small for all 
contributions and can be safely neglected. For the higher Q'^ value of Fig. 3.7 (b) they 
would be completely invisible, so we only show the =10 MeV results (=MS). Since the 
finite mass corrections within the ACOT scheme turn out to be neghgible as compared 
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Figure 3.7: The charm production contribution to the charged current structure function 
F2 for a wide range of the factorization scale using GRV94. The curves are as for the 
neutral current case in Fig. 3.6. In Fig. 3.7 (a) the thicker curves have been obtained with 
a (purely regularizing) strange mass of 10 MeV which according to Eq. (3.14) (and to the 
analogous limit for the subtracted GF term [20]) numerically reproduces MS. For the thinner 
curves a strange mass of 500 MeV has been assumed. In Fig. 3.7 (b) all curves correspond 
to rus = 10 MeV (= MS). 
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to massless MS it is not surprising that we confirm the findings of Chapter 2 concerning 
the importance of quark initiated corrections. They are -in the case of CC production of 
charm- not suppressed with respect to gluon initiated corrections for all reasonable values 
of the factorization scale. Only for small choices of //^ ~ + ra^ can the quark initiated 
correction be neglected. In this region of /x^ also gluon initiated corrections are moderate 
and Born approximation holds within ~ 10%. For reasons explained in Section 3.3.1 the 
absolute value of both corrections -gluon and quark initiated- become very significant 
when large factorization scales like p^"-^ are chosen. This can be inferred by looking at 
the region indicated by the arrow in Fig. 3.7 which marks the scale ^ = 2p™"'^ which was 
used in [17]. Analyses which use ACOT with a high factorization scale and neglect quark 
initiated corrections therefore undershoot the complete 0{al) result by the difference 
between the solid and the dot-dashed curve, which can be easily as large as ~ 20 %. For 
reasons explained in the introduction to this section we have used the radiative strange 
sea of GRV94. When larger strange seas like CTEQ4 are used the inclusion of the quark 
initiated contributions is even more important. 

3.4 Conclusions 

In this chapter we have calculated and analyzed DIS on massive quarks at 0{al) within the 
ACOT scheme for heavy quarks. For NC DIS this contribution differs significantly from its 
massless MS analogue for — Q"^. Both give, however, a very small contribution to the 
total charm structure function such that the large relative difference is phenomenologically 
immaterial. At higher values of the factorization scale jjF' ~ (pjp'^^y the contributions 
become significant and their relative difference vanishes. The QS^^^ contribution of Section 
3.2.2 can therefore be safely approximated by its much simpler MS analogue at any scale. 
For CC DIS quark initiated corrections should always be taken into account on the same 
level as gluon initiated corrections. Due to the smallness of the strange quark mass ACOT 
gives results which are almost identical to MS. 



Chapter 4 

Charm Fragmentation in Deep 
Inelastic Scattering 

This chapter generahzes the theoretical considerations of Chapter 3 to the semi-inclusive 
deep inelastic production of heavy flavored hadrons. We apply these results to the CC 
and NC production of charm where we investigate the universality of charm fragmentation 
functions and the charm production dynamics, respectively. The work in this chapter is 
based on [4]. 

4.1 Introduction 

In the preceeding chapter we have analyzed heavy quark initiated contributions to fully 
inclusive deep inelastic (DI) structure functions. Towards lower values of the Bjorken 
variable x, heavy (charm) quarks are produced in about 20% of the neutral current (NC) 
[11-13] and charged current (CC) [43] deep inelastic events in lepton-nucleon collisions. 
Therefore in this kinematical range, heavy quark events contribute an important compo- 
nent to the fully inclusive DI structure functions of the nucleon. However, due to accep- 
tance losses this component can usually not be measured directly by inclusively tagging on 
charm events and more differential observables have to be considered hke Ed [16, 17, 45], 
Pt or T) [11-13, 15] distributions, where Ed, Pt and rj are the energy, transverse momen- 
tum and pseudorapidity of the charmed hadron produced, i.e. mainly of a D*^*^ meson. In 
this chapter we consider Ed spectra represented by the usual scaling variable z defined 
below. Within DIS the charmed hadron energy spectrum is the distribution which is most 
sensitive to the charm fragmentation process and may give complementary information 
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to one hadron inclusive e'^e~ annihilation which is usually chosen to define fragmentation 
functions (FFs) [29, 62]. A well understanding of charm fragmentation is essential for any 
charm observable, e.g. the normalization of pr and rj distributions in photoproduction is 
substantially influenced by the hardness of the FF [63,64]. The z distribution of charm 
fragments is directly measured in CC neutrinoproduction [16, 17,45] and may give insight 
into details of the charm production dynamics in NC electroproduction. It has, e.g., been 
shown in [11,12] that the energy spectrum of D(*^-mesons produced at the ep collider 
HERA may be able to discriminate between intrinsic and extrinsic production of charm 
quarks. 

Intrinsic heavy quark densities may arise due to a nonperturbative component of the 
nucleon wave function [55] or due to a pcrturbative resummation [20,22-24] of large 
quasi-coUinear logs [ln((5^/m^); Q and m being the virtuality of the mediated gauge 
boson and the heavy quark mass, respectively] arising at any order in fixed order extrinsic 
production (or fixed order perturbation theory: FOPT). Here we will only consider the 
latter possibility for inducing an intrinsic charm density c{x, Q^) which is concentrated 
at small x and we will ignore nonperturbative components which are expected to be 
located at large x [55] . Technically the resummation of large perturbative logs proceeds 
through calculating the boundary conditions for a transformation of the factorization 
scheme [23, 47, 51], which is switched from n/ to n/ + l active, massless flavors, canonically 
at Q"^ — m^. For fully inclusive DIS the Kinoshita-Lee-Nauenberg theorem conflnes all 
quasi-coUinear logs to the initial state such that they may be absorbed into c{x, Q^). For 
semi-inclusive DIS (SI DIS) also flnal state coUinearities arise which are resummed in 
perturbative fragmentation functions D^{z, Q"^) (parton i decaying into charm quark c) 
along the fines of [51, 63]. The scale dependence of c{x, Q"^) and D1{z, Q"^) is governed by 
massless renormalization group (RG) evolution. 

Besides this zero mass variable flavor number scheme^ where mass effects arc only taken 
care of by the boundary conditions for c{x^Q'^) and D1{z^Q'^)^ variable flavor number 
schemes have been formulated [20,22-24], which aim at resumming the quasi-collinear 
logs as outlined above while also keeping power suppressed terms of 0[{m'^ /Q'^Y] in the 
perturbative coefficient functions. Our reference scheme for this type of schemes will be 
the one developed by Aivazis, Collins, Olness and Tung (ACOT) [19-21]. In the ACOT 
scheme full dependence on the heavy quark mass is kept in graphs containing heavy 
quark lines. This gives rise to the above mentioned quasi-collinear logs and to the power 
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suppressed terms. While the latter are regarded as mass corrections to the massless, 
dimensionally regularized, standard coefficient functions (e.g. in the MS scheme), the 
former are removed numerically by subtraction terms, which are obtained from the small 
mass limit of the massive coefficient functions. 

The outline of this chapter will be the following: In Section 4.2 we will shortly overview 
the relevant formulae for SI DIS for general masses and couplings including quark scat- 
tering (QS) and boson gluon fusion (GF) contributions up to 0{al). We will thereby 
present our ACOT based calculation for the QS^^^^ component of SI structure functions. 
In Section 4.3 we will analyze the charm fragmentation function in CC and NC DIS. In 
Section 4.4 we draw our conclusions and some uncomfortably long formulae are relegated 
to Appendix B. 

4.2 Semi— Inclusive Heavy Quark Structure Func- 
tions 

This section presents the relevant formulae for one (heavy flavored) hadron inclusive DIS 
structure functions. The contributions from scattering events on massive quarks are given 
up to 0{al) in Section 4.2.1 and GF contributions are briefly recalled in Section 4.2.2. 
Section 4.2.3 presents all subtraction terms which render the structure functions infrared 
safe and includes a discussion of these terms. 

4.2.1 Scattering on Massive Quarks 

We consider DIS of the virtual Boson B* with momentum q on the quark Qi with mass 
mi and momentum pi producing the quark Q2 with mass m2 and momentum p2. The 
latter fragments into a heavy quark flavored hadron Hq,^, e.g. a \Q2(li) meson, qi being any 
light quark. Phenomenologically most prominent are of course charm quarks fragmenting 
into -mesons which are the lightest heavy flavored hadrons. 

We will strictly take over the formulae and notations of our inclusive analysis in 
Chapter 3 whenever possible and merely extend them for SI DIS considered here. In 
particular we take over the definition of the structure functions J^i given in terms of the 
usual experimental structure functions Fj in Eq. (3.8). 

^Bracketed upper indices count powers of a^. 
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Since we want to investigate the energy spectrum of charm fragments, we introduce 
the Lorentz-invariant z = Phq^ ■ Pn/q ■ Pn which reduces to the energy Ehq^ scaled to 
its maximal value 1/ = qq in the target rest frame. Therefore in contrast to Eq. (A. 4) we 
do not integrate the tensor u;^^'^ over the full partonic phase space but keep it differential 
in the corresponding partonic variable z' = p2 ■ Pi/q ■ Pi or the mass corrected variable z 
which is defined below. In order to obtain hadronic observables we have to extend the 
Ansatz of Eq. (17) in [19] such that it includes a nonperturbative hadronization function 
Dq^. In the limit of vanishing masses the massless parton model expressions have to be 
recovered. Our Ansatz will be 



where /i is the factorization scale, z = z' j z'^^q with z'i^q = and = (v^ + v^)/ \/2 

for a general vector v. W^'^ is the usual hadronic tensor and Cj^^ is its partonic analogue. 
Eq. (4.1) defines the fragmentation function Dq^ to be a multiplicative factor multiplying 
inclusive structure functions at LO/Born accuracy, i.e. 



where the 2 sl'"^! Q^) ^^^^ defined and given in (3.8). The scale dependence of Dq^ is 
bracketed here and in the following because it is optional; a more detailed discussion on 
this point will be given at the end of Section 4.2.3. We do not construct our Ansatz in Eq. 
(4.1) for the convolution of the fragmentation function along light front components for 
the outgoing particles which would only be Lorentz-invariant for boosts along a specified 
axis. Since the final state of DIS is spread over the entire solid angle it has no preferred 
axis as defined for the initial state by the beam direction of collider experiments. Note 
that Eq. (4.1) is in agreement with usual factorized expressions for massless initial state 
quanta, as considered, e.g., in Chapter 2, since there mi = such that z'j^q — 1. 

Up to 0{a\) the hadronic structure functions for scattering on a heavy quark read 




(4.1) 



(4.2) 



'^1=1,2,3 




(4.3) 



with [19] 



X 



+ A) 



(4.4) 
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As throughout this thesis we set the renormahzation scale equal to the factorization scale 
In Eq. (4.3). The kinematical boundaries of the phase space in Eq. (4.3) are introduced 
by the theta function cut Qq. In the massless limit 6q — > 1. The precise arguments of Qq 
are set by the kinematical requirement 



with 



^ max 



±A[s, ml -Q^] {s-ml) + (Q^ + mf + g) (g + mj) 
2s{s -mj + Q2) 



where s — {pi + q)^- Note that Eq. (4.5) also poses an implicit constraint on ^' via 

z' < z' 

mm ^max' 

The Hf for nonzero masses are obtained in exactly the same way as outlined for 
fully inclusive structure functions in Appendix A if the partonic phase space is not fully 
integrated over. They are given by 

HHli',z') = Cf [ (S. + K) S{1 - e) S{1 - i) 

where the normalization factors Ni are given in Appendix A. 3 and the Mandelstam vari- 
ables read 



s,{0 = iPi + lf 



m 



= 1-^[(A-E,_)e + A + E,_] 

= [5i(0 + s,-](^'-4), 

where z'r. — — is the would-be pole of the ^-channel propagator and we use A' = 

si + 

A[s,m\, —Q"^]. Si and Vi are the soft real and virtual contribution to iff, respectively. 
They can be found in Appendix A. 3 whereas the f^{si,ii) are listed in Appendix B. In 
the massless limit the H!^ reduce to the MS coefficient functions in [29,30] up to some 
divergent subtraction terms which we will specify in Section 4.2.3. 
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4.2.2 Gluon Fusion Contributions at 0{al) 

We have already obtained the semi-inclusive coefficient functions for GF production of 
massive quarks for general masses and couplings in Chapter 2. They are given by 

F,^^{x,z,Q\i,') = ^ x'g{x',t,') M^,CQ') Dq,C-,[i^']) 6, (4.7) 

Jax J z S ^ S 

with the fi from Eq. (2.2) and ax — [1 + (mi + m2)^/(5^]x. The Qg cut guarantees that 
Cmin < C < Cmax whcrc Cmin,max are givcu in Eq. (2.3). Similarly to Eq. (4.5) Cmin < Cmax 
may also constrain the phase space available for the x' integration. 



4.2.3 Subtraction Terms 



It requires three MS subtraction terms to render the double convolutions in Eqs. (4.3), 
(4.7) infrared safe: 

/2 



i±i:fin.^' 
1-e 



ml 



l-21n(l-0 



+ 



(4.8) 



1 + /i^ 



X Dq, (|:,)"') (3i(x,/^') , 



(4.10) 



where Pig\0 = 1/2 + (1 - C')^]- Note that SUBg in Eq. (4.9) differs slightly from 
Eq. (2.6) in Chapter 2 because we are allowing for a nonzero initial state parton mass mi 
here which we did not in Chapter 2. 

The subtraction terms define the running of the initial state quark density {SUBg, 
SUBg) and the final state fragmentation function [SUBd] in the massless limit. They 
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remove coUinear logarithms and scheme defining finite terms from the convolutions in 
Eqs. (4.3), (4.7) and they are constructed such that the massless MS results of [29, 30] are 
recovered in the limit 



lim 

"11,2— >0 



(4.11) 



where SUBq and SUBn regularize whereas SUBg regularizes Contrary 

to the fully inclusive SUBg term in Eq. (3.17) there is no need to include an additional 
~ Indji^/ml) subtraction in Eq. (4.9) because the ('^ it-channel singularity of the massless 
hmit of GF^^) is located at C = and is outside the integration volume of Eq. (4.7). If only 
initial state subtractions, i.e. SUBg and SUBg, are considered and final state subtractions, 
i.e. SUBd, are not performed one reproduces, in the limit mi — > the results in [31,50] 
for producing a heavy quark from a light quark. Note that in this case the fragmentation 
function Dq^ should be taken scale-iridependent, say of the Peterson form [37] . In the limit 
where also the final state quark mass m2 approaches zero and the final state subtraction 
term SUB^ is subtracted from the results in [31, 50]^ the massless quark results in [29, 30] 
are obtained and a running of Dq^ is induced via a RG resummation of final state coUinear 
logs as formulated for one hadron inclusive e'^e~ annihilation in [51, 63]. 

Apart from removing the long distance physics from the coefficient functions the sub- 
traction terms set the boundary conditions for the intrinsic heavy quark density Qi [47] 
and the perturbative part of the heavy quark fragmentation function Dq^ [51]: 



Qi{x,Ql) 



271 



m 



27r 



Up 



1 

z' 



19 
/2 



o Qo 



(4.12) 



1 + z 
1-z' 



ln^-l-21n(l-^') 



mi 



+ 



(4.13) 



where Qo, Qo are the transition scales at which the factorization scheme is switched 
from n/ton/-|-l, n/ + 2 active flavors, respectively (assuming here for simplicity that 
mi < m2; a generalization to mi > m2 is straightforward). For general Qo also the gluon 
density and ag undergo a scheme transformation. Canonically Qo is set equal to the heavy 
quark mass mi which guarantees [47] up to two loops a continuous evolution of ag and 
2See Equation (B.ll) in [50]. 
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the light parton densities across Qq. All available heavy quark densities are generated 
using Qq — TOi and we will therefore follow this choice here although a variation of Qo 
might substantially influence the heavy quark results even far above the threshold [18] 
as was found in [65]. Note that at three loops a continuous evolution across Qq can no 
longer be achieved, neither for the parton distributions [23] nor for ag [66]. Analogously 
to Qo = mi we use Qq — m2 throughout. In Eq. (4.13) we have made the distinction 
between the scale dependent FF Dq.^{z, Ql) and the scale independent FF Dq.^{z) explicit. 
Following the terminology in [51,63] the latter corresponds to the nonperturbative part 
of the former and describes the hadronization process at the end of the parton shower 
which is described perturbatively by the massless RG evolution. Alternatively, Dq^{z) 
corresponds to a scale-independent FF within FOPT where no collinear resummations 
are performed. These two points of view may induce a scheme dependence if Dq^{z) is 
fitted to data. In principle, the massless evolution equations generate nonzero FFs also 
for light partons to decay into heavy fiavored hadrons. These light -^heavy contributions 
are important at LEP energies [63] but can be safely neglected at the scales considered 
here^ and we will assume D^^q^^z, fi^) = throughout. 

4.2.4 SI Structure Functions at 0{al) 

In the next section we will consider three types of 0{al) VFNS structure functions. The 
first two are constructed at full 0{al) 

^Q5(o+i)+GF _ pSUBg+SUBg+[SUBo] j^^^ 

where the inclusion or omission of the bracketed SUBd term corresponds to a running 
or scale-independent fragmentation function, respectively, as discussed in the previous 
section. It is somewhat unclear whether QS^^^ contributions (and the corresponding 
subtractions) should be considered on the same perturbative level as GF^^\ see the in- 
troduction to Chapter 3 for a more detailed discussion on that point. In the original 
formulation of the ACOT scheme [20] QS'^^^ contributions are neglected at the level we 
are considering here and we therefore do also consider this option via the partial 0{a\) 
structure function 

pQSm+GF _ pSUB, _ 



^We could therefore, in principle, restrict the evolution of the charm FF to the nonsinglet sector. 
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For obtaining the numerical results of the next section the general formulae of this section 
have to be adjusted by choosing masses and couplings according to the relevant NC and 
CC values as in Section 3.3.1 and 3.3.2. For the CC case where Qi should be identified 
with strange quarks the boundary condition in Eq. (4.12) is inadequate since ~ ^qcd 
is below the perturbative regime of QCD. We will have recourse to standard strange seas 
from the literature [41,42] instead. 

4.3 The Charm Fragmentation Function in SI DIS 

We will investigate the charm fragmentation function in CC and in NC SI DIS. In CC 
DIS the charm production mechanism is undebated since charm is dominantly produced 
by scattering on light strange quanta. Our reasoning will therefore be that Dc is directly 
accessible in CC DIS at relatively low spacelike momentum transfer. An extracted Dc can 
then be apphed to NC DIS, where it might give insight into the details of the production 
dynamics [11]. Also a test of the universahty [67] of the charm FF measured in CC DIS 
and e'^e~ annihilation [63] would be an important issue directly related to the factorization 
theorems [68] of perturbative QCD (pQCD). 

All Ec parameters discussed below refer to a Peterson type [37] functional form given 
in Eq. (2.8) 

4.3.1 CC DIS 

In CC DIS at fixed target energies one does not expect to gain much insight into the charm 
production process since charm is dominantly produced in scattering events on strange 
quarks^ in the nucleon, thereby permitting an experimental determination of the strange 
quark content of the nucleon [16,17,45,46]. On the other hand the well understood 
production mechanism makes a direct determination of the charm fragmentation function 
feasible by measuring the energy spectrum of final state charm fragments. This is obvious 
in leading order accuracy where^ 

daLO oc s{x)Dc{z) (4.16) 

^Wc assume a vanishing Cabibbo angle. Our results remain, however, unehanged if the CKM sup- 
pressed c background is included. 

^We will suppress some obvious scale-dependences in the following formulae and in their discussion. 
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is directly proportional to D^. More precisely, it is not z but the closely related energy of 
the c ^ /lu decay muon which can be observed in iron detectors [16, 17, 45]. The smear- 
ing effects of the decay comphcates the determination of the precise shape of Dc but only 
weakly influences an extraction of (z) [45] which is valuable information if physically mo- 
tivated one-parametric Ansatze [37, 69] for D^. are assumed. At NLO the production cross 
section is no longer of the simple factorized form of Eq. (4.16) and double convolutions 
(symbol <S) below) of the form of Eqs. (4.3), (4.7) have to be considered. However, to a 
reasonable approximation 



holds also at NLO accuracy within experimental errors and for the limited kinematical 
range of present data on neutrinoproduction of charm. In Eq. (4.17) the approximate 
multiplicative factor V absorbs the precise i^-factor K{x, z, Q^) obtained from a full 
NLO QCD calculation. V is not a simple universal fragmentation function but a nontrivial 
process-dependent functional which is, however, mainly sensitive on Dc and shows little 
sensitivity on the exact parton distributions considered. The occurrence of x, and z 
in Eq. (4.17) as indices and as a functional argument, respectively, reflects the fact that 
the dependence on x and is much weaker than is on z. Eq. (4.17) tells us that s{x) 
fixes the normalization of da once K is known. On the other hand K (or V) can be 
computed from Dc with little sensitivity on s{x), such that s{x) and Dc{z) decouple in 
the production dynamics and can be simultaneously extracted. This point can be clearly 
inferred from Fig. 4.1 where it is shown that the wide spread of CC charm production 
predictions which were obtained in [31] using GRV94 [41] and CTEQ4 [42] strange seas 
can be brought into good agreement by a mere change of the normalization given by the 
ratio sgrv{x) / ^cteqa{x) ■ The remaining difference is not within present experimental 
accuracy which can be inferred from the shaded band representing a parametrization 
[31] of CCFR data [16]. High statistics neutrino data therefore offer an ideal scenario 
to measure Dc complementary to an extraction from LEP data on e~^e~ DX [63, 
70,71]. This has been first noted in [67] where also a successful test of the universality 



da NLO — 



oc 



{[s (g) d&s + fif <8) dag] (g) Dc) (x, z, Q^) 
daio K{x, z, Q"^) 
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Figure 4.1: The charm production cross section obtained in [31] for two x, points in the 
CCFR [16] kinematical regime. Up to a constant of normalization which was conveniently 
chosen in Eq. (4) of [31] Sg/f represents the triple-differential cross section d^a / dxdydz 
where x and y are standard and z = po ■ Pn/q ■ Pn- Shown are the predictions using 
GRV94 [41] (solid) and CTEQ4 [42] (dashed) partons and a curve (dot-dashed) where the 
normalization of the CTEQ4 prediction is changed by multiplying with the ratio of the strange 
seas scRvix) / scteq{x) ■ Foi' curves a scale-independent Peterson FF with €c = 0.06 has 
been used. 



of the charm FF has been performed. With new data [16,17,45]^ at hand and with a 

sounder theoretical understanding of neutrinoproduction of charm it would be desirable 

to update the analysis in [67]. Nowadays one can in principle examine the possibility 

of a uniform renormalization group transformation from spacelike momenta near above 

the charm mass {uN —>■ DX) to timelike momenta at the peak (e^e^ DX). In 

[63, 70, 71] charm fragmentation functions extracted from LEP data have been tested 

against px and t] distributions measured in photoproduction at HERA. However, we 

believe that a comparison to z differential neutrinoproduction data is worthwhile beyond, 

^Data from NuTeV [72,73] are to be expected in the near future; /z^ events observed at NOMAD 
await further analysis [74]. 
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since the latter measure directly the fragmentation spectrum whereas pt and rj shapes 
are rather indirectly influenced by the precise hardness of the FF via their normalization 
[63, 64]. Unfortunately up to now no real production data are available but only strange 
sea extractions [16, 17, 45] resulting from an analysis of the former. We therefore strongly 
recommend that experiments produce real production data such that the above outlined 
program can be executed with rigor. Here we can only flnd an Sc parameter which lies 
in the correct ball park and examine a few points which will become relevant for an 
extraction of once data will become available. 

An outstanding question is the possible effect of a finite strange mass on the full semi- 
inclusive charm production cross section including 0{al) quark scattering contributions. 
By comparing the thick and the thin solid curve in Fig. 4.2 (a) it is clear that the effect 
of a finite rUg can be neglected even at low scales and for a maximally reasonable value of 
THg — 500 MeV. For the larger scale of Fig. 4.2 (b) the effect of choosing a finite would 
be completely invisible. A further question which might influence the extraction of a uni- 
versal FF from neutrinoproduction is the one of the scheme dependence in handling flnal 
state quasi-coUinear logarithms ln((5^/m^). If these are subtracted from the coefficient 
functions as discussed in Section 4.2.3, the subtraction defines a running of the charm 
FF which becomes scale dependent^ according to Eq. (4.13). In Fig. 4.2 we examine 
such resummation effects for CCFR kinematics [16]. We use the same Peterson FF with 
£c = 0.06 once for a fixed order calculation (solid lines) and once as the nonperturbative 
part on the right hand side of the entire c — > FF on the left hand side of Eq. (4.13) 
(dashed curves). We note that towards intermediate scales around Q'^ ~ 20 GeV^ one 
begins to see the softening effects of the resummation which are enhanced as compared 
to FOPT. However, as one would expect at these scales, the resummation effects are 
moderate and could be compensated by only a slight shift of the Ec parameter which is 
therefore, within experimental accuracy, insensitive to scheme transformations. We note 
that -as was already shown in [31]- according to Fig. 4.1 an Sc of around 0.06 which 
we took from an older analysis in [38] seems to reproduce the measured spectrum quite 
well. For (E^) = 80 GeV, (Q^) = 20 GeV^ a value Sc ~ 0.06 gives an average (z) ~ 0.6 
consistent with (z) = 0.68 ± 0.08 measured by CDHSW [45]. In [63]^ a distinctly harder 

''Evolving the charm FF we adopt for consistency the evolution parameters mc,b and A^^^ of the 
CTEQ4(M) [42] parton distribution functions and we use Qo = rric. 

similar Ec has been obtained in [70] in a related scheme. The Sc value in [71] has no connection to 
a massive calculation and cannot be compared to the values discussed here. 
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Figure 4.2: The same quantity as in Fig. 4.1. All predictions have been obtained for the 
CTEQ4 parton distributions. The solid curves result from a scale-independent Peterson FF 
with Ec = 0.06. In (a) the thin solid curve has been obtained using the formulae of Section 
4.2 and choosing a finite strange quark mass of rus — 500 MeV whereas the thicker curve 
corresponds to the asymptotic ^ limit (=MS). In (b) these two options would be 
completely indistinguishable and we only show the MS result. For the dot-dashed curves the 
final state quasi-collinear logarithm has been absorbed into a running of the charm FF. A 
Peterson FF with Sc — 0.06 has been used as the nonperturbative part of the input as given 
in Eq. (4.13). 



value of Sc — 0.02 was extracted from LEP data on e+e" — > D*X. If the latter fit is 
evolved down to fixed target scales it is -even within the limited experimental accuracy- 
incompatible with the CCFR neutrino data represented in Fig. 4.1. From Fig. 4.2 it is 

clear that the difference cannot be attributed to a scheme dependence of the Sc parameter 
which is too small to explain the discrepancy. It would of course be interesting to know 
how much the above mentioned smearing effect of the c ^ fiu decay might dilute the dis- 
crepancy. In any case, charm fragmentation at LEP has been measured by tagging on D*s 
whereas neutrinoproduction experiments observe mainly Ds through their semileptonic 
decay-channel (dimuon events). ARGUS [75] and CLEO [76] data at ^/s ~ lOGeV indeed 
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show [77] a harder energy distribution of D*s compared to Ds. It seems therefore to be 
possible within experimental accuracy to observe a nondegeneracy of the charm fragmen- 
tation functions into the lowest charmed pseudoscalar and vector mesons. We note that 
an Ec value around 0.06 which is in agreement with neutrino data on D-production is also 
compatible with the D energy spectrum measured at ARGUS where the evolution may 
be performed either via FOPT using expressions in [62] or via a RG transformation along 
the lines of [51, 63]. If forthcoming experimental analyses should confirm our findings the 
lower decade mc(~ 1 GeV) ARGUS(10 GeV) may be added to the evolution path 
ARGUS(10 GeV) ^ LEP(Mz) paved in [63] for the charm FF. 

4.3.2 NC DIS 

The fragmentation function of charm quarks observed in NC DIS is of special interest 
since it allows for directly investigating [11,12] the charm production mechanism which 
is a vividly discussed issue in pQCD phenomenology [18,23,78]. Whereas for intrinsic 
heavy quarks one expects to observe a Peterson-like hard spectrum attributed to the 
dominance of the leading order quark scattering contribution, one expects a much softer 
spectrum for extrinsic GF since the gluon radiates the cc pair towards lower energies 
(z) during the hard production process before the nonpcrturbativc hadronization takes 
place. Experimental analyses have been performed in [11,12] and the steep spectrum^ 
(best visible in Fig. 6 in [11]) of observed D-mcsons together with the missing of a hard 
component at larger z seem to give clear evidence for the dominance of extrinsic GF over 
QS which was excluded at the 5% level [11]. In a complete VFNS the QS'^'^^ component 
makes, however, just one part of the 0{al) structure functions in Eqs. (4.14), (4.15). 
Especially, there also exists a GF^^^ component, albeit with the leading log part of it 
subtracted. Since the subtraction term in Eq. (4.9) is proportional to Dc and therefore 
only removes a hard component from GF^^^ one expects the rise towards lower z to 
survive the subtraction. Furthermore a perturbative evolution of the charm fragmentation 
function Dc{z, Q^) might soften somewhat the hard QS term. 

These expectations can be quantitatively confirmed in Fig. 4.3 where we show for 

HERA kinematics [11] the total (solid hue) normahzed 0{al) production cross section 

for transverse virtual photons {Fl — 0) on protons. We also show the individual compo- 

^The variable xd considered in [11, 12] differs slightly from z in definition. The variables are, however, 
identical at the 2% level [79]. 



4-3 The Charm Fragmentation Function in SI DIS 



43 




Figure 4.3: The normalized charm production cross section da/dz {z = po ■ Pn/q ■ Pn) for 
HERA kinematics (^i = 314 GeV). For a comparison with HI data [11] the cross section 
has been integrated over 10 GeV^ < Q"^ < 100 GeV^ and 0.01 < y < 0.7. Following the 
experimental analysis [11] only contributions from transverse photons {F^ = 0) are considered. 
Shown is the total 0{al) result (solid line) and the individual contributions. Details to the 
calculation of the total result and the individual contributions are given in the text. The charm 
mass has been kept finite at the CTEQA value of nic = 1.6 GeV everywhere except for the 
thin dotted curve where the rric limit has been taken. 

nents contributing to it: The processes 7*5' cc (dot-dashed) and 7*c eg (dotted; inch 
virtual corrections) correspond to the GF^^^ — SUBg and QS^^^ — SUBq — SUBo terms, 
respectively, subtracted at fi = Q. They are not physically observable and only sensible 
if they are added to the QS^^^ Born term (dashed) as in Eqs. (4.14), (4.15). We have per- 
turbatively resummed all logarithms of the charm mass via massless evolution equations 
starting at the charm mass [Qq = Qo = itIc] and using the standard boundary conditions 
in Eqs. (4.12), (4.13) for Ec = 0.06. Finite charm mass effects on the subtracted QS'^^^ 
contribution can be inferred by comparing the thick and the thin dotted curves, where the 
THc limit has been taken for the latter. As has been theoretically anticipated in [21] 
the charm mass can be safely set to zero in QS^^^ and the involved convolutions in Eq. 
(4.1) may be replaced by the massless expressions in [29, 30] which simplifies the numerics 
essentially and which we will therefore do for the n = 2mc curve in Fig. 4.4 below. As 
also stressed in [21] it is, however, essential to keep the charm mass finite in the GF^^^ 
contribution since rric tempers the strength of the z'~^ u-channel propagator singularity. 
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Figure 4.4: A Comparison of the total 0{al) result (solid lines) to HI data [11] on (circles) 
and D^* (triangles; both including charge conjugation) production. Shown are results for two 
choices of the factorization/renormalization scale fx. Also shown is the outcome of a fixed 
order 0{al) GF calculation for comparison (dashed line). The dot-dashed line follows the 
suggestion in [20] and neglects quark initiated contributions at 0{a\), i.e. the difference 
between the solid (/x = Q) and the dot-dashed line is given by the (thick) dotted line in Fig. 
4.3. For the dot-dashed curve as well as for the fixed order calculation (dashed) only fJ, — Q 
is shown since the scale dependence is completely insignificant. 

Obviously the 0{al) result of an ACOT based calculation deviates essentially from the 
naive Born term expectation and it seems by no means legitimate to treat GF^^^ as a 
higher order correction here. Contrary to the corresponding inclusive results in Chapter 
3 and to the expectations in [20] also the subtracted QS'^^^ term is numerically significant 
in the semi-inclusive case considered here. In the light of the huge 0{al) corrections 
it seems, however, undecidable as to whether include or omit QS^^^ at 0{al) without 
knowing 0{a'l) corrections within ACOT. 

In Fig. 4.4 the resummed result (VFNS: sohd lines) can be compared to unsubtracted 
0(aiy^ GF^^^ (fixed order, dashed line). We show the total 0(al) VFNS result for 
the choices /i — Q,2mc. The modest scale dependence arises exclusively through the 
QS^^^ term. For any of the other curves a variation of fj, is completely insignificant 
and we therefore only show ji — Q. A full VFNS calculation (sohd lines) seems hardly 
^°An 0{al) NLO calculation [80] within FOPT gives very similar results [79]. 
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distinguishable from fixed order perturbation tfieory (dasfied line) within experimental 
accuracy. The two approaches are even closer if one follows the suggestion in [20] and 
does not include (dot-dashed line) the subtracted QS^^^ term at the level of QS^^^ + 
GF^^\ The data points in the figure correspond to HI measurements [11] of (circles) 
and D*^ (triangles; both including charge conjugation) spectra. The measurement is 
restricted to r^^ < 1.5. Extrapolating to the full phase space gives rise to large acceptance 
corrections which are, however, quite uniform [79] over the kinematical range considered 
and therefore have a minor effect on the normalized spectrum. Since FOPT and ACOT 
based calculations are very close it seems improbable that an experimental discrimination 
between the two approaches will be possible. The Born term in Fig. 4.3 is far from being 
the dominant contribution and an intrinsic c{x, Q^) stemming from the resummation 
of perturbative logs can therefore not be excluded. The tendency of the data appears 
somewhat softer than any of the calculations and the tendency seems to be confirmed by 
preliminary data in a lower z bin [81]. The rcsummed calculation appears to be too hard 
at larger z around 0.6 if the QS^-^^ component is included. As already mentioned, at the 
present stage of the calculations it cannot be decided whether this hints at an intrinsic 
problem within VFNS calculations or whether this may be cured by 0{a^) corrections. 

4.4 Conclusions 

In this chapter we have performed an ACOT [20] based analysis of heavy quark fragmen- 
tation in DIS including a calculation of semi-inclusive scattering on massive quarks at 
0{al) for general masses and couphngs. As in the inclusive case of Chapter 3, effects 
from finite initial state quark masses can be neglected for practical applications to charm 
production in CC and NC DIS. The involved convolutions in Section 4.2.1 can therefore 
safely be replaced by their analogues in [29,30] and in Appendix B of [50]. Neutrinopro- 
duction is an ideal environment to extract the charm FF within DIS and a Peterson [37] 
type FF with Sc ^ 0.06 seems to lie in the correct ball park, where the sensitivity on 
the choice of scheme is small and finite effects are irrelevant. The Sc value above is 
compatible with e'^e~ data if a nondegeneracy of charm quarks fragmenting into D and 
D* mesons is allowed for. For NC DIS it seems unhkely that a discrimination between 
fixed order and resummed calculations will be possible at HERA. Both approaches give 
similar results which show a spectrum that is somewhat harder than the tendency of the 
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data [11, 12,81]. The resummed calculation is made worse if the 0{al) quark scattering 
contribution is included at the perturbative level considered here. 



Chapter 5 



Summary 



In the first part of tliis tliesis wc liavc studied deep inelastic production of lieavy quarks 
in neutral current and charged current processes. For this purpose we have calculated the 
relevant partonic subprocesses to order 0{al) (BGF, QS) for general masses and couplings 
taking into account massive initial state quark-partons as needed in the variable flavor 
number scheme of ACOT [20]. Our calculation of the vertex correction with general 
masses and couplings is, to our best knowledge, new. 

The integrated partonic results could be used to formulate and investigate heavy quark 
contributions to inclusive deep inelastic structure functions within the ACOT variable 
flavor number scheme. By the calculation of the before missing radiative corrections to 
scattering amplitudes on massive quark partons (including virtual corrections) the ACOT 
scheme could be completed to full order 0{al). Furthermore, we utihzed the unintegrated 
partonic structure functions to extend the original ACOT scheme [20] from its inclusive 
formulation to (^-differential) one-hadron-inclusive leptoproduction. 

With help of these results we demonstrated in the charged current sector that the 
effects of a finite strange mass are small and rapidly converge towards MS in the limit 
rUs 0. Furthermore, it turned out that the effects of finite initial state masses can 
be neglected in all cases, such that the (massive) QS*^^-* contribution of Section 3.2.2 can 
safely be approximated by its much simpler MS analogue at any scale. Nevertheless, it 
should be stressed that radiative corrections to quark initiated processes are in general 
not negligible and should be included in complete NLO analyses employing heavy quark 
parton densities. 

Finally, we have performed an ACOT based analysis of heavy quark fragmentation in 
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DIS led by the observation that semi-inclusive CC DIS is well suited to extract the charm 
fragmentation function. We reassured that theoretical uncertainties due to scheme choices 
and finite nis effects are irrelevant and found out that a Peterson type fragmentation 
function with Ec ~ 0.06 seems to he in the correct ball park. The £c value above is 
compatible with e'^e~ data if a nondegeneracy of charm quarks fragmenting into D and 
D* mesons is allowed for. In the case of semi-inclusive NC DIS it seems unhkely that a 
discrimination between fixed order and resummed calculations will be possible at HERA. 
Both approaches gave similar results exhibiting a spectrum that is somewhat harder than 
the tendency of the data [11,12,81]. The resummed calculation is made worse if the 
0{al) quark scattering contribution is included at the perturbative level considered here. 



Part II 



The Structure of Real and Virtual 

Photons 



Chapter 6 

Introduction and Survey 



Photon physics is an active field of research as is documented by and in a large number of 
reviews [82-91]. Particularly in the past few years there has been much progress due to 
the wealth of experimental results from the e^e~ collider LEP and the ep collider HERA 
[90,91]. Among the reactions initiated by high energy photons such processes providing 
a hard scale are of particular interest since they can be described (at least partly) by 
means of perturbative QCD. More precisely, in this thesis we are interested in processes 
which can be described within the well known framework of the parton model and in 
this case the 'hadronic nature of the photon' is quantitatively described by the partonic 
structure of the photon. The classical way of measuring the photonic parton content is 
deep inelastic electron-photon scattering (DlSg^) which is the main source of information. 
The DlSe^, data on the photon structure function F2{x, Q"^) are mainly sensitive to the 
up-quark density u'iyX^Q'^) as can be seen from the parton model expression for F2 (in 
LO), F2 oc Au^ + + s'*', whereas the gluon distribution g"'[x,Q'^) is only indirectly 
constrained at small values of x due to the evolution. Complementary information has 
become available in the last years due to 'resolved photon processes' (e.g. production of 
(di-)jets or hadrons with large pt (-Et), heavy quark production, isolated prompt photon 
production) in 7p and 77 collisions at HERA and LEP, respectively, which are mainly 
sensitive to the gluon distribution in the photon. However, the experimental constraints 
on the gluon density are still weak, especially at x < 0.1, and one has to resort to model 
assumptions about the parton distributions of the photon. 

It is the virtue of the phenomenologically successful radiative parton model (GRV 
model) [26,41,57,92,93] to predict the small-x behavior of parton distributions by pure 
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DGLAP-evolution of valence-like input distributions at a low input scale jj? ~ 0.3 GeV^. 
Recently, the parton distributions of the proton have been revised [26] due to precision 
measurements of the proton structure function F| at HERA leading to slightly changed 
parameters (low input scale, as) and a less steep small-x increase compared to the previous 
GRV94 [41]. In order to test quantitatively (within the GRV approach) if protons, pions 
and photons show a similar small-x (i.e. high energy) behavior it is necessary to update 
also the parton distributions of pions and photons. 

At present e^er and ep collider experiments the photon beams consist of 
bremsstrahlung radiated off the incident lepton beam resulting in a continuous spec- 
trum of target photons 7(-P^) where is the photon's virtuality. The bremsstrahlung 
spectrum is proportional to such that the bulk of target photons is produced at 

~ -Pmin — 0- The parton content of such (quasi-)real photons is well established both 
experimentally and theoretically and it is quite natural to expect the parton distributions 
to decrease smoothly with increasing P^ and not to vanish immediately. In this sense the 
real photon 7 = 7(P^ — 0) is just a 'primus inter pares' and unified approaches to the 
parton content of virtual photons 7(-P^) which comprise the real photon case in the limit 

— > are highly desirable. This is also reflected by the fact that measurements of the 
real photon structure function in single-tag events integrate over the bremsstrahlung 
spectrum from P^j^ up to a P^ax which depends on the experimental details. For instance, 
at LEP1(LEP2) P^^ is as large as P^^ ~ 1.5 GeV2(4.5 GeV^), cf. Section 2.2 in [90]. 
AUthough the bulk of photons is produced at P^ ~ P^j^ the amount of ignorance of the 
P^-dependence, mainly in the range P^ < A^, feeds back on the determination of the 
structure function P^ (parton distributions) of (quasi-)real photons. 

It is the central goal of the second part of this thesis to perform LO and NLO 
analyses of the parton content of pions and real and virtual photons within the latest 

finc-tuncd/improved setting of the GRV model [26]. 

The outline of Part II will be as follows: 

• In Chapter 7 we provide the basic kinematical background for studying the structure 
of real and virtual photons in two-photon scattering events. Structure functions 
for virtual photons arc defined in a general (model independent) way which will be 
studied in the following chapters either in fixed order perturbation theory (Chapters 
8, 12) or within the framework of the QCD-improved parton model (Chapters 9, 



. . . , 12). Furthermore, we demonstrate the factorization of the cross section for the 
process ee eeX into a flux of photons times the cross section for deep inelastic 
scattering (DIS) on these "target" photons for arbitrary virtuahties of the target 
photon in the Bjorken hmit <^ Q^. It should be noted that the factorization is 
essential for a theoretical description of two-photon processes in terms of structure 
functions of target photons 7(-P^) which can be measured in deep inelastic electron- 
photon scattering. 

In Chapter 8 we calculate the photon photon cross sections according to the dou- 
bly virtual box 7*((5^)7*(-P^) — > qq in lowest order perturbation theory. The general 
expressions are casted in a form which easily allows to read off various important 
limits, e.g., the quark-parton model (QPM) results for the structure functions of 
real and virtual photons and the heavy quark contributions to the photon struc- 
ture functions. For deeply virtual target photons the perturbative results make 
reliable predictions and we compare them (Fes) with present e'^e~ virtual photon 
data. On the other hand, for quasi-real target photons these results are plagued by 
mass singularities which have to be subtracted and afford the introduction of parton 
distributions which will be done in the next chapter. 

In Chapter 9 the complete theoretical framework necessary for our phenomenological 
study of the parton content of real and virtual photons in Chapter 11 is provided 
where special emphasis is laid on a unified treatment of real and virtual photons 
also in NLO. 

Chapter 10 is devoted to an analysis of the parton content of the pion in LO and 
NLO QCD. Since only the pionic valence density v'^{x, Q"^) is experimentally rather 
well known at present, we utilize a constituent quark model [94, 95] in order to 
unambiguously relate the pionic light sea and gluon to the much better known 
(recently updated) parton distributions of the proton [26] and Q^). These 

results will serve via vector meson dominance (VMD) as input for the hadronic 
component of the photon in the next chapter. 

In Chapter 11 we turn to LO and NLO analyses of the parton content of real and 
virtual photons within the framework of the radiative parton model. Apart from 
utilizing the latest refined parameters of the radiative parton model [26] there are 



some novelties as compared to the original GRV^ approach [96]. At first, the bound- 
ary conditions for the real photon are based on a coherent superposition of vector 
mesons which maximally enhances the up-quark contribution to the photon struc- 
ture function F^(a;, Q^) as is favored by the experimental data. As a result no extra 
normalization factor for the photonic boundary conditions is needed. Furthermore, 
in order to remove model ambiguities of the hadronic light quark sea and gluon 
input distributions of the photon (being related to the ones of the pion via VMD), 
inherent to the older GRV^ [96] and SaS [97] parametrizations, we employ, as al- 
ready mentioned, predictions for the pionic light quark sea and gluon which follow 
from constituent quark model constraints. The resulting predictions for the real 
photon structure functions F2{x,Q'^) will be compared with all presently (January 
1999) available relevant data. Most recently the OPAL collaboration [98] at the 
CERN-LEP collider has extended the measurements of the photon structure func- 
tion -F^(a;, Q'^) into the small-x region down to a; ^ 10^^, probing lower values of x 
than ever before and we include a comparison with these data as well. Finally, we 
construct LO and NLO boundary conditions for the virtual photon which allow for 
a smooth transition to the real photon case and perform a careful study of the re- 
sulting predictions for the partonic content and the structure function F2^^ \x, Q^) 
of virtual photons. 

In Chapter 12 we test our model for the parton content of virtual photons in more 
detail and compare these QCD-resummed results with the fixed order box expres- 
sions of Chapter 8. It is demonstrated that present e+e" and DIS ep data on the 
structure of the virtual photon can be understood entirely in terms of the standard 
'naive' quark-parton model box approach. Thus the QCD-resummed (renormal- 
ization group (RG) improved) parton distributions of virtual photons, in particular 
their gluonic component, have not yet been observed. The appropriate kinematical 
regions for their future observation are pointed out as well as suitable measurements 
which may demonstrate their relevance. 

Finally, in Chapter 13 we summarize our main results and discuss open questions. 

Various limits of the doubly virtual box calculation and parametrizations of our 
pionic and photonic parton distributions have been relegated to the Appendices C 
and D, respectively. 



Chapter 7 

Photon— Photon Scattering 



In this chapter wc provide the basic kinematical background for studying the structure of 
real and virtual photons in two-photon scattering events. 

In Section 7.1 we give a short introduction into so-called two-photon processes, in- 
troduce the usual kinematical variables to describe them and derive the cross section for 
such events. The notation follows mainly the report of Budnev et al. [82] . Here we only 
present a selection of the topic needed throughout the thesis. Additional information can 
be found in [82] and, e.g., in the more recent reviews [85, 90]. 

Next, in Section 7.2 we define structure functions for virtual photons and relate them 
to the photon-photon scattering cross sections. These structure functions will be studied 
in the following chapters either in fixed order perturbation theory (Chapters 8, 12) or 
within the framework of the QCD-improved parton model (Chapters 9, . . . , 12). 

Finally, in Section 7.3 we demonstrate the factorization of the cross section for the 
process ee eeX into a flux of photons and the cross section for deep inelastic scattering 
(DIS) on these "target" photons. While the factorization is well known for (quasi-)real 
target photons with virtuality ~ wc deal with the general ^ case in the Bjorken 
limit. It should be noted that the factorization is essential for a theoretical description of 
two-photon processes in terms of structure functions of the (real or virtual) photon. 
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7.1 Kinematics 



The kinematics of particle production via photon-photon scattering in e~^e coUisions 



e (pi)e+(p2) e (p'i)e+(p2)7*(g)7*(p) e (Pi)e+(p2)X(px 



(7.1) 



is depicted in Fig. 7.1. The momenta of the incoming and outgoing leptons are denoted 



Pi 



P2 




W"^ = {q+ p) 



Figure 7.1: Two-photon particle production. The solid lines are the incoming and outgoing 
leptons and the wavy lines are virtual photons which produce a final state X consisting of 
hadrons (or leptons). 

by Pi = {Ei,pi) and p[ = {E[,pi) {i = 1, 2) respectively and the momenta of the photons 
are given by 

Q = Pi-Pi , Q = -q =Qi , , ^ 

(7.2) 

p = P2-p'2 , P^= -P^ = Ql . 

In general both photons have spacelike momenta and refers to the photon with smaller 
virtuality (P^ < <5^)- X denotes the final state produced in the 7*(g) + 7*(p) — > X 
subprocess. For later use, we define the following variables: 

r P-Q P- 1 

u =p.q ^ a; = — , = — , yi = , y2 



z/ 2z/ p ■ Pi Q ' P2 

W^ = {p + qf = 2z/(l -x-6) = Q^^—^ - P2 . (7.3) 

X 
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Pi 
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Pi 












^ ^ P,^' 










^ P2 


P2 



^(2) 



Figure 7.2: Squared matrix element of the process (7.1). Integration over the phase space of 
the system X is implied as indicated by the vertical cut. 

The cross section for the process (7.1) is given by 

1 



da = —\M\'dQ^ 

p.p. 



with the invariant matrix element M, the M0ller flux factor of the two incoming leptons 



4v^ (pi ■ P2Y — mlml 



and the Lorentz-invariant {n + 2)-particle phase space 



(7.5) 



^^^^(2vr)45(4) [q + p- px) dV . 



(7.6) 



Here px = ^ihji G X is the total momentum and dT = Yli 2^(2^)^ ' ^ ^ phase 
space volume of the produced system X. 

The cross section can be expressed in terms of the amplitudes M'^'^ of the 7*(g) + 
7*(p) ~^ X subprocess as follows (see Fig. 7.2): 



da 



d^p[d% (47ra)' 1 
2E[2E'^{2txY Q^P^ p^^^^l)^i2)Vy^.'u',,u 



(7.7) 
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with 

ly/^V',/.. _ 1 y M"'^'"'M^'^(27r)^5(^)(g + p- px)dT . (7.8) 
For unpolarized leptons the (leptonic) tensors and L(2) (see Fig. 7.2) are given by 

^(.? = iTr[(^, + me)7"(/, + me)7^]. (7.9) 

The factor 1/2 is due to a spin average over the incoming leptons. In addition, we 
introduce the dimensionless quantities [82] 



(2p. - g.)"(2p, - g,)/^ 



which have the interpretation of (unnormahzed) density matrices for the corresponding 
virtual photons. 

7.1.1 The Hadronic Tensor 1^^'^ '^^ 

According to the optical theorem W'^''^''^'^ is the absorptive part of the virtual 77 forward 
amplitude shown in Fig. 7.3. 




Figure 7.3: The photon-photon forward scattering amplitude T^''^''^'^; m, n, m' , and n' are 
helicity indices. The tensor W^" defined in (7.8) is the absorptive part of the 77 forward 
amplitude: W^'''''^^''' = ^ImT'''''''^'' . 

' TV 



Taking into account P- and T-invariance (symmetry /xV ^ fiu) and gauge invariance, 

i.e. 
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the Lorentz-tensor W^'''''^'^ can be expanded in terms of a basis of 8 independent tensors 
constructed from the vectors q, p and the metric tensor g. The choice of the tensor basis 
in which the expansion is carried out is clearly arbitrary and different forms are discussed 
in the literature [82, 99, 100]. In the following we will stick to the expansion given in [82]: 



TT 



- (i?^"gf g^' + R^'^Qt'Q'^ + R^'^'QtQl + R'''''Q2QiW:i^ 

- {R^^QI'q^^' - R''''Q'iQ2 + R^''''Q1Q2 - i^'^'^ga'gO W^TL 



(7.11) 



The tensor structures used in (7.11) are connected to the photon polarization vectors, see 
App. B and C in [82]: 



(7.12) 



with = 1 — Ax6 and 



R 



■ap _ _gap i^ip^q^ + q'^p^) - p^q'^q^ - q^p'^p^ 



1/2/32 



g? 



The photon momenta gf,p, the unit vectors Qi,Q2, and the symmetric tensor R'^^ satisfy 
the following (orthogonality) relations: 

g • gi = p • g2 = , g^ 2 = i , 

q'^R^p = p'^R^f, = Ql^R^p = , R^pR'^f = 2 , = -i?"^ . (7.14) 

With help of these relations it is easy to see that the various tensors in front of the 
invariant functions Wab in Eq. (7.11) are mutually orthogonal and therefore can be used 
to project out the invariant functions. With obvious notation {Pwtt'^'^' ''^'^^ix't^' ,1x1^ — W^tt 



(7.13) 
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etc.) the projectors read: 



rwxT 




4 








p m'*^' 






P M'*^' 




= Ql Q1Q2 Q2 5 


p m'*^' 






p m'^^' 




= ^(i?^-gf g^' + R^'^'Q^Ql + R^^'^'Q'^Q'^ + R^'^Q'^'Qt) 




p m'*^' 




= h^Ri''''' Ri"" - R'^'i'Ri'''') , 


p m'^^' 




= -i(i?^-gf g^' - R-^'^Qt'Q'^ + R^''''QtQ2 - R^'^'Q^Qi) 




(7.15) 



The dimensionless invariant functions IVafe depend only on the invariants M^^, g^ and 
and are related to the 77-helicity amplitudes^ Wm'n',mn in the 77-CMS via 



W^TT = -(W^++,++ + W^+-,+-) , 

Wlt = Wo+,o+ , 



LL 



00,00 



(7.16) 



-(W^++,oo - W^o+,-o) ■ 



The amphtudes W^r^, ^tl) ^tl correspond to transitions with spin flip for each of 

the photons (with total helicity conservation). As we will see in the next section only 

6 of these amplitudes (VFtt, W^tl, W^lt, W^ll, W^tt) ^tl) enter the cross section for 

unpolarized lepton beams because the tensors in (7.10) are symmetric whereas the tensor 

structures multiplying W^rj. and W^j^ in (7.11) are anti-symmetric such that these terms 

do not contribute when the leptonic and the hadronic tensors are contracted. Only if the 

initial leptons are polarized, can the amphtudes W^rj. and W^i^ be measured as well [82]. 

^The photon helicities can adopt the values m',n',m,n = +1,-1,0. Total hehcity conservation 
for forward 77-scattering impUes m' — n' = m — n and due to P- and T-invariance {Wm'n',mn = 

/ 1 \m'— n'+m— nTT/ Tl/ 21 TT/ 

V ^) —ra' —n' ,—'m — n — ^* —m' — n' ,—'m — n *^mn,m'n 



) there exist 8 independent helicity ampli- 



tudes. See App. C in Ref. [82] and Ref. [100] for further details. 
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7.1.2 Derivation of the Cross Section 

Using Eqs. (7.9), (7.10) and (7.11) one obtains by a straightforward (but tedious) calcu- 
lation 



-^(1) -^(2) f^'i^',f^u — V -r 



M^pt^^TT + 2|p+-p+- 1 COS 20 + 2pt+pfWTL 

(7.17) 

+ 2p?V2++iyLT + pTpfWi^L - 8|p+V^°|W^^L COS 



where is the angle between the scattering planes of the e~ and the e"*" in the center- 
of-mass system (CMS) of the colliding photons and the pi's are elements of the photon 
density matrix: 

9^++ _ 9^— _ ^"/3 D _ (2pi -P-P- qf ml 
2p, -2p, -p, R^,- (^.^)2_g2p2 +1-4^' 

2p2 -2p2 -P2 RaP~ (^.^)2_g2p2 V' 

2 

00 _ /o _ o ^++ - ■ 



P?° = pfQiagi, = 2p^-2 + 4^, 

.2 



P2 = PrQ2„Q2;3 = -2 + 4-^ 

2|p+-pM COS20 = ^ - 2(p++ - l)(p++ - 1) , 



8|p+Vncos 



4C (2pi • p - p • ?) {2p2-q-p- q) 



yQ2p2 (p . g)2 _ g2p2 

with C = (2pi - 5)"(2p2 - pfRap = -(2pi - g) • (2p2 - p) 
+ {p.q)2-Q2p2 C^Pi ■ P - P ■ (l) i'^P2-Q-p-q) , 

\pr\ = pr-h 



\pr\ 



i^ri- (7-18) 



Note that all these quantities are expressed in terms of the measurable momenta Pi,P2 
and p'i,P2 (respectively q,p) only and therefore are entirely known. 

With help of Eqs. (7.7) and (7.17) we easily find the fully general final result for the 



1.2 Photon Structure Functions 



61 



ee — > eeX cross section [82, 85, 90]: 



dV(ee ^ eeX) = 



a 



77 



E[E'2 167r4g2p2 p^^ 

„00 ^00 , 



^00 ^++, 



++^00, 



+ pTP2''(^hL + 2|p+-p+-|TTTCOs20 - 8|p+V^Vtlcos( 
with 



77 



1/2 



(7.19) 



Here the cross sections aab (used as a shorthand for cttt, ctl, clt, ctll, ttt, ttl, ^tt, 
t-pl) are identical to the corresponding structure functions Wab up to a division by the 
appropriate flux factor of the two incoming photons^, i.e. 

1 1 



C^ab 



W, 



ab 



2iy(3 



-W, 



ab 



(7.20) 



The cross section in (7.19) considerably simphfies in certain kinematical regions [82, 
85,90]. For instance, if both photons are highly virtual Eq. (7.19) can be evaluated in 
the limit Q'^,P^ » ml [90] in which some relations between the elements of the photon 
density matrix exist. Of special interest in this thesis is the case where one of the lepton 
scattering angles becomes small leading to a small virtuality f^i of the corresponding 
photon while the other photon provides a hard scale ^ 1 GeV^. In this hmit the cross 
section factorizes into a product of a flux of quasi-real target photons times the cross 
section for deep inelastic electron-photon scattering, see for example [90] . This process is 
the classical way of measuring the structure of (quasi-real) photons. The findings in the 
latter limit can be generalized to the case of photons with non-zero virtuality 7^ as 
we will see in Section 7.3. This allows to study the structure of virtual photons in deep 
inelastic e7(P^) scattering processes in a continuous range of the scale P^. 



7.2 Photon Structure Functions 



It is the aim of this section to relate the structure functions of a (virtual) target photon 
to the invariant functions Wab- The defining relations are generally vahd for arbitrary P^. 
However, they only have a meaningful interpretation as structure functions of a target 
^Note that a factor \ has already been absorbed into the definition of " in Eq. (7.8). 
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photon probed by a deeply virtual photon 7(Q^) in the limit <C Q^. The following 
expressions are simplified if one introduces the transverse components of a four-vector 
and of g^u 

1 



T _ 



q ■ X 

^^^^ 



T _ 



where 'transverse' refers to g: q ■ x^ = 0, q'^gj^i 



q''9^ 



q- 
0. 



(7.21) 



7.2.1 Structure Functions for a Spin— Averaged Photon 

Usually one introduces structure functions for a spin-averaged target photon. The cor- 
responding structure tensor can be obtained by contracting W^'j^'^/j^i, given in Eq. (7.11) 
with the metric tensor g"'^' . With the help of Eqs. (7.13) and (7.14) one obtains^: 



— Rn'fi 



TT 



+ Qifi'Q 



(7.22) 



~9^ 



I T T Q 



2T 



2L 



where W2T = Wtt + M^lt and W2L = H^tl + M^ll- 



Alternatively the spin-averaged tensor can be expressed in standard form in terms of 
the structure functions Fi = Wi and F2 = UW2: 



Comparing Eqs. (7.22) and (7.23) we find 

1 Q' 



(7.23) 



2xF; 



<7> 



W2T - -VTsL 



;7.24) 



These relations can be re-expressed in terms of the photon-photon cross sections Uah 
Wab/{2uP) (a, 6 = 2, L, T) [82, 85, 90] 



2xF; 



<7> 



A 2 ^ 

47r^a 

1 

47r%/3 



ctt — 2'^'^^ 



0"2T — 2'^2L 



(7.25) 



•^Of course the virtual photon has three (+,—,0) degrees of freedom. The factor 1/2 guarantees the 
conventional normaUzation in the real photon limit with only two (+, — ) transverse degrees of freedom. 
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Finally, Fl satisfies the usual relation 



(7.26) 



This can be seen by contracting W'^,'^^ with the polarization vectors of longitudinal 
probe photons given in Eq. (7.12) thereby employing again the orthogonality relations 
in Eq. (7.14) 



followed by the appropriate normalization: 



7.2.2 Longitudinal and Transverse Target Photons 



Since the fluxes of transverse and longitudinal virtual photons will turn out to be different 
(see Eq. (7.37) below) it is most convenient to introduce structure functions of transverse 
respectively longitudinal target photons (instead of spin-averaged target photons). The 
procedure is completely analogous to the one in the previous section using Eqs. (7.12)- 
(7.14) and for this reason the description will be brief. 

I. Transverse Photons 



With the help of Eq. (7.12) we can construct the structure tensor for a transverse photon 
target which can can be cast again into different forms 

K', = ^[<(p)e;(p) + 6*r(p)6::(p)]iy,v',,. = Ir^'^'w,,^',,^ 



(7.27) 



„1 p7T I TPl 



i7t 



and we can directly read off the structure functions: 



2xF^ 



7t 



1 

1 1 



TT 



w- 



2T 



fJaTT 

47r^Q! 

1 

47r2a/3 



(7.28) 



C"2T ■ 
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Repeating the steps in Sec. 7.2.1 to determine F^'^^, we obtain W^'^ — Wur implying 

= ^^2xWl^ = ^^F^-^ - 2xF^-^ . (7.29) 

II. Longitudinal Photons 

The structure tensor for a longitudinal photon target is given by (using again Eq. (7.12)) 
W;^^ = er'(p)e^(p)W^MV',^. = Q2'Q2W^^v',^. 

,1 



and we find the following result for a longitudinal target photon: 

1 q2 



(7.31) 



Finally, we have (as could be expected) 

ly^L ^ ^ ^7L _ _J_2a;^r7L _ ^2^7L _ 2xF^^ . (7.32) 

Further inspection of Eqs. (7.22)-(7.32) reveals a relation ["< 7 >= 7t— |7l"] between 
the spin-averaged, transverse and longitudinal target photons 

F<^>^p^.__p^ (i = l,2,L) 

which is a consequence of the completeness relation for spacelike photons (cf. [82], 
Eq. (B.l)): 

er{p)e\{p) + e''-'{p)e''_{p) - el^{p)e'^{p) = -g^^' + ^ . 

^ ^ ' ^ ' ' > ^ (7.34) 

27t - 7l = 2 < 7 > 



7.3 QED-Factorization 

It is well known that for the general cross section for the process ee — > 6677 — > eeX 
factorizes into a product of a flux of target photons (radiated off the electron) with the 
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Figure 7.4: Factorization of the ee 6677 — > eeX cross section into a flux of "target" 
photons radiated off the lower lepton line times the cross section for deep inelastic electron- 
photon scattering (black part). The cut in the photon line indicates a time order between the 
two subprocesses (photon emission followed by deep inelastic 67 scattering) and implies also 
that these two factors are independent of each other. 



deep inelastic electron-photon scattering cross section [85,90], see Fig. 7.4 for a graphical 
representation: 



da{ee ^ eeX) _ . daje^ ^ eX) ( p2 ^ r.. 



(7.35) 



dxdQMzdP^ .7T/.= v~- . ^^^Q2 

where z = E^/E ^ 1/2 is the fraction of the lepton energy carried by the photon (in the 
e"'"e~-CMS). The cross section for deep inelastic electron-photon scattering in Eq. (7.35) 
reads 

da{e-f ^ eX) Air'^ay a [1 + (1 - y)^ 1 



dxdQ'^ 



y 



(7.36) 



[(1 + (1 - y f)2xF;' + 2(1 - y)F^] 



with the usual variable y 



p-pi 



yi. Furthermore, /^^/e denotes the flux factor of 



transversely (or circularly) polarized photons with virtuality P^. For use below we also 
provide the flux factor /^^/e of longitudinally polarized photons: 

'1 + (1 - zY 1 2mh' 



a 
2^ 
a 
2^ 



Z P2 

2{l-z) 1 

i P2 



pi 



(7.37) 



The factorization in (7.35) is essential for relating the concept of the structure of a 
(real) photon to experimental measurements of two-photon processes. For this reason we 
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want to generalize Eq. (7.35) for photons with virtuahty 7^ and show that, in the 
Bjorken hmit, factorization holds for virtual "target" photons as well. For definiteness, 
as the Bjorken limit we consider 

Q^ = Ql^oo , u^oo , Q'^/2v = fixed . (7.38) 

Practically, this means P"^ = Q\ <^ Q^,^ such that 5 = P^/2i/ = xP'^/Q'^ is a small 
quantity which can be neglected. 

The starting point is the general cross section in Eq. (7.19). Employing Eqs. (7.28), 
(7.29), (7.31), and (7.32) and rearranging the terms inside the square brackets it can be 
written as 



(fa 



^(2pt^2pt^[2xFr+eFl-] 



+ p^°2p++[2xF7^ +2|p+-p+-|TTTCOs20-8|p]^V2 VtlCOS0 



(7.39) 



,nOO 

with 



Pi 



++ 



2pt 

The general strategy will be to demonstrate that P2 ^ and ^-^e proportional to the 
flux factors of transverse and longitudinal photons, respectively, radiated off an electron 
and that the interference terms disappear after having performed an appropriate angular 
integration. 

In the Bjorken limit it is useful to perform a "light cone decomposition" of the 
4-momenta of the two photons [82] : 

Q = C^Pl + ^qP2 + q± P = ^pPl + ^pP2 + P± 

with (7.40) 

p1 = 0, Pi-q±=Pi-p± = 0, pi-p2 = S/2 

where S is the square of the ee-CMS energy. The momentum fractions ^+ and 
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and the transverse momenta can be easily calculated: 

2pi • ? = = -Q' = -Q'/S 

2p2-q^ SC^ => = 2p2 ■ q/S yi) 

2pi • P = = 2pi • p/S y^) 

2p,.p^SC;^-P' it^-PVS (7.41) 

= -^CC - ?i ^ ?i = -^'(1 - C) 
= -siti- - pI Pi = -p'ii - 

g± -Px = -yf^^cos(p = Q^P^ {1 - - COS (p . 

Here is the angle between the scattering planes of the e" and e+ in the e+e~-CMS. 
Obviously ^+ is neghgibly small such that we can use^ 

P^C;P2+P±. (7.42) 

(On the other hand cannothe neglected in the Bjorken limit.) In the e+e~-CMS the 4- 
momenta of the incoming leptons can be written as pi = {E, 0, 0, E) and p2 = {E, 0, 0, —E) 
where E = ^ (neglecting terms of the order 0{^)). Since the transverse 4-vector is 
given by p±^ = {0,p±^,p±y,0) we can infer from Eq. (7.42) that E^^ = ^pE, i.e., in the 
e"'"e~-CMS is the energy fraction of the lepton energy transferred to the photon. For a 
real (P^ = 0) photon we recover the familiar relation p = ^~P2 between the 4--f^omenta p 
and p2 of the coUinearly radiated photon and its (massless) "parent" lepton, respectively 

Before turning to the photon density matrix elements in the Bjorken limit it is helpful 
to relate the variable v = p ■ qio , ^~ and the transverse momenta: 

2^ = ^C^,+(l + p) with p = ^ -2^x5(1 - e+)(l - COS0 , (7.43) 

where p oc VS is small in the Bjorken limit. Employing Eq. (7.43) we find in addition 

yi = ^ = e,+(i + p), z/2 = ^ = e;(i + p) • (7.44) 

Introducing the variables ui = q ■ (pi + P2)/VS — ^{^q + ^q) and a;2 = p • (pi + 

P2)/\fS = ^((7 + C,t) ^'^'^^ ])hase space can be written as [82] [Eq.(5.15b)] (up to terms 

"'Of course, in the calculation of quantities which are themselves small of the order (e.g. p]^ in the 
e+e~-CMS or p2j_ in the 77-CMS) ^+ must be taken into account. 
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of the order 0(^)) 



E[ E'2 S 



^dQ^dP^dujiduj2d(l) = -dQ^dC^dP'^dC;d(l) 



TT 



(7.45) 



-dQ'dyidP'd^-d(p{l + 0{p)) 



where the third equahty can be understood with the help of Eq. (7.44). 

In the Bjorken hmit the photon density matrix elements in Eq. (7.18) can be cast in 
a very compact form using the symmetric notation Qf = and Q2 = P^'- 



2pr = -Q 



pT = -Qi 

Vi 
Vi 



y^ Ql Qi 

2(1-1/.) 1 1 , 
2{l-yi) 1 2mly, 



(7.46) 



+ 0{5) . 



Vi QI Qt 

These results have to be expressed by the independent integration variables Q^, yi, P^, 
and Furthermore, from here on we identify yi = y. 

Obviously P2 ^ and pf are proportional to the flux factors of transverse respectively 
longitudinal photons in Eq. (7.37) 



pf = lp'^iM^;^n+o{p) 



(7.47) 



and similarly we can write 



= -Q' 
y 

. pT 



y Q\ 

2(1-?/) 



pT = -Q' 
y 



2(1 -y) 1 
y 



(7.48) 



2p++ 1 + (1 - yf 
where we have discarded the mass terms due to Q"^ ^ ml . 

Inserting Eqs. (7.47), (7.46), and (7.45) into Eq. (7.39) one straightforwardly obtains 
(using ^ = 2uP/S = ^-y(3) 



(icr(e7T eX) 



+ 0{p) + interference-terms 



s=^-s + [Tt ^ 7l] 



(7.49) 
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with the cross sections for deep inelastic electron-photon scattering given by [cf. 
Eq. (7.36)] 



(i(7(e7T,L eX) , An'^a a 



y 



(7.50) 



where e can be found in Eq. (7.48). Note that only two of the variables x, y, and Q"^ are 
independent since they are related via — Sxy. The terms proportional to p vanish 
after ^-integration like terms of the order 0{5). 

Unfortunately, the interference terms are proportional to cos and cos 20 where is 
the angle of the electron scattering planes in the 77-CMS while is the angle of the 
electron scattering planes in the e''"e~-CMS (= laboratory system). However, we show 
below that 

COS0 = COS 0(1 + 0{p)) . (7.51) 

Therefore, we can also get rid of the interference terms (proportional to cos ~ cos 0) by 
integrating over 0, leaving a remainder of the order 0{\/5). The latter becomes obvious 
if we remember that the variable p defined in (7.43) is proportional to a/5 cos such that 
terms pcos0 occurring for example in Eq. (7.51) are proportional to -\/^cos^0 which do 
not vanish by integrating over 0. 

Before calculating cos let us state the final factorization formula: 

da{ee eeX) _ ^ daiejT-^eX) ,nf;;\ 

A few comments are in order: 



,2.c(a(e7L_^eX) ^ ^^^^^ (7.52) 

> dydQ^ '^=«^"^ 



• In the limit — the contributions from longitudinal target photons have to 

vanish since a real photon has only two transverse physical degrees of freedom and 

we recover Eq. (7.35).^ 

^More precisely, assuming that the structure functions of a longitudinal target photon vanish like 
P'^ /A? where A is a typical hadronic scale, say A = rup, we obtain a non-zero result because J^^/e oc 1/P^ 
which is, however, negligible compared to the contribution from transverse target photons. 
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• The factorized result in Eq. (7.52) is valid up to terms 0{\/S) which formally go 
to zero in the Bjorken limit. However, for practical purposes it is not clear when 
S — xP'^/Q'^ is small enough for Eq. (7.52) to be a good approximation of the 
exact cross section in Eq. (7.19). Here, a numerical comparison of the factorization 
formula with the exact cross section would be interesting. 

To complete our derivation we still have to show that cos = cos + 0(p) where cos 4> 
is given by Eq. (A.4) in [82] 

cos0= with Pi'i^-Pi.R'''{q,p) (7.53) 

and where R'^^{q,p) has been defined in Eq. (7.13). Using the decomposition in Eq. (7.40) 
it is straightforward to obtain 

Pi. ■ ?± 



Pi±-P2± 

Pii = -S(l-C) + ^(P) (7.54) 

P2l = -^(l-ep) + 0(p) ■ 

Inserting these relations into Eq. (7.53) and comparing with the definition of cos0 in 
Eq. (7.41) we find the above stated result 

cos 4> = + 0{p) = cos + 0{p) . (7.55) 



Chapter 8 

The Doubly Virtual Box in LO 



In this chapter we calculate the invariant amplitudes Wat (or the cross sections (Tab) in 
lowest order perturbation theory and compare them with present e~^e~ virtual photon 
data. These expressions are usually referred to as box results due to the diagram repre- 
senting the tensor W^''^''^'^, where the photons are attached to a fermion box, see Fig. 8.1. 




Figure 8.1: Box-diagram for 7*(g, /i)7*(p, z/) 7*(g, /i')7*(j9, z/') .There are 4 possibilities 
to attach the photons to the vertices. (2 for the initial state times 2 for the final state.) 

The calculation is a straightforward application of Feynman rules and will only be 
shortly outlined here. In order to obtain W^''^''^'^ according to Eq. (7.8) one has to build 
up the amplitude M^^" of the process 7*(g, ^)'~f*{p, v) — i{ki)l{k2) shown in Fig. 8.2 where 
f is either a lepton or a quark of mass m. The kinematics can be described with the help 
of the 4-momentum conservation relation q + p = ki + k2 and the usual Mandelstam 
variables 



s = {q + pf, t = {q-kif = {p-k2f , u={p-kif = {q-k2f , 



8.1 Unintegrated Structure Functions 



72 




Figure 8.2: Amplitude M^^" for the process 7*(g, Ai)7*(p, i^) f(/ci)f (/i;2) where f is a fermion 
of mass m. 

Si = 2q-p = s + Q'^ + P'^ , ti=t-m^ , tii = m - (8.1) 

satisfying si + ti + Mi = 0. 

The external fermion lines are on-shell, i.e. k\ = k\ = rri^, whereas the two photons 
are virtual (space-like): q"^ = —Q^ < 0, = —P^ < 0. The Dirac traces occurring 
in W^'"''^'^ (due to the closed fermion loop) have been evaluated with the help of the 
Mathematica [101] package Tracer [102]. Finally, the individual structure functions Wab 
have been projected out using the projection operators given in Eq. (7.15). 
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The general structure of the boson-boson fusion cross section (in lowest order) is given 
byi 



ah 



dzi 



levr^a^e^AT^ 

q L 



+ 



B 



^1 



1-^1 



+ 



D 



w„,. 



where zi = 1 + ti/si, A'^c = 3 is the number of colors and Cq is the quark charge. The QED 
case can be obtained from Eq. (8.2) by setting N^Cg —>■ 1. Note that zi is very similar to 
the fractional momentum variables used in Part I of this thesis to describe energy spectra 
of heavy quarks (mesons): 



ki -p 

q • p 2q ■ p 



^The i— channel amplitude has a propagator 1/ti and therefore has the form (a + bti)/ti. The 
u— channel is obtained by ii <-> ui = — (si + ti) and thus reads (a + hui)/ui = (c + dti)/{si + ti). 
Squaring the amplitudes generates the structure in Eq. (8.2). 
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On the other hand, Zi — —Ui/si (or 1 — zi — — ti/si) can be viewed as a natural 
dimensionless combination of the Mandelstam variables of the problem. Furthermore, it 
is convenient to use the dimensionless variables 

4m^ 



Q2 Q2 p2 p2 

Si Iv S\ Iv 



(8.3) 



in order to write the coefficients in a form which is manifestly symmetric under a; ^ 5 ^ 
and which easily allows to read off the important "massless limits" ^ or vr? ^ to 
be discussed in Appendix C. 

With /9 = Vl — 4x5 the coefficients read: 
-1 



^Wtt ~ 

CvFtt ~ 



327r/35 
-1 

4x5(x' + 5") + 8x^5" [1 + (1 _ ^ _ 5)2] ^ 

^i- [l - 2x{\ -x)- 25(1 - 5) + 4x5(1 - 2x)(l - 25) 
47rp^ L 

5(1 - x - 5) 



2x/32 - (1 - X - 5) (4x5 + A/?^) 25/3^ _ (i _ ^ - 5) (4x5 + A/?^) 
A2^^(1 - X - 5)2 - 2A^2(1 - X - 5)2 - 2 fl - 2x(l - x) - 25(1 - 5) 





47r/?5 




5(1 


— X — 


^) 




-47r/^5 




5(1 


— X — 





X 



2x/52- (l-x-5)(4x5 + A/52) 



2x^52(1 - X - 5)2 



A/52(|l + 2x(-l + x-5)) + 4x(|- 1 + 2X + 25- 2x5(1 + X + 5)) 
(1 - 2x)2 



2x5(1 - X - 5)2 



8x5(1 - X - 5)2 



1 + 2x5) 



-1 



327r/55 



;i-x-5)2[4x5 + A/3' 



212 



^Note however that neither Wtl nor Wlt but only the sum Wtl + Wlt is invariant under exchanging 
X and i5. 
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6xSj^ 
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{x + 6 - ixdf 



- 2x){l - 26){4x6 + A/?^) 

Vx5{l — X — 5) 



x6{l — X — 5) 



Ax6{-3 + x + S) + 2x + 26- X(3^{1 -x- 6) 
1- Ax -4:6 + 12x6) 



El 



W4 



- 2x6(3"^ + (1 - X - 6){4x6 + A/3^) 



167r/?3 
_(l-2a;)(l-25) 

(l-2x)(l-25) 



W2 



Cvi/^. — 



— X — 5) 
x5(l — X — 5) 



-4a;5 + X(3^) 



x6 



j:5(l — a; — 5) 
27r/g3 



(8.4) 



The coefficients of W^lt can be obtained from the corresponding ones of Wtl by exchanging 
X 6: Awj^T = Awt~j^[x ^ 6], Cw^^ = Cwt:A^ ^ 6], and Ew^^ = -Ew^tlI^ ^ 

It is noteworthy that our results in Eq. (8.2) generahze the ^-differential expressions 
for the (QED) structure functions F], F^, and F^ of real photons given in [103] to the 
P V case. 
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8.2 Inclusive Structure Functions 

The desired inclusive structure functions are obtained by integrating over the kinemati- 
cally allowed range in Zi. The boundaries for the -integration are given by 

z^,^^{l±m/2 (8.5) 

with = 1 - 4mVs = 1 - A and ^ 1 - Ax^P^/Q'^ = 1 - Ax5. 
Noticing that 

zi,+ - ^1,- = , 1 - ^1,+ = ^1,- , 1 - ^1,- = ^1,+ , 

= (1 - - ^1,-) = ^ (8.6) 

the required integrals can be immediately obtained 

13^ 4/3/3 





dzi 


/"zi,+ 




1 Zl,- 


{^-zif 








dzi 






Jzi,- 


/■zi,+ 




^1 



2 



- 4a;(5 + A/3' 
[n^ = lnl±^ (8.7) 

dzi — (3(3 . 

Now Eq. (8.2) can be integrated using Eq. (8.7) and we arrive at the following result 
for the inclusive structure functions expressed by the coefficients given in the previous 
Section 8.1: 

Wa, = lQ>^a'e%Q{(3') \2Cw^^L + /3^ {^^w^. ^ + E^^^) ] (8-8) 

with 

L = l.l±M. (8.9) 
1-/3/3 ' ' 

The ©-function guarantees that the physical threshold condition s > Am^ is satisfied.^ 

Recalling the relation aab = '^^ab we can rewrite Eq. (8.8) for the photon-photon 
cross sections 

f - / 4 \ 1 IQn'^a^N 

... = « + /3/3 (2^».. + £»..) I . iv - . 

(8.10) 



^It is easy to see that > also implies /J^ > 0. 
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Inserting the coefficients given in Eq. (8.4) into Eq. (8.10) one obtains the final result for 
the doubly virtual box 7*((5^)7*(-P^) — > qq in leading order: 



_ ^ 1 



1 - 2x{l -x)- 25(1 -5)- Ax5{x^ + 5^) + Sx^S"^ [l + (1 - a; - 5f] 



+ X(3\l -x-Sy- -X'(3\l -x-6) 



L + (3p 



4a;(l - x) -1 



+ 45(1 -6)- 8x5(1 -x'- S') - (4x5 + Xf3'){l - x - Sy 



4x5 + A ^2 



4 



-Xp^(^l-2S{l + x-S)^ - 25(^-1 + 2x + 25 



2x5(1 + X + 5)) 



x(l - 65 + 65^ + 2x5) + 5/3^ 



4x5 



4x5 + A /?2 



_ ^ 1 

'^^^ ~ 4^^ 



^x5(^l - 2x(l - x) - 25(1 - 5) + 2x5(1 - 2x - 25 - x^ - 5^ + 6x5)) 
2Xp^{l - X - 5) (^x + 5 - 2x5(1 + x + 5)) - ^A2;34(1 - x - 5)' 



(1 - X - 5)2(4x5 + X(3^) + 2(x + 5 - 4x5)^ 



N 2 r^. 



2x + 25 - 4x5(3 - x - 5) - Xp'^{l - x - 5) 



+ /3/3(l-3x-35 + 8x5) 
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N 1 



AxS + A/32 



+ 3 -4x- 45 + 4x6 



— 

'T'T 



N 4 



4x6 



(8.11) 



This recalculation is in agreement with the results of Ref. [82] (with iV^e^ 1) with 
exception of a relative sign between the part containing the logarithm L and the part 
proportional to (3(3 in t^l^- 

A derivation of various important limits of the doubly virtual box expressions in (8.11) 
can be found in App. C. 



8.3 Comparison with Present e+e Virtual Photon 
Data 



We will now compare our LO-box expressions with present e+e virtual photon data. 
The physically measured effective structure function in the Bjorken limit is [82,85,90] 



Feff(x; Q^yi; ^2) 



1 



CTtt + £(Z/l)0-LT + £(?/2)0-TL + £(?/l)£(|/2)o-LL (8.12) 



4'K'^a(3 

where the kinematical variables have been given in Eqs. (7.2) and (7.1) and where e{iii) 
are the ratios of longitudinal to transverse photon fluxes, 

£(y,) = 2(l-y,)/[l + (l-y^)1- (8.13) 

Furthermore, the photon-photon cross sections — crab{x,Q^, P^) with a — 

(L, T), b = (L, T) have been defined in (7.20). In the following we shall consider the 
kinematical region <S 1 relevant for double-tag experiments [105, 106] performed thus 
far where Eq. (8.12) reduces to 



Feff(x,Q^P^)~ 



Q- 1 
4'K^a P 



(7tT + Clt + ctl + Cll 



(8.14) 



"'This relative sign has also been noted in [104] where in addition the overall sign in ttl is different. 
Concerning the latter, we agree with the results of [82] . 



8.3 Comparison with Present e'^e Virtual Photon Data 



78 



Recalling the general definitions of the structure functions of spin-averaged, transverse 
and longitudinal target photons in Sec. 7.2 we can finally write 

= F<-(^^)>(a;,g^) + ^F-(^^)(a;,g^). (8.15) 

So far, our results are entirely general. 

We shall furthermore introduce the decomposition 

= 5 + (8.16) 

with afl,'^'* denoting the light (heavy) quark q = u, d, s {h = c, b, t) contributions. The 
light u, d, s contributions to o"^;, are obtained from Eq. (8.11) by setting m = nig = 
(A = 0) and summing over q = u, d, s. (Note that the box expressions involving a 
real photon, 7*(Q^)7(P^ = 0) ^ qq, require on the contrary a finite regulator mass 
m = rriq ^ 0; here one usually chooses to be, somewhat inconsistently, a constant, 
i.e. Q^-independent effective constituent mass, ~ 0.3 GeV.) For each heavy quark 
flavor h — c, b, t the heavy contribution cr^;, in (8.16) is obtained from Eq. (8.11) with 
Cq = Cfi and m = ruh- Only charm gives a non-negligible contribution for which we choose 
rUc — lA GeV throughout. 

Finally, it is instructive to recall the asymptotic results of our virtual (P^ 7^ 0) box 
expressions for the light q — u, d, s quarks derived from (8.11) in the Bjorken limit 
PVg2 ^ I App. C.2: 



4t ^ A^c(Se^)^x + (1 - xf] In + 4x(l - x) 



4l ^ <^It ^ A^c(SeJ)^x [Axil - x)] 

- . (8.17) 

We now turn to a comparison of the effective structure function Feff obtained in lowest 
order perturbation theory with all presently available e'^e~ data of PLUTO [106] and the 
recent one of LEP-L3 [105]. In addition to the full LO-box calculation of Feflf (solid hue) 
we show the 'asymptotic box' results (dotted line) where the light quark contributions 
have been calculated according to (8.17). As can be seen the 'full' and the 'asymptotic' 
curves are very similar and coincide in Fig. 8.4 over the entire x-range. Furthermore, 
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Figure 8.3: Predictions for Fes as defined in (8.14). The liglit (u, d, s) and heavy (charm) 
contributions in (8.16) of the 'full box' expressions in (8.11) are calculated as explained in 
the text below Eq. (8.16). The 'asymptotic box' results refer to the light quark contributions 
being given by (8.17). The PLUTO data are taken from [106]. 



the charm contribution (dash-dotted line) is shown separately to demonstrate its relative 
importance due to the charge factor in Eq. (8.11) suppressing contributions from down- 
type quarks by a factor of 16 relative to the up-type quarks. Of course, in Fig. 8.3 the 
charm contribution is restricted by the available phase space for producing a cc pair. 
Generally, the LO-box predictions for Fcs in (8.14) shown in Figs. 8.3 and 8.4 are in 
agreement with the present low statistics data.^ This is not unexpected in the case of the 
L3 data shown in Fig. 8.4 since the target photon is deeply virtual, = 3.7 GeV^, such 
that perturbation theory is applicable. On the other hand the PLUTO data (Fig. 8.3) with 
P^ = 0.35 GeV^ are just in the transition region from deeply virtual (P^ ^ A^) to real 
(P^ = 0) photons where non-perturbative effects are expected to become increasingly 
important, especially in the small-x region. Clearly, more precise data for Peff in this 
'transition region' with intermediate target virtualities are highly desirable. 

Two questions naturally arise: 

^Note also, that the rightmost data point in Fig. 8.4 is already close to the boundary of the phase 
space XmaK = (1+ P^/Q^)"^ 
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Figure 8.4: As in Fig. 8.3, but for Q"^ = 120 GeV^ and P^ = 3.7 GeV^ appropriate for the 
LEP-L3 data [105]. 



(i) Firstly, it would be interesting to see if NLO corrections to the doubly virtual box 

[107, 108] can further improve the description of the data. 

(ii) The other interesting question is whether for » -P^ ^ it is necessary to 
resum 'large' coUinear logarithms InQ^/P^ occurring in Eq. (8.17)^ to all orders in 
perturbation theory via renormalization group (RG) techniques. This issue will be 
addressed in Chapter 12 where we compare RG-resummed parton model expecta- 
tions for with the fixed order box results presented here. 

It should be emphasized that, at > > A^, the InQ^/P^ terms in (8.17) 

are finite which therefore need not necessarily be resummed but only if a resummation 

is phenomenologically relevant. This is in contrast to the situation for the real photon 

(p2 _ y^itii its well known mass singularities In Q'^/m'^, see Eq. (C.14), which afford the 

introduction of scale dependent (RG-improved) parton distributions which are a priori 

unknown unless one resorts to some model assumptions about their shape at some low 

^It has, however, already been noted in [104] that the term containing InQ^/P^ dominates only if 
(Q^/P^) is very large. 
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resolution scale (see, e.g., [6] and the recent reviews [89,90]). 

In the next chapter wc present the necessary theoretical background for describing 
(virtual) photon structure functions within the QCD-improved (RG-improved) parton 
model. 



Chapter 9 



The Partonic Structure of Real and 
Virtual Photons: Theoretical 
Framework 

The following chapters are mainly devoted to an analysis of the structure of the (real and 
virtual) photon within the framework of the renormalization group (RG) improved parton 
model. It is the main purpose here to provide the necessary theoretical background for 
these phenomenological studies where special emphasis is laid on a unified treatment of 
both real (P^ = 0) and virtual (P^ ^ 0) photons. 

After motivating our unified approach in Sec. 9.1 we turn to a detailed description of 
the unpolarized photon structure functions within the QCD improved parton model in 
Sec. 9.2. This section provides the complete technical framework for our analysis of the 
parton content of real and virtual photons in Chapter 11. In Sec. 9.3 we deal with the 
RG-based evolution equations for photonic parton distributions and their solutions. For 
completeness, we discuss boundary conditions for a deeply virtual target photon and from 
the evolution equations we derive leading and ncxt-to-lcading order expressions for the 
momentum sum carried by the photon's quark and gluon densities. Finally, in Sec. 9.4 
we briefly discuss how to include contributions from longitudinal target photons. 

9.1 Introduction 

We have seen in Chapter 7 that the cross section for the process ee — > eeX factorizes in the 
Bjorken limit into fluxes of transversely (longitudinally) polarized photons and the cross 
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Figure 9.1: Deep inelastic electron-photon scattering (cf. Fig. 7.4) in the parton model where 
the target photon is produced by the bremsstrahlung process. According to the factorization 
theorems [68, 109] the 67 cross section factorizes into 'long distance' (non-perturbative) parton 
distributions (grey oval) describing the internal structure of the photon and 'short distance' 
(perturbatively calculable) electron-parton subprocesses (black oval). 

sections for deep inelastic scattering on these transverse (longitudinal) target photons. 
The deep inelastic scattering cross sections can in turn be expressed in terms of two 
independent structure functions, e.g. Fj'"^^'^ \x,Q'^) and F^^^^ \x,Q'^) {F"^^^^ \x,Q'^), 
Fe^(^')(x,Q2)). 

Furthermore, in Chapter 8 we have obtained these structure functions (or the effective 
structure function Fcs) according to a perturbative calculation of the doubly virtual box 
7* ((5^)7* (P^) — > qq in lowest order perturbation theory. As has been discussed there in 
the deeply virtual case ^ (where A is a hadronic scale) perturbation theory is 
expected to be reliable and we have compared these fixed order 'box results' F^^^ with 
presently available e^e~ data leaving the open question if resummations of the collinear 
logarithms InQ^/P^ might further improve the situation. On the other hand, for < A^ 
the perturbatively calculated structure functions (cross sections) are invalidated by non- 
perturbative (long-distance) mass singularities which have to be subtracted. This problem 
can be consistently solved by turning to a parton model description (see Fig. 9.1): Since 
we are considering the Bjorken limit, the structure functions can be described in terms 
of photonic parton distributions, quite similar to the hadronic case where factorization 
theorems [68, 109] state that the structure functions are given by convolutions of par- 
ton distributions with Wilson coefficient functions. The parton distributions contain the 
long-distance part of the cross section and therefore can (in general) not be calculated 
perturbatively and have to be fixed by experimental information or model assumptions. 
On the other hand the coefficient functions involve only short-distance physics and are 
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calculable in perturbation theory. The main value of this approach lies in the fact that 
the parton distribution functions (PDFs) are universal, i.e. do not depend on the partic- 
ular hard scattering subprocess. Once fixed by experimental information from, e.g., deep 
inelastic electron-photon scattering the PDF's can be used to make predictions for other 
processes, for example photo-and electroproduction of jets or heavy quarks. 

Due to QED gauge invariance the contributions (structure functions, cross sections) 
from longitudinal target photons have to vanish in the real photon limit (P^ — > 0). More 
specifically, they are suppressed by a factor divided by a typical hadronic scale, say 
P'^/nip^ [82], for not too large virtualities of the target photon. Therefore, we will 
mainly focus on the structure functions of transverse target photons. The contributions 
from longitudinal photon targets will be considered at the end of this chapter. In the 
following we can therefore drop the index T and use instead 7(-P^) to denote transverse 
target photons. Our main interest will lie on F2^^ \x,Q'^) since up to now the cross 
section has been measured only for small values of y where the contribution of F^{(x y^) 
is neghgible. 

The following parton content can be assigned to a transverse photon target: 
q^^^'\x,Q^), g'^^^'\x,Q^), and r^^^^\x,Q^) where q and g are the photonic analogues 
of the usual quark and gluon parton distributions inside hadrons. In addition the tar- 
get photon also contains an elementary photon-parton F'^^^'^^ reflecting the fact that the 
photon is a genuinely elementary particle which can directly enter a partonic subprocess. 
As usual in a massless parton model approach (justified by the factorization theorems 
[68, 109]) the partonic subprocesses are calculated with on-shell partons partly due to the 
P^/ Q'^ power suppression of effects invoked by the off-shellness of the partons (which is 
of the order 0(P^)). While this is quite familiar for the quarks and gluons the latter 
rule also refers to the photon-parton V'^^^'^\x,Q'^) in the target photon, i.e., the 'direct' 
subprocesses, e.g., 7*((5^)r'''(^^^ — >■ qq have to be calculated with P^ = (cf. rule (ii) in 
[6]). Consequently we can use the same (massless) short distance coefficient functions 
in our calculations irrespective of P^. All effects due to a non-zero target virtuality are 
entirely taken care of by the P^-dependence of the parton distributions (and by 

the Weizsacker-Williams flux factors) which of course have to be given in the same fac- 
torization scheme as the coefficient functions in order to obtain a physically meaningful, 
i.e. factorization scheme independent, result. By construction a physically smooth behav- 
ior in P^ (including a smooth transition to the real photon case P^ = 0) of the observable 
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structure functions can be automatically achieved by demanding a smooth behavior of 
the boundary conditions for the parton distributions. 

The heavy quark contributions deserve special care due to the additional scale nih 
provided by the heavy quark mass. In this case terms P"^ /m\ are possibly not negligible 
and have to be kept in the heavy quark coefficient functions. We will come back to this 
point in Sec. 9.2.2. 



9.2 Photon Structure Functions in the QCD- 
Improved Parton Model 



We shall introduce the decomposition 

F?^'''\x,Q') = Flf\x,Q')+F?p{x,Q') (9.1) 

with -^i'^'^h)^ denoting contributions with light quarks q — u,d,s (heavy quarks h — c,h,t) 
in the final state. In NLO(MS) F^'^^^\x, Q^) and F^f^\x, Q^) are given by the following 
expressions: 



X 



^^Q") - E e;|g-(^^)(a;,Q^)+g--(^^)(x,Q^) + 

q=u,d,s 



27r 



X [C,,2 (g + g)^^^') + 2 C,,2 ® ^^^^')] + ^ I (9.2) 



Flp{x,Q^) =1 Yl ej|g-(^^)(x,g^) + g-(^^)(x,g^) + 



27r 



+ -e C^,i(xJ 

TT 



(9.3) 



q=u,d,s 

X [Cq,i^{q + qr^'''^ + 2Cg,i®g^^'''' 

where (8) denotes the usual convolution integral of two functions / and g defined on the 

interval [0, 1]: 

{f<S)g){x)^ I dxi I dx2f{xi)g{x2)S{x - X1X2) = [ —f{z)g{x/z). (9.4) 

Jo Jo Jx ^ 

Here q^^^^\x,Q'^) = q'^^^^\x,Q'^) and g'~^^^^\x , Q"^) provide the so-called 'resolved' con- 
tributions of 7(P2) to Py^^^^ with the MS coefficient functions [30, 48, 54] 

C,,2(a:) = C,,i(x) + ^2x 



4 
3 



l + x^ (1-x 3\ 1 ,^ 



9.2 Photon Structure Functions in the Q CD-Improved Parton Model 



86 



Cg,2{x)=Cg,,{x) + -Ax{l-X) 



(x^ + (1 - X?) In + 8x(l -x)-l 

X 



(9.5) 



while C-y,2,i in (9.2) provides the 'direct' contribution as calculated according to the 'box' 
diagram 7* ((5^)7 — > qq [HO]: 

3 



C^^i{x) 



(9.6) 



(1/2) 

with i = 1,2. The convolution with the []+ distribution can be evaluated using, for 
example, Eq. (A. 21) in Ref. [57]: 



Jx y \y 



X 



9{y) — 9{x) 

y 



9{x) / dyf{y) . 
Jo 



(9.7) 



The coefficient functions C^x and ^01 the longitudinal structure function 



7{P2) 



2xF^^^ ■* may be deduced from Eqs. (9.5) and (9.6) and are given by [111]: 



4 1 3 

Cq,L(^) = 3 ' Cg,i^{x) = - 4x(l - x) , C^,L(a^) = JYJ2) ^^'"^^^^ ■ 



9.2.1 Scheme Choice 



In order to avoid the usual instabilities encountered in NLO(MS) in the large-x region due 
to the ln(l — x) term in C^,2(a^) in Eq. (9.6), we follow Ref. [96, 112] and absorb such terms 
into the photonic MS quark distributions in Eq. (9.2): this results in the so-called DIS^ 
factorization scheme which originally has been introduced for real photons by absorbing 
the entire 'direct' C^,2 term appearing in Eq. (9.2) into the NLO(MS) quark densities 



{q + qV£,' = {q + qy^^' + el^C,,,{x) 

yois^ - 



y(p2) 



a 

TT 



(9.9) 



with C^^2{x) given by Eq. (9.6). How much of the 'finite' terms in Eq. (9.6) is absorbed 
into the MS distributions in Eq. (9.9), is of course arbitrary and a matter of convention 
[113]. Since such different conventions [113, 114] differ by terms of higher order and turn 
out to be of minor importance for our quantitative results to be discussed in Chap. 11, we 
prefer to stick to the original DIS^ scheme [96, 112] as defined in Eq. (9.9). Furthermore, 
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the redefinition of the parton densities in Eq. (9.9) imphes that the NLO(MS) sphtting 
functions k'ij^l {x) of the photon into quarks and gluons, appearing in the inhomogeneous 
NLO RG evolution equations for p^^^^XjQ"^), have to be transformed according to [96, 
112,115] 

^«|Dis.=^W-2^e^p(;)0C,, (9-10) 

where 

(a;) = 1 1 1 _(i _ 2x) In^ x - (1 - 4x) In x + 4 ln(l - x) - 9x + 4 



2 In^ a; + 2 ln^(l - + 4 In x - 4 In x ln(l - x) 
2 



+ {x^ + (1 - xf\ 

- 4 ln(l - a;) + 10 - ^ TT^ 

A;W(x) = 3 Ve^^ |-2(l + x) ln2x-(6 + 10x) Inx + ^ + ^x^ + Sx-iel (9.11) 
</ ^3 1 3x 3 I 

q ^ ^ 

with k^q^ referring to each single (anti)quark flavor. The LO splitting functions are given 
by ^i? = I (^)^ and Pi? = I [1 + (1 - xf\ /x. 

In NLO the expression for F^^^ ^ in the above DIS-y factorization scheme is given by 
retaining the Cq^2-i Cg^2 terms while dropping the destabilizing 2 term in Eq. (9.2), which 
has already been absorbed into the quark densities according to Eq. (9.9). Similarly, 

DIS 

Cq^i and Cg^i remain unchanged while C^,i has to be modified according to ^ ^ — 
C^-C^=-12x(l-x). 

The LO expression for i^^*-^ ■* (i = 1,2) are obviously entailed in Eq. (9.2) by simply 
dropping all NLO terms proportional to Cg^i, Cg^^ as well as C-y^i. 



9.2.2 Heavy Flavor Contributions 

The photonic quark distributions discussed thus far and which appear in Eq. (9.2) are 
adequate for the / = 3 light -u, rf, s flavors. Since heavy quarks h = c, b, t will not be 
considered as 'light' partons in the photon (as in the case of the proton [26] and the pion 
[5,116]), their contributions to F^^^ ^ (i = 1,2, L) have to be calculated in fixed order 
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Figure 9.2: Direct (a) and resolved (h) contribution of a heavy quark h in deep inelastic 
e7(P^) scattering. 



perturbation theory. The heavy quark contribution consist of two parts, a 'direct' one 
and a 'resolved' one as depicted in Fig. 9.2: 



(9.12) 



The 'direct' contribution derives from the box-diagram 7*(Q^) 'j ^ hh expression, i.e. the 
usual Bethe-Heitler cross section [117] 
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(9.13) 



Q4 

where = 1 — Aml/W"^ = 1 — Amlx/ (1 — x)Q'^ . A similar expression holds for the longi- 
tudinal structure function [111] F^ = F2 — 2xFi. The 'resolved' heavy quark contribution 
to F^^^ ^ has to be calculated via 'y*{Q'^)g — > hh and is given by [111] 

dz 



F-p''-ix,Q'] 



-Z^7^^-^^(z,/i^)/2 



p7(P'),res 



2\ ii* iQ'^)g~^hh (X ^2 



z 



(9.14) 
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where I f^^^'^^^^~^ix,Q^) and 2ff^'^'^^^''~^{x,Q^) are given by Eq. (9.13) with e^a 
e1as{np)/Q, -^min = x{l + Am1/Q'^) and fip ~ Ami [18]^- To ease the calculations we 
shall keep these LO expressions in Eqs. (9.13) and (9.14) also in NLO, since the full 
NLO expressions for heavy quark production [119] turn out to be a small correction to 
the already not too sizeable (at most about 20%) contribution in LO. Notice that such 
small corrections are not larger than ambiguities due to different choices for rrih and for 
the factorization scale fip. For our purposes it is sufficient to include only the charm 
contributions which will be calculated using rric = 1.4 GeV. 

As already mentioned in the introduction the heavy quark contributions deserve special 
care due to the additional scale ruh provided by the heavy quark mass. For general 
P^{<^ Q^) terms P'^/m\ are not always power-suppressed and have to be kept in the 
heavy quark coefficient functions. The corresponding 'direct' and 'resolved' heavy quark 
contributions can be found in Eqs. (18) and (19) of Ref. [111]. 



9.3 Q2_Evolution 
9.3.1 Evolution Equations 

As motivated in Sec. 9.1 a transverse hadronic photon can be described by its partonic 
quark, gluon and photon content denoted by qj^^ \x,Q'^) (i = 1, . . . , /) with qi — u, 
q2 — d, Qs — s and so on^, g'^^^^\x, Q^), and V^^^^^x, Q^). The general evolution equations 
for these parton densities are given by (suppressing the x- and Q^-dependence) 



dg; 



ding 



2 



k=l 



g = P„ ® r^<-=' + 2 f: P,, » + P,, ® 9-"-) (9.15) 



fc=l 

^^''"^^ = p,, r-(-^) + 2 x: p.,. + p.. ^ 9^^''^ 



dlnQ2 



k=l 



^An alternative choice would be 1^% = + 4m^ [118] which satisfies the requirement /ti|, ^ also 
for large P^. 

^Note that due to charge conjugation invariance q^^^ \x, Q^) = q]^^ \x, Q^). 
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where / denotes the number of active (massless) quark flavors. The evolution kernels 
are generalized splitting functions 



P.,(x,a,a.)= j:--^pg-)(x), (9.16) 

l,m=0 ^ ' 

with Pq.q^ being the average of the quark-quark and antiquark-quark splitting functions. 
Note that terms due to photon radiation from the quark line are taken into account in 
Eq. (9.15) as it must be in order to guarantee momentum conservation (see the erratum 
to [120]). 

Since the electromagnetic coupling constant a ^ 1 we can safely neglect terms of 
order 0{a^) in Eq. (9.15)^. Noticing that q'^'^^''^ and g"'^^^^ are already of order 0{a) we 
can take / = in all generalized splitting functions in (9.16) which are multiplied by g'^'-^^^ 
or g"'^^^^ such that the functions Pij reduce in this case to the conventional QCD splitting 
functions Pij{x, a^). On the other hand in P^.-y, P^-y and P^^ the contributions with / = 
obviously vanish such that the only contribution to order 0{a) comes from setting I — 1 
in Eq. (9.16) in this case. Furthermore, the terms proportional to P-yg^. and P-yg are of 
the order 0{a'^) and can be dropped such that the evolution equation for the photon 
distribution decouples and can be solved separately as will be discussed in Sec. 9.3.5. 
Here it suffices to know the photon density at order q;° only in order to obtain the quark 
and gluon densities in Eq. (9.15) and the structure functions in Eq. (9.2) consistently to 
order 0{a). At 0th order, however, no splitting is possible and the photon carries its full 
momentum with probability one: V^^^^x, Q"^) — 5{1 — x) -\- 0{a). 

Therefore, Eq. (9.15) can be reduced to the following inhomogeneous evolution equa- 
tions for the quark and gluon content of the photon [120]: 

1 7(p2) / 
^ k=l 

where we have utilized the standard notation kg. = Pg^-y, k^ = P^,^ for the photon-quark 
and photon-gluon splitting functions. The splitting functions Pjj (,x, a^) are presently 
known to next-to-leading order in as and can be found in [112, 121, 122]. 
^For this reason we can also neglect the running of a{Q^) and use a ~ 1/137. 
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The Q^-evolution is most conveniently treated in the MeUin moment space where the 
convolutions turn into ordinary products: (a ® b){x) — > (a ® = a"6" where the n-th 
MeUin moment of a function f(x, Q'^) is defined by 



Jo 



(9.18) 



Hence the evolution of n-moments is governed by a set of ordinary coupled differential 
equations which can be solved analytically [112]. 

In order to decouple the evolution equations as far as possible we introduce fiavor 
non-singlet (NS) quark combinations: 



7(P2) 



,i=l 



y(p2) 



3<l3 



7(P2) 



(9.19) 



with i, J = 1, . . . , / and the group theoretical index I — j"^ — 1. Thus, for three active 
fiavors (m, d, s) we have two non-vanishing NS combinations Qns,3 = 2(m — d) and Q'ns.s ~ 
2{u-\-d — 2s). Furthermore, we define a singlet (S) vector which helps to write the singlet 
evolution equations and their solutions in a compact matrix form: 




,(P2) 



(9.20) 



i=l 



We wish to make an expansion of the solutions of the evolution equations and even- 
tually of the structure functions in terms of a^. For this reason the evolution equations 
will be rewritten as differential equations in ct^ with help of the renormalization group 
equation for the strong coupling constant: 



da. 



dlng2 



(9.21) 



with /3o = 11 — (2/3)/ and /3i = 102 — (38/3)/. Expanding the splitting functions 
in as in a form which can be generically written as P = §(P(°) + ^P^^^ + . . . ) and 
k — -^{k^^^ + ^k^^^ + ■■■)' respectively, and neglecting terms of the order 0{al) the 
inhomogeneous evolution equations can be brought into the following form: 



d ^7{P^),n _ _ a 2 (o)n 



a 



da 



ajpo 
J_Ap(o)n , 



as TTpo 



(l)n 



A , (0)n 

Wo''' 



n(3o 



p{l)n _ Pi p(0) 
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da, 



■5s 



(P2). 



Oil 

J_^p(0)n ^ J_ 



^J_f ^(l)n _ A g(0)« 

as ttPo V 2/3o 

p(l)n _ A p(0)n 

2/9o 



9s 



.7(p2),n 



(9.22) 



Here and denote the LO and NLO photon-to-parton splitting functions in the 
singlet sector: fc" = (/eg, k^)'^ with /eg = 2^f^j^ kg. while /c^°'" and /c^^^" are the LO and 
NLO parts of the non-singlet combination /c" = 2 Yli=i ~ J^q ■ Furthermore, p(°)" 



and p(^)" refer to LO and NLO 2x2 matrices of one- and two-loop splitting functions 
whereas pI^-*" and Pgg''" occurring in the NS equation arc scalar functions. The complete 
set of LO and NLO (non-)singlet splitting functions can be found, e.g., in [112, 123]. 



9.3.2 Analytic Solutions 



The general solutions of the inhomogencous evolution equations (9.22) are formally iden- 
tical to the ones of a real photon [112] and can be written, as usual, as a sum of a 
('pointlike') particular solution gpi of the inhomogencous problem and a general solution 
?had of the homogeneous ('hadronic') equations : 



_ 7(P^),n ^(P^),n ^(P^),n 

y ~ ^Ns,^ ' ys ' y ' 



(9.23) 



The particular solution is not unique and we fix it by the condition qpi{Ql) = at the 
reference scale Q^. 

The NLO 'pointhke' (inhomogencous) flavor-singlet solution is given by [111,112] 



27r 



1 + ^f/" 



2n 



1 



1 



a 



1 + d„ 



+ 



1-L 



dn 



1 a 



k^n _ A^(0)n _ ^n^(0)n 



(9.24) 



where 4 = -{2/ /3o)P^^>, L = as{Q'^)/as{Ql) and the 2 x 2 NLO evolution matrix 
U is fixed through the commutation relation [C/",P(°)"] = f [/" + P(^)" - [30, 
112, 123]. The NLO 'hadronic' (homogeneous) solution depending on input distributions 
^iia?'"(<5o) = Q^^^^^'^'iQl) which will be specified in Chapter 11 is given by [111, 112] 



2n 



/ r d„TTn 



2% 



yS,had Wo) 



(9.25) 
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The LO results are obviously entailed in these expressions by simply dropping all higher or- 
der terms (/3i, k^^^'^, U). The non-singlet solutions can be easily obtained from Eqs. (9.24) 
and (9.25) by the obvious replacements P^^^ — > P'^qq , P^^^ — > k ^ ke together with 

Note that the solutions (9.24) and (9.25) are valid for a fixed number f of active flavors. 
Evoluting beyond an MS 'threshold', set by the heavy quark masses nic, mh and mt, the 
number of active flavors has to be increased by one, i.e., / — > / + 1 in the calculation of 
as- The evolution of al {Q'^), corresponding to a number of / active flavors, is obtained 
by exactly solving in NLO(MS) 

numerically [26] using af\M^) = 0.114, rather than using the more conventional approx- 
imate solution 



a 



^^\Q^) 1 I3[^^ lnln(gVA2) 



s 



(9.27) 



47r ln(QVA2) {f3i^^f[HQV^')? 

which becomes sufficiently accurate only for ^ ~ 2 GeV^ with [26] A^'''^'^'' = 
257, 173.4, 68.1 MeV, whereas in LO {(3i = 0) A[Jf'^^ = 175, 132,66.5 MeV. Furthermore, 
P^/^ = 11 - 2//3 and (3[^'^ = 102 - 38//3. For matching a, at the MS 'thresholds' Q = 
Qf ~ ruf, i.e. aF^^^(m^^J = Q;F\m^_^J, we have used [26] nic = 1.4 GeV, mj, = 4.5 GeV 
and rrit = 175 GeV. On the other hand, we fix / = 3 in the splitting functions Pj^j'^^ in 
Eqs. (9.24) and (9.25) for consistency since we treat the heavy quark sector (c, b, . . .) by 
the perturbativcly stable full production cross sections in fixed-order perturbation theory, 
i.e. 7*(Q2)rT(^') ^ cc and Y{Q'^)9"'^^^^ ^ cc, etc., keeping ^ 0. 

Thus the full solutions can be obtained iteratively by the following procedure written 
out in detail in [124]: Between thresholds we can use Eqs. (9.24) and (9.25) for the 
corresponding value of /. Crossing a threshold Q"^ — ml we have to replace / by / -|- 1 in 
the following expressions occurring in Eqs. (9.24) and (9.25): 

and the full solution at the threshold evolved with / active flavors so far then serves 
as input at the new 'input' scale Qq = m\ which subsequently has to be evolved using 
Eqs. (9.24) and (9.25) with / + 1 active flavors. 
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9.3.3 Numerical Mellin— Inversion 

The solutions in x-space are then obtained by numerically inverting the Mellin transfor- 
mation in Eq. (9.18) according to 

fix, Q') = ^ l^dn x-y-{Q') (9.28) 

where the integration contour c has to enclose all poles in the complex n-plane. This 
can be realized, e.g., if the contour ranges from negatively to positively infinite imaginary 
values and lies to the right of all singularities of In the case of parton distributions 
and structure functions all singularities are located on the real axis and hence we can 
parameterize the contour as follows: 

c{z) = Co + e^'t'z , Co e R . (9.29) 

For such contours the transformation (9.28) can be written as [112] 

1 /"°^ 

f{x,Q^)^- dz Im {e'-^x-"^'^ r=<'\Q^)} . (9.30) 
Jo 

Employing > | the factor 3;-2^exp(i(/>) (dampens the integrand for increasing values of z 
allowing for a smaller upper limit z^s.x in the numerical calculation of (9.30) as compared 
to the standard integration contour with = ^. We choose the following parameters in 
all practical apphcations of Eq. (9.30) [112, 123]^: 

3 10 J 0.8 for non-singlet inversions 

4 ln(l/a;) M.8 for smglet mversions 

where the difference between non-singlet and singlet inversions is due to the different 
singularity structures of the corresponding splitting functions the rightmost pole lying at 
n = and n — 1, respectively. 



9.3.4 Boundary Conditions for a Deeply Virtual Target Photon 

For a deeply virtual target photon 7(-P^) with <^ <^ the photon structure func- 
tions are perturbatively calculable according to the doubly virtual box j*{Q^)j*{P^) — > qq 

"^The upper limit of integration Zmax refers to the photon case only. For inversions of pionic parton 
densities in Chap. 10 we have used Zmax{x < 0.001) = 24, Zmax(0.001 < a; < 0.05) = 40, 2:max(0.05 < a; < 
0.2) = 56, ^;max(0.2 < a; < 0.4) = 72, 2;max(0.4 < a; < 0.7) = 88 and 2:max(0.7 < a; < 1) = 136. 
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(and higher order corrections) either within fixed order perturbation theory (FOPT) or 
by resumming large coUinear logarithms InQ^/P^ occurring in the fixed order calculation 
with help of renormalization group (RG) techniques to leading order [125] and next-to- 
leading order [126, 127] accuracy. 

More specifically, in Ref. [126] the framework of the RG-improved operator product 
expansion (OPE) has been adopted in order to calculate predictions for the unpolarized 
(twist-2) structure functions F^^^^ ^^{x,Q'^) of a spin-averaged target photon. These 
(OPE) results can be translated by a one-to-one correspondence into the framework of 
the RG-improved parton model (PM) by utilizing 'technical' boundary conditions given 
at the input scale — ^ such that the condition 

FO™(a;, P^ Q^) = F<^(^')>'^"'(a;, Q') (9.31) 

is fulfilled for ~> P^. The appropriate boundary conditions are given by [111, 126, 127] 

/<j(^^>(x,g^ = P^) = (9.32) 

in LO and by 

P')-^^1P''{x,Q'-P') 

= Kel + (1 - ^f] In ^ + 6x(l -x)-2^ 

P^) = 0. (9.33) 

in the NLO(DIS^) factorization scheme. It should be noted, that these 'technical' bound- 
ary conditions have no direct physical interpretation at — P^ since the virtual target 
photon 7(P^) is not resolved by the scale = P^. Furthermore, at Q'^ — P^, the ar- 
range is kinematically restrained to x < 1/2 due to < x < (1 -l-P^/Q^)"^. However, the 
boundary conditions (9.32) and (9.33) generate the correct parton content at » P^ 
needed to reproduce the purely perturbative results of [126] in the Bjorken limit over the 
whole x-range < x < 1.^ 

It is instructive to derive Eqs. (9.32) and (9.33) from our fixed order calculation of the 

doubly virtual box 7* ((5^)7* (P^) qq in Chapter 8. The universal asymptotic (Q^ S> P^) 

^To be precise, Eq. (9.2) supplemented with the boundary conditions in (9.33) reproduces the OPE 
result of [126] up to spurious terms of the order 0{aas) which are of next-to-next-to-leading order 

(NNLO) in a perturbative expansion of the photon structure function Fj*^^^^ in powers of ag. 
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leading log (LL) part of the box to the structure function F^^^^ '^^{x, Q"^) in Eq. (C.IO) 
(or Eq. (C.8)) may be used to define light (anti-)quark distributions in the virtual photon 
target 

F^iP" i^, Q') luniv ^N,J2et^x[x' + {l- xf] In |! 

^ E ^4 [Cr'^'i^, Q') + CP'-i^. Q')] (9-34) 

q=u,d,s 

with 

CP'-i-^ Q') = Cr'^i^, Q') = N^el^ [x' + (1 - xf] In |^ . (9.35) 

Note that these fixed order box expressions do imply to this order a vanishing gluon 
component in the virtual photon, g^^^^ ''^ {x, Q^) — 0. 

Higher powers of the finite but asymptotically large logarithm InQ^/P^ occurring in 
(9.35) may be resummed to all orders in perturbation theory with help of the renormal- 
ization group by imposing leading order boundary conditions^ 

/<j(^^>(a;, = P') = /<r^>(x, Q' = P')=0, (/ = q, g) . (9.36) 

In order to derive the NLO input distributions in (9.33) we have to take into account 
the full asymptotic box structure function F^^^^ ^^(x, Q^) in Eq. (C.IO) following from 
Eq. (8.11) in the hmit <^ for the light quarks (m = 0) which can be used to define 
scheme-dependent light (anti-) quark distributions in the virtual photon target in NLO 

fSP'^ {X, Q') ^ F<->'^-'^^«^ {x,Q') + 2Y: elx CS> {x, Q^) (9.37) 

where the 'direct' term F^''^''^^^{x, Q"^) vanishes in the DIS-y scheme implying 

+ [x^ + (1 - xf] In ^ + 6a;(l - x) - 2^ . (9.38) 

The fixed order box' distributions in the DIS-y scheme can again be resummed employing 
the NLO(DIS-y) evolution equations with the input distributions 



Ncel ^[[x^+{l- xf] In ^ + Qx{l 



X] 



®This is quite similar to the case of heavy quark production studied in Part I of this thesis where 
logarithms In Q^/m^ of the heavy quark mass nih may be resummed to all orders using perturbatively 
calculable boundary conditions for the heavy quark distribution function in a hadron. 
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^Dit = P') = g^Sil: Q' = P')-o (9.39) 

which coincide with the results in Eq. (9.33). 

So far, the discussion has been for a spin-averaged target photon. However, the bound- 
ary conditions for a transverse target photon can be obtained for free from F^^^^ \x, Q'^) 
given in Eq. (C.8) by changing 6x(l — x) to 8x(l — x) in (9.39): 

= A^ceJ ^{ [x^ + (1 - x)^] In 1 + 8x(l - x) - 2} 

9'X'M" = P") = ^- (9-40) 

The corresponding results in the MS scheme can either be obtained by the factorization 
scheme transformation in (9.9) or by repeating the above described steps in the MS scheme 
with the 'direct' term 

^<7>,d^r,MS(^^ ^N,Y,4 [x^ + (1 - xf] In + 8x(l - x) - l| . (9.41) 

Finally, let us stress that the finite coUinear logarithm In Q^/ a priori needs not to 
be resummed. This issue has to be decided by phenomenological relevance. In addition, 
the predictive power of observables obtained from perturbatively calculated boundary 
conditions generally appears to be weakened by the freedom in choosing the input scale 
Qq [65] (at least at low orders in perturbation theory). Albeit Qq = is a 'natural' 
choice for the input scale due to the vanishing of the coUinear logarithm InQ^/P^ it is 
by no means compelling and a different input scale Qq P^ would result in equally well 
justified boundary conditions 

^ro''^'Ml) = N.el^[x- + {l-xf]\nf^ 

iZV^^^ Ql) = M + (1 - ^)'] 1- % 

+ [x' + (1 - xf] In ^ + 8x(l -x)-2^ 

gVSkMl)-^ (9-42) 

as long as IuQq/P^ is not large. Writing InQ'^/P^ — InQl/P^ + InQ'^/Ql we see in this 
case that a part of the coUinear logarithm is kept in fixed order (In Ql/ P^) while the other 
part (InQ^/Qo) resummed to all orders. 
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9.3.5 Momentum Sum 

Due to energy-momentum conservation the momentum fractions carried by the individual 
partons must add up to one^: 

-i 



dx X 



J2 r^'''\x,Q') + r'r(^'\x,Q') 

.f=g,q,9 



= 1 . (9.43) 



However, as we have seen before r^^^^\x, Q^) drops out of the general evolution equations 
in (9.15) leaving inhomogeneous evolution equations, see (9.17) and (9.22), for the (re- 
solved) partonic content of the photon and their solutions require boundary conditions at 
some reference scale (input scale) Q^. A sum rule for the total momentum of the quarks 
and gluons in the photon 

M2^(^'Hg2)= J2 f dx xP^^'\x,Q^) (9.44) 

would provide useful information for constraining the input distributions. Eqs. (9.43) and 
(9.44) can be rewritten in terms of Mellin moments (suppressing the obvious Q^- and 
P^~-dependence) 

To order a the evolution equation for V^^^^^ reduces to 

^L_^P„«n(-' (9.46) 

where the photon-photon splitting function P^^ can be constructed from the second mo- 
ments of ks = 2 ^ kg and kg as follows: Using the evolution equations (9.17) for the quark 
and gluon densities and (9.46) for the photon and exploiting in addition the conservation 
of the hadronic energy-momentum tensor one can write 
d 

dhTg^ 

Furthermore, P-y-y oc S{1 — x) to all orders in as because in order 0{a) only virtual 
diagrams contribute to photon-photon splittings while radiation of real photons starts at 
order 0{a'^) and we can write 

P^^ = 5(1 - x)F^=' = -5(1 - x)(kr' + K=^) ■ (9.48) 



''This is certainly true for a transversely polarized target photon considered here whereas the longitu- 
dinal degrees of freedom deserve possibly some attention. 
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The n-moments of the photon-gluon and photon-quark sphtting functions can be found, 
e.g., in [112, 123]. The second (n = 2) moments then read 

j^(0)n=2 ^ ^o)n=2 ^ j^(0)n=2 ^ 2 J] 6^, kg)"=' + k^^^^=' = 4 ^ (9.49) 

where the latter equation is vahd both in the MS and the DIS-y scheme as follows from 
Eq. (9.10) with + P(?"=' = 0. 

RecaUing the expansion k = ^(k^°^ + ^k*^^) + . . . ) we arrive at the following expression 
for the photon-photon sphtting function (to order a)^: 



p 



--S{l-x)^J24 \l + ^ + 0{al)] . (9.50) 

TT ^ — ' ^ L TT J 



Eq. (9.46) can be solved analytically either in n-moment space or directly in x-space 
since the convolution becomes trivial due to the 5(1 — x) in (9.50) and one easily finds 



r{x, Q') = r{x, Ql) exp [-^124 (l^ §1 + ) 



5{l-x) 



V 9 q 



+ 0{a^) (9.51) 



with ci being a constant (depending on and P^) which automatically enters the descrip- 
tion due to the boundary condition r(a;, Ql) = S{l-x) [l + 0{a)] and h = J ^4?^ In Q'^. 
Employing the renormalization group equation for the strong coupling in (9.21) Ii can be 
trivially integrated and one obtains 

' Po asm f^o 47r/5o + A«.(Q^) ' ^ ' ^ 

where the first logarithm in (9.52) dominates over the second one which is approximately 
'^^{'^siQ'^) — C(s{Qo)) and hence parametrically of NNLO. Note that contributions pro- 
portional to asiQ"^) — as{Qo) are generated also from terms of the order 0{al) in (9.50). 

Eq. (9.43) holds order by order in a, so that (9.51) imphes for the total momentum 
carried by the (resolved) photonic partons in LO-QCD [114] 

M,-(^^)(Q^) = ^ [e 1- f + '^(Qo, P')] (9.53) 



^Alternatively, P^-y can be inferred from the Abelian {Ca = 0), diagonal (independent of n) part of 
Pgg with the appropriate replacement Tfi = ^ ^ Nc J2q ^q- 
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and in NLO(DIS^) 

4 V- 21 47r/3o + /3ia,(g^) 
Po^""' ^47r/3o + /3ia,(g2) 

where the sum Ylq always extends over the 'light' u, d, and s quarks in the scheme 
adopted here (with / = 3 fixed in the splitting functions) whereas Pq, 0i and depend 
on the number of active flavors (and therefore the heavy quark thresholds have to be 
taken into account by an iterative procedure as described in Sec. 9.3.2). 

The constant ci is tightly related to the momentum sum of the hadronic input dis- 
tributions at = Ql: M2(Q'^ = QD/a = ci/n. While ci is perturbatively calculable 
for large » from the boundary conditions in (9.32), (9.33) or (9.40) this is gener- 
ally not the case for A^, especially not for real P'^ — photons and an important 
constraint on the parton densities seems to be missing. 

In the past few years some attempts were undertaken [97, 128, 129] to infer Ci from 
elsewhere by relating Eq. (9.51) to the hadronic part of the photon vacuum polarization 
which in turn can be determined via a dispersion relation from the well measured cross 
section (Th = (T{e~^e~ hadrons). However, it turns out that the usefulness of such an 
approach is spoiled by higher-twist contributions present in the experimental quantity 
[129]. We will return to this issue when we discuss boundary conditions for the real 
photon in Chapter 11. 

9.4 Longitudinal Target Photons 

Before turning to our analysis of the parton content of pions within a constituent quark 
model needed as basic ingredient for our analysis of the parton content of real and trans- 
versely polarized virtual photons within the framework of the radiative parton model in 
Chapter 11 let us discuss for completeness also the case of longitudinally polarized virtual 
target photons. While these contributions vanish like, say, P^/m^ in the real photon limit 
due to gauge invariance they yield a finite contribution if P^ ^ A^ and should be taken 
into account in Eq. (7.52). 
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As already mentioned in Sec. 9.3.4, the structure functions of (transverse or longitu- 
dinal) deeply virtual target photons 7(-P^) with <^ <^ are reliably calculable in 
fixed order according to the doubly virtual box 7*(Q^)7*(-P^) qq- Nevertheless, it has 
been demonstrated recently [130, 131] that a partonic treatment of 7l(P^) is phenomeno- 
logically useful and relevant. The quark and gluon distribution functions q'^^^^^\x , Q^) 
and g'^^''^^\x, Q^) of a longitudinal taiget photon satisfy /lomo^'eneous evolution equations 
[132] in line with the expectation that (to order 0{a)) there is no (on-shell and there- 
fore) transverse photon-parton inside the longitudinal target photon, r^'^^^''\x , Q'^) = 0. 
In any case, the complete theoretical framework for the longitudinal photon target 7l is 
obtainable from Sec. 9.2 and 9.3 by formally setting = 0, i.e by dropping all 'direct' 
contributions (k^ (g) V^'^ = 0, T'^'^ ®C^ = etc) and by using the perturbatively calculable 
boundary conditions at 3> A^ [132] 



q-y-^^'Xx, = P') = ?^(^')(x, = P') = N4 ^Ax{1 - x) 

In 

g^-^^'\x,Q' = P') = Q, 



(9.55) 



which follow directly from the asymptotic box expression in Eq. (C.9). 



Chapter 10 

Pionic Parton Densities in a 
Constituent Quark Model 



Before specifying our boundary conditions for the parton distributions of real and virtual 
photons we derive pionic parton distributions from constituent quark model constraints. 
Being, of course, also interesting in their own rights, they will serve in the next chapter as 
input for the hadronic component of the photon. The analysis presented in this chapter 
is almost identical to Ref. [5]. 

The parton content of the pion is poorly known at present. The main experimental 
source about these distributions is mainly due to data of Drell-Yan dilepton production in 
7r~-tungsten reactions [133-137], which determine the shape of the pionic valence density 
f '^(.T, Q^) rather well, and due to measurements of direct photon production in n^p 'jX 
[133, 138-140] which constrain the pionic gluon distribution g'^{x, Q^) only in the large-a; 
region [141]. In general, however, present data are not sufficient for fixing g'^ uniquely, in 
particular the pionic sea density q^{x, Q^) remains entirely imconstraincd experimentally. 
Therefore in [116] a constituent quark model [142] has been utilized to relate q'^ and to 
the much better known radiatively generated parton distributions f^{x, Q^) of the proton 
[41]. Therefore, previously [116] a constituent quark model [142] has been utilized to 
relate q'^ and g'^ to the much better known radiatively generated parton distributions 
P{x,Q'^) of the proton [41]. These relations arise as follows: describing the constituent 
quark structure of the proton p — UUD and the pion, say n'^ — UD, by the scale (Q^) 
independent distributions C/^''^^(x), Dp{x) and D'^'^{x), and their universal (i.e. hadron 
independent) partonic content by Vc{x,Q^), gc{x,Q^) and qc{x,Q^), the usual parton 
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content of the proton and the pion is then given by 

nx,Q') = f^^[U^{y) + D^{y)]fJl,Q^) (10.1) 

FM') = f-\u''\y) + D-\y)\fc{-.Q') (10.2) 
J X y \y / 

where f = v, q, g with = + = {u^ + rfP)/2, t;^ = + (i^^, g'^ = 

{u'"^ + d'"'^)/2 and zZ'^^ = o?'^^ due to ignoring minor SU(2)flavor breaking effects in the 
pion 'sea' distributions. Assuming these relations to apply at the low resolution scale 
g2 ^ ^2 ^^2^ ^ 0.23 GeV2, 

A*NLO — 0-34 GeV^) of [41] where the strange quark content 
was considered to be negligible, 

(x, /x^) = sP{x, i/) = s'^ix, n^) = s'^ix, n^) = 0, (10.3) 

one obtains from (10.1) and (10.2) the constituent quark independent relations [116] 

T;^(n,/x^) _ g"(n,/i^) _ ^"(n,/i^) 
v'P{n,p?) qP{n,iJ?) gP{n,iJ?) 

where for convenience we have taken the Mellin n-moments of Eqs. (10.1) and (10.2), i.e. 
f{n,Q^) = x"^"^ f{x,Q'^)dx. Thus, as soon as v^{x,fi'^) is reasonably well determined 
from experiment, our basic relations (10.4) uniquely fix the gluon and sea densities of the 
pion in terms of the rather well known parton distributions of the proton: 

a'^in,!.') = ^3^/(n,/.^), r{n,l^') ^ ^-^^^^ q^n, (10.5) 
Furthermore, the sum rules [116] 



"1 



I 



[ v''{x,Q'^)dx = 2 (10.6) 

1 /•! 



xv''{x,Q^)dx = f xvP{x,Q^)dx (10.7) 
Jo 



impose strong constraints on v'^{x,ii'^) which are very useful for its almost unambiguous 
determination from the ttN Drell-Yan data. Notice that Eq. (10.7), together with (10.4), 
implies the energy-momentum sum rule for to be manifestly satisfied. In addition, 
Eq. (10.7) implies that the valence quarks in the proton and the pion carry similar total 
fractional momentum as suggested by independent analyses within the framework of the 
radiative parton model [41, 141]. 
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The relations in Eq. (10.5) imply that any updating of f^{x, ji'^) yields a corresponding 
updating of /'^(x, /x^). Recently an updating of fP{x,fj,'^) within the framework of the 
radiative (dynamical) parton model was undertaken [26] utilizing additional improved 
data on F^{x,Q^) from HERA [143-148] and a somewhat increased as{Mz) = 0.114 
resulting in a slight increase in /i^ (fJ-to = O-^S GeV^, //-nlo ~ ^-^^ GeV^). An improved 
treatment of the running as{Q'^) at low was furthermore implemented by solving in 
NLO(MS) Eq. (9.26) numerically [26] rather than using the approximate NLO solution 
in Eq. (9.27) as done in [41, 116, 141], which is sufficiently accurate only for >, ~ 2 
GeV^ [26]. The LO and NLO evolutions of f'^{n,Q'^) to Q"^ > /x^ are performed in 
Mellin n-momcnt space, followed by a straightforward numerical Mellin-inversion [92] to 
Bjorken-x space as described in Chapter 9. It should be noted that the evolutions are 
always performed in the fixed (light) / = 3 fiavor factorization scheme [18, 26, 41, 116], i.e. 
we refrain from generating radiatively massless 'heavy' quark densities h'^{x,Q'^) where 
h = c, b, etc., in contrast to [141]. Hence heavy quark contributions have to be calculated 
in fixed order perturbation theory via, e.g., g^g^ hh, W^u^ hh, etc. (Nevertheless, 
rough estimates of 'heavy' quark effects, valid to within a factor of 2, say, can be easier 
obtained with the help of the massless densities c^{x, Q^) and b'^{x, Q^) given in [141].) 

Using all these modified ingredients together with the new updated [26] f^{x, ji"^) in 
our basic predictions in Eq. (10.5), the present reanalysis of the available Drell-Yan data 
[134-136], closely following the procedure described in [116], yields 

<o(a;>/"Lo) = 1.129a;-o-49^(l-a;)0-349(l + 0.153V^) (10.8) 
^Slo(^,/^nlo) = 1.391x-°-^^^(l-x)°-^26 (10.9) 

where [26] ji^Q = 0.26 GcV^ and /iNLo ~ ^.40 GeV^. These updated input valence 
densities correspond to total momentum fractions 

[ xvlQ{x,filQ)dx = 0.563 (10.10) 
Jo 

f xv^^o{^,pi^^^o)dx = 0.559 (10.11) 
Jo 

as dictated by the valence densities of the proton [26] via Eq. (10.7). 

Our new updated input distributions in Eqs. (10.8), (10.9) and (10.5) are rather differ- 
ent than the original GRV^ input [141] in Fig. 10.1 which is mainly due to the vanishing 
sea input of GRVtt in contrast to the present one in Eq. (10.5). On the other hand, our 
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Figure 10.1: The valence and valence-like input distributions xf'^{x,Q'^ — /x^) with / = 
V, q, g as compared to those of GRV^r [141]. Notice that GRV^ employs a vanishing SU(3)flavor 
symmetric q'^ input at /Xlo = 0-25 GeV^ and /^nlo = 0-3 GeV^ [141]. Our present SU(3)flavor 
broken sea densities refer to a vanishing s'" input in (10.3), as for GRVtt [141]. 
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Figure 10.2: Comparison of our NLO valence distribution at = 20 GeV^ with the one 
of GRV^ [141] and GRS [116]. This density plays the dominant role for describing presently 
available ttN Drell-Yan dimuon production data. For illustration, the gluon and sea densities 
are shown as well. The SU(3)flavor symmetric GRV^ sea q'^ — is not shown, since it is 
similar to s"^ of our present analysis and of GRS which are all generated from a vanishing input 
at Q'^ = jj?, cf. Eq. (10.3). The SMRS [137] results refer also to a SU(3)flavor symmetric sea 
q"" = u'"^ = d'"^ = s'^ = s'". 
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updated input in Fig. 10.1 is, as expected, rather similar to the one of [116]. In both cases, 
however, the valence and gluon distributions become practically indistinguishable from 
our present updated ones at scales relevant for present Drell-Yan dimuon and direct^ 
production data, = — 20 GeV^, as illustrated in Fig. 10.2. Therefore our 

present updated pionic distributions give an equally good description of all available t:N 
Drell-Yan data as the ones shown in [116]. Notice that the different gluon distributions 
presented in Fig. 10.2 can not be discriminated by present direct-photon production data 
[138-140] due to the uncertainty of the theoretically calculated cross section arising from 
variations of the chosen factorization scale and from possible intrinsic kx contributions, 
cf. for example L. Apanasevich et al. [138-140]. 

For completeness let us mention that our basic predictions (10.5) for the valence-like 
gluon and sea densities at = /x^, as shown in Fig. 10.1, can be simply parametrized in 
Bjorken-x space : in LO at = /^lq — 0.26 GeV^ 

xg^'ix^filo) = 7.326x^-^^^(1 - 1.919 1.524 x)(l -x)^-^^^ 
xq''{x,iJ,lo) = 0.522 a;°-^'^°(l - 3.243 + 5.206 a;) (l-x)^-^°, (10.12) 

whereas in NLO at — //nlo — 

0.40 GeV^ we get 

a;c/"(x,^^Lo) = 5.90 a;i-2™(l - 2.074 + 1.824a;)(l-a;)^-^^° 
a;g^(x,/x^Lo) = 0.417 x°-2°^(l - 2.466 + 3.855 a;)(l - a;)^-^^^ (10.13) 

Finally, Fig. 10.3 shows our resulting predictions for xg'^{x,Q^) and xq'^{x,Q'^) as 
compared to the former GRVjr results [141]. The GRV^ results for xq'^ are significantly 
steeper and softer for x >, 0.01 due to the vanishing SU(3)flavor symmetric (light) sea input 
X q'^{x, /J,"^) = 0, in contrast to our present approach [116] based on a more realistic finite 
light sea input in Eq. (10.5). The valence-like gluon and sea inputs at = fj,'^. which 
become (vanishingly) small at a; < 10~^, are also shown in Fig. 10.3. This illustrates again 
the purely dynamical origin of the small-x structure of gluon and sea quark densities at 

> /i^. Our predictions for s'^ — s'^, as evolved from the vanishing input in Eq. (10.3), 
are not shown in the figure since they practically coincide with q'^{x, Q^) of GRVtt shown in 
Fig. 10.3 which also results from a vanishing input [141]. Simple analytic parametrizations 
of our LO and NLO predictions for p{x, Q'^) are given in Appendix D. 

To conclude let us recall that an improvement of /'^(x, Q'^) is particularly important 
in view of its central role in the construction of the photon structure function and the 
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Figure 10.3: The small-x predictions of our radiatively generated pionic gluon and sea-quark 
distributions in LO and NLO at various fixed values of Q'^ as compared to those of GRV^ 
[141]. The valence-like inputs, according to Eq. (10.5) as presented in Fig. 10.1, are shown 
for illustration by the lowest curves referring to /x^. The predictions for the strange sea density 
_ g|,g similar to the GRV^ results for q'^ . The results are multiplied by the numbers 
indicated in brackets. 

photonic parton distributions [96,97,112-114,149]. Furthermore, recent (large rapidity 
gap) measurements of leading proton and neutron production in deep inelastic scattering 
at HERA [150] allow, under certain (diffractive) model assumptions, to constrain and test 
the pion structure functions for the first time at far smaller vales of x (down to about 
10~^) than those attained from fixed target TriV experiments. 



Chapter 11 



Radiatively Generated Parton 
Distributions of Real and Virtual 
Photons 

Having presented the general theoretical framework in Chapter 9 we now turn to an anal- 
ysis of the partonic content of (virtual) photons within the phenomenologically successful 
framework of the radiative parton model [26, 41, 57, 92, 93]. This analysis follows Ref. [6]. 
The results in Sec. 11.3 have been taken from [7]. 

11.1 Introduction 

Modern theoretical QCD studies [96, 112, 113, 149, 151] of the parton distributions of real, 
i.e. on-shell, photons /'''(x, Q^), f = q, q, g, agree surprisingly well with measurements 
of the (anti-) quark and gluon contents of the resolved real photons as obtained from 
e^e^ and ep reactions at collider energies (for recent reviews, see [89,114,152]). For 
clarity let us denote the resolved real target photon with virtuality = —p^ ~ by 
7 = 7(P^ ~ 0) which is probed by the virtual probe photon 7*((5^), = —q"^, via the 
subprocess 7*(Q^)7 ^ X as in e~^e~ e^e~X. Here, p denotes the four-momentum of 
the photon emitted from, say, an electron in an e+e^ or ep collider. In the latter case it is 
common to use instead of for denoting the photon's virtuality, but we prefer P^ for 
the subprocess ^{P'^)p X according to the original notation used in e^e^ annihilations. 
(Thus the factorization scale in f^^^^\x,Q^) refers now to some properly chosen scale of 
the produced hadronic system X, e.g. Q ~ p^^* in high-pT jet events, etc.). 
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In general one expects [97, 104, 111, 125-127, 153-155] also a virtual photon 7(P^ ^ 0) 
to possess a parton content f^^^^\x, Q"^). It is a major problem to formulate a consistent 
set of boundary conditions which allow for a calculation of p^^'^^x, Q"^) also in the next- 
to-leading order (NLO) of QCD as well as for a smooth transition to — 0, i.e. to the 
parton distributions of a real photon (see Refs. [97,111,156] for a detailed discussion). 
Indeed, experimental studies of the transition of the deep inelastic (di-)jet cross section 
from the real photon to the virtual photon region at HERA point to the existence of a 
nonvanishing, though suppressed, parton content for virtual photons [157-161]. These 
measurements have triggered various analyses of the dependence of the ep jet production 
cross section on the virtuality of the exchanged photon [162, 163] and experimental tests of 
such predictions will elucidate the so far unanswered question as to when a deep inelastic 
scattering (DIS) ep process is eventually dominated by the usual 'direct' 7* = 7(-P^) 
induced cross sections, not contaminated by the so far poorly known resolved virtual 
photon contributions. More recently, NLO calculations of the (di-)jet rate in ep (and 
67) scattering, which properly include the contributions of resolved virtual photons, have 
become available [164] and the resolved virtual photon contributions have already been 
included in a Monte Carlo event generator [165] as well. 

It is the main purpose of the present chapter to formulate a consistent set of boundary 
conditions, utilizing valence-like input parton distributions at the universal target-mass 
independent [41, 57, 93, 96, 112, 141] low resolution scale Qo = /^^ — 0.3 GeV^, which allow 
for a calculation of p^^'^\x, Q^) also in NLO-QCD as well as for a smooth transition to 
the parton distributions of a real photon, = 0. We shall furthermore employ the pionic 
parton distributions of Chapter 10, p{x, Q^), which are required for describing, via vector 
meson dominance (VMD), the hadronic components of the photon. It should be recalled 
that the pionic gluon and sea densities, g'^{x,Q^) and q'^{x,Q^), have been uniquely 
derived from the experimentally rather well known pionic valence density v'^{x,Q'^) and 
the (recently updated [26] dynamical) parton distributions f{x, Q^) of the proton. Thus 
we arrive at essentially parameter-/ree predictions for p^^^\x, Q^) which are furthermore 
in good agreement with all present measurements of the structure function of real photons. 

Sec. 11.2 is devoted to an analysis of the parton distributions and structure functions 
of real photons and the resulting predictions are compared with recent experiments. In 
Sec. 11.3 we compare our unique small-a; predictions -characteristic for the radiative 
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parton model- with very recent small-x measurements of the photon structure function 
F2{x,Q'^) by the LEP-OPAL collaboration. Sec. 11.4 contains the formulation of our 
model for the parton distributions of virtual photons, together with some quantitative 
predictions for structure functions as well as a comparison with very recent data extracted 
from DIS dijet events. Our conclusions are drawn in Sec. 11.5. In App. D we present simple 
analytic parametrizations of our LO and NLO predictions for the parton distributions of 
real and virtual photons. 

11.2 The Parton Content of Real Photons 
11.2.1 Boundary Conditions 

The picnic parton distributions presented in the previous chapter will now be utilized 
to construct LO and NLO input distributions for the real photon via a vector meson 
dominance ansatz at the low resolution scale Qq = /^^ ~ 0.3 GeV^. 

In NLO the nonperturbative hadronic VMD input P{x, Ql) refers to the partons in 
the DIS^ scheme which guarantees the perturbative stability of the resulting F^{x, Q'^) 
provided this input is given by the NLO f'^{x, Ql) of Chap. 10, while the corresponding 
input in LO is given via VMD by the LO f'^(x, Ql) of Chap. 10. We shall assume that 
the input resolution scale Ql — ij!^ c:^ 0.3 GeV^ for the valence-hke parton structure is 
universal, i.e. independent of the mass of the considered targets p, tt, 7, etc. [41,57,93, 
96,112,116,141]. 

The hadronic VMD ansatz for p{x, jj?) is based on a coherent superposition of vector 
mesons [97] 

|7)M^had - y \p)^fi + ^ (11-1) 

where the 0-meson contribution is considered to be strongly suppressed at /i^ ^ m|. As 
usual, we identify the parton distribution functions in the parton model with matrix ele- 
ments of local twist-2 operators within the operator product expansion (OPE) formalism 
and we obtain using (11.1) 

= fUf^') = (7|0/|7),^had 

= 7^(p|C>/Ip)m2 + ——{{p\0f\u;)^2 + (u;|0/|p)^2) + — (a;|0/|u;)^2 

Jp JpJijJ Jul 
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with f — u,d, s, g. Assuming, within the OPE, 



{P\0q\p)^fl 




\0q\uj)i^2 


= (7r° 




(p|0,|a;)^2 = 2 73,e-^^(7r°|0,|7r°)^2 


(p|0sIp)m2 




\0g\uj)^2 






{p\Og\u) = 



where the last equation holds due to isospin conservation, one obtains 

{u + uy{x,fj,'^) = a{gl + gl + 2 g^g^ cos 9){u + uY\x,ix'^) 

{d + dy{x,iJ,^) = a{gl + gl- 2 gpg^ cos 9){d + dy\x, 11^) 

{s + sy{x,t,') = a{gl + gl){s + sr\x,t,')^0 

g-^ix^f,') = a{gl + gl)g^\x,pi'). (11.2) 

Here Oq^g refer to the leading twist-2 quark and gluon operators and gy = Att/ fy with 

gl = 0.50, gl = 0M3, (11.3) 

i.e. /p/47r = 2.0 and /J/47r = 23.26, as obtained from a zero-width calculation of the 
relevant leptonic widths r{V — > i'^£~) — a'^mygy/S presented in [166]. The omission of a 
finite-width correction for g^ is due to the central role [166] of the precise results in [167] 
which do not require such a correction in contrast to the situation for the less precise 
resonance analysis at e+e~ colliders [85, 168, 169]. 

For the a priori unknown coherence factor (fit parameter) cos^ in Eq. (11.2) wc take 
cos^ = 1, i.e. we favor a superposition of u and d quarks which maximally enhances 
the contributions of the up-quarks to in Eq. (9.1). This favored value for cos^ is 
also supported by fitting cos^^ in (11.2) to all presently available data on F^(a;, Q^), to 
be discussed below, which always resulted in cos^^ ~ 1 in LO as well as NLO. This is 
also in agreement with the LO results obtained in Ref. [97]. The LO and NLO input 
distributions /'^(x, /x^) of the pion in (11.2) arc taken from the analysis in the previous 
chapter which correspond to [5, 26] /Xlo — 0.26 GcV^ and /U^lq = 0-40 GcV^ in LO and 
NLO, respectively. Since by now all free input quantities have been fixed in Eq. (11.2), 
we arrive at rather unique parameter-/ree predictions for p{x, Q'^) and ^2(0;, Q"^). 

The calculation of p{x, Q"^) at > //^ follows from the well known inhomogeneous 
renormalization group (RG) evolution equations in LO and NLO which we solve, as 
usual, analytically for the n-th MeUin moment of p{x, Q^), followed by a straightforward 
MeUin-inversion to Bjorken-x space as described in detail in Chapter 9. The general 
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structure of these solutions is 

Fix, Q') = Q') + fUx, Q') . (11.4) 

Here /Jj denotes the perturbative 'pointhke' solution which vanishes at = fi^ and is 
driven by the pointlike photon splitting functions kf_g\x) appearing in the inhomogeneous 
evolution equations, while depends on the hadronic input p^x^^ii^) in Eq. (11.2) and 
evolves according to the standard homogeneous evolution equations. 

The prescription for the VMD ansatz in Eq. (11.2) at the input scale /i^, together with 
cos ^ = 1 as discussed above, yields a simple expression for the general Q^-dependence of 

r(x,g2): 



r(x,g2) = /;,(x,Q2) + a 



G)nx,Q') + 5f\{Gl-Gl)s%x,Q-') 



(11.5) 



with 5u — —1, 5d — +1 and 5s — — 0, and where the index tt obviously refers to 7r° and 

Gl = (^7p + ^7j'^ 0.836 
Gl = ((?,- ~ 0.250 

Gl = Gl^gl + gl^QM?,. (11.6) 

Simple analytic LO and NL0(D1S^,) parametrizations for the pointlike piece /Jj(a;, Q^) 
are given in App. D.2, whereas the ones for /'^(x, Q"^) can be found in App. D.l. 



11.2.2 Quantitative Results 

Having outlined the theoretical basis for our photonic parton distributions, we now turn 
to the quantitative results. First we apply our parameter-free predictions for f'^{x^Q'^) 
to the structure function of real photons which, according to Eq. (9.1), is finally given by 

If:^{x,Q') = 2 elL\x,Q')+'^^[C,,,®q' + C,,,®g<]\ 

q=u,d,s 

+ -FUx,Q')+-FC{x,Q') (11.7) 

X ' X ' 

where p{x, Q"^) refers to the DIS^ factorization scheme defined in (9.9) and the charm 
contributions and are given by Eqs. (9.13) and (9.14), respectively. 
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Figure 11.1: Comparison of our radiatively generated LO and NLO(DIS-y) predictions for 
F2 (x, (5^), based on the valence-like parameter-free VMD input in Eq. (11. 2), with the data 
of Ref. [170-179]. For our comparison the GRV-y [96, 112] results are shown as well. In both 
cases, the charm contribution has been added, in the relevant kinematic region W > 2mc, 
according to Eqs. (9.13) and (9.14). 



11.2 The Parton Content of Real Photons 



115 




Figure 11.2: Comparison of our predicted LO and NLO(DIS-y) distributions m''^ = u^, cP = d'^ 
and g'^ at = 10 GeV^ with the LO/NLO GRV^ densities [96, 112], the LO SaS ID and 2D 
[97] and the NLO AFG [113] distributions. The 'hadronic' (pionic) components of our total 
LO and NLO results in Eq. (11.4) are displayed by the dashed curves. 
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In Fig. 11.1 we compare our LO and NLO predictions with all available^ relevant data 
[170-179] for F2 of the real photon. Our present new NLO results are rather similar to 
the ones of AFG [113], but differ from the previous (GRV-y) predictions [96,112] which 
are steeper in the small-x region, as shown in Fig. 11.1, because the dominant hadronic 
(pionic) sea density q'"{x,Q'^) is steeper since it has been generated purely dynamically 
from a vanishing input at Q"^ = //^ [96, 112, 141]. Similarly the SaS ID [97] expectations, 
for example, fall systematically below the data in the small to medium region around 

~ 5 GeV^, partly due to a somewhat different treatment of the hadronic coherent 
VMD input as compared to our results in Eqs. (11.2), (11.5), and (11.6). 

The relevant LO and NLO photonic parton densities are compared in Fig. 11.2 at — 
10 GeV^. For illustration we also show the purely 'hadronic' component (homogeneous 
solution) in (11.4) of f"' which demonstrates the dominance of the 'pointlike' component 
(inhomogeneous solution) in (11.4) for u'^ and (P in the large-x region, x > 0.1 . 

In Fig. 11.3 we show our predictions for xu^{x,Q'^) and xg'^{x,Q'^). The parton dis- 
tributions of the photon behave, in contrast to the ones of a hadron, very differently in 
the limits of large and small x. In the former case, the purely perturbative pointlike 
part in (11.4) dominates for x ^ 0.1, especially for the quark distributions. On the other 
hand, this uniquely calculable contribution amounts at most to about 20% at very small 
X where the hadronic VMD component in (11.4) dominates, giving rise to a very sim- 
ilar increase for x — > as observed in the proton case. In Fig. 11.3 we also show our 
valence-hke inputs at = /Xlo,nlo which become (vanishingly) small at x < 10"^. This 
illustrates the purely dynamical origin of the small-x increase at Q'^ > /x^. Also note- 
worthy is the perturbative LO/NLO stability of u'~'{x, Q"^) which is almost as good as the 
one required for a physical quantity like F2{x, Q^) in Fig. 11.1. The situation is, as usual 
[26,41,57,93,96, 112], different for g"'{x,Q'^). Nevertheless, despite the sizable difference 
between the LO and NLO gluon distributions in Fig. 11.3 in the small-,x region, the di- 
rectly measurable and the gluon-dominated heavy quark contribution in Eq. (9.14) 
shows a remarkable perturbative stability [26]. 

Next, we compare in Fig. 11.4 our predictions for xg'^{x, Q"^) at = (p^^*)^ = 75 GeV^ 
with recent HERA (HI) measurements [180,181]. Our somewhat flatter results for xg'^ 
in the small-a; region, as compared to the older GRV-y expectations [96, 112], is caused by 
the recently favored flatter gluon distribution in the proton [26] which determines g'^ via 
^January 1999 



11.2 The Parton Content of Real Photons 



117 



100 = (GeV^) X g\x,Q^) la 




Figure 11.3: Detailed small-x (as well as large-x) behavior and predictions of our radiatively 
generated = and g'^ distributions in LO and NLO(DIS^) at fixed values of Q^. The 
dashed-dotted curves show the hadronic NLO contribution f^^^ to /'"' = + f^^^ in Eq. 
(11.4). The valence-like inputs at = A*lo,nlO' according to Eq. (11.2), are shown by 
the lowest curves referring to jj?'. For comparison we show the steeper NLO GRV,, [96,112] 
expectations as well. The results have been multiplied by the number indicated in brackets. 



g"" [5], cf. Eq. (11.2), at small x. 

Finally it is interesting to consider the total momenta carried by the photonic partons, 



M^^iQ^)^ V [\r{x,Q')dx. 

r _ Jo 



;ii.8) 



f=g,g,9 



Inspired by the ideas and suggestions put forward in Refs. [97, 128], it has been conjectured 
recently [129] that this leading twist-2 quantity M2 should satisfy, in LO-QCD, 



(11.9) 



where the well known dispersion relation relates the hadronic part of the photon's vacuum 
polarization 

-i2 roo 



(11.10) 
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Figure 11.4: Comparison of our LO and NLO predictions for xg'^ at = {{p'x)'^) = 75 
GeV^ with HERA(Hl) measurements [180, 181]. The GRV^ and SaS expectations are taken 
from Refs. [96, 112] and [97], respectively. 

to ah = a{e'^e'' — > hadrons). It should be noted that Ilh{Q^), being an experimental 
quantity, includes, besides the usual twist-2 term, all possible nonperturbative higher- 
twist contributions. The 'consistency' relation (11.9) is, however, expected to hold already 
at ^ 2 to 4 GeV^ to within, say, 20 to 30% where the twist-2 component in Ilh{Q'^) may 
become dominant, as possibly indicated by DIS ep processes. Indeed, our LO results imply 
M^{2 GeV^)/a ~ 0.976 and M^{4: GeV'^)/a ~ 1.123 which compares favorably with [182] 
Uh{2 GeV^)/a = 0.694 ± 0.028 and 0^(4 GeV^)/a = 0.894 ± 0.036, respectively. 

11.3 The Photon Structure Function at Small— x 

Recently the OPAL collaboration [98] at the CERN-LEP collider has extended the mea- 
surements of the photon structure function F2{x,Q'^) into the small-a; region down to 
X ~ 10~^, probing lower values of x than ever before. The observed rise of towards 
low values oi x, x < 0.1, is in agreement with general QCD renormalization group im- 
proved expectations. It has, however, been noted that the rising small-x data at lower 
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scales ~ 2 — 4 GeV^ lie above the original QCD expectations anticipated almost a 
decade ago [96, 97]. 

In the following we will demonstrate that our updated paramcter-/ree QCD predictions 
for F2{x,Q'^) [6], discussed in Sec. 11.2, are in general also consistent with the OPAL 
small-x measurements at all presently accessible values of Q^. 

Before presenting our results it is instructive to recapitulate briefly the main differences 
between the original GRV-y [96] approach to the photonic parton distributions and our 
parameter-free predictions of Sec. 11.2. In the latter approach a coherent superposition of 
vector mesons has been employed, which maximally enhances the li-quark contributions 
to F^, for determining the hadronic parton input f^{x,Ql) at a GRV-hke [26] input 
scale Ql = /^lo — 0-26 GeV^ and Ql = /x^lo — O-^O GeV^ for calculating the (anti)quark 
and gluon distributions p{x, Q"^) of a real photon in leading order (LO) and next-to- 
LO (NLO) of QCD. Furthermore, in order to remove the ambiguity of the hadronic 
light quark sea and gluon input distributions of the photon (being related to the ones of 
the pion, P{x,Qq), via vector meson dominance), inherent to the older GRV-y [96] and 
SaS [97] parametrizations, predictions [5,116] for q'^{x,Q'^) and g'^{x,Q'^), cf. Chapter 
10, have been used in our present approach which follow from constituent quark model 
constraints [94,95]. These latter constraints allow to express q'^ and g'^ entirely in terms 
of the experimentally known pionic valence density and the rather well known quark- 
sea and gluon distributions of the nucleon, using most recent updated valence-like input 
parton densities of the nucleon. Since more recent DIS small-x measurements at HERA 
imply somewhat less steep sea and gluon distributions of the proton [26], the structure 
functions of the photon will therefore also rise less steeply in x than the previous GRV^ 
[96] ones as was already discussed in the previous section and as will be seen in the figures 
shown below. In this way one arrives at truly parameter-free predictions for the structure 
functions and parton distributions of the photon. 

In Figs. 11.5 and 11.6 we compare our predictions of Section 11.2, denoted by GRSc 

[6] and the older GRV-y results [96] with the recent small-x OPAL measurements [98] and, 

for completeness, some relevant L3 data [179] are shown as well. Our parameter-free LO 

and NLO expectations are confirmed by the small-x OPAL data^ at all (small and large) 

experimentally accessible scales Q^. This is in contrast to the GRV-y and SaS results 

^Note that the new OPAL data [98] at = i g 3 7 QeV^ supersede the older OPAL data shown 
in Fig. U.l at Q2 = i.ge and 3.76 GeV^, see [98]. 
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which at LO are somewhat below the data at small in Fig. 11.5 and seem to increase 
too strongly at small x in NLO, in particular at larger values of as shown in Fig. 
11.6. The main reason for this latter stronger and steeper x-dependence in LO and NLO 
derives from the assumed vanishing (picnic) quark-sea input at Qq = /^lo.nlo 
anti(quark) distributions of the photon as well as from relating the hadronic gluon input 
of the photon directly to its (pionic) valence distribution [96, 141]. This is in contrast to 
the more realistic (input) boundary conditions employed here. 

Clearly these small-x measurements imply that the photon must contain [98] a dom- 
inant hadron-like component at low x, since the simple direct 'box' cross section (based 
on the subprocess 7* ((5^)7 qq) yields F^^^^^ ^ as a: — 0, in contrast to the data 
for X < 0.1 in Figs. 11.5 and 11.6. The QCD RG-improved parton distributions of 
the photon arc thus essential for understanding the data on F2{x,Q'^), with its domi- 
nant contributions deriving from q"'{x,Q'^) — q^{x,Q^). It would be also interesting and 
important to extend present measurements [89, 90, 180, 181] of the gluon distribution of 
the photon, g'^{x,Q^), below the presently measured region 0.1 ^ x < 1 where similarly 
g'^{x < 0.1, Q^) is expected to be also somewhat flatter, see Fig. 11.4, in the small-x 
region than previously anticipated [96] . 
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Figure 11.5: Comparison of our parameter-free predictions from Sec. 11.2 denoted by GRSc 
[6], the previous GRV^ [96] and SaS [97] results for F^{x,Q'^) with the recent OPAL 
(1.9 GeV^) small-,T measurements [98] at two fixed lower scales Q^. Some relevant L3 
data [179] are displayed as well. 
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Figure 11.6: As in Fig. 11.5 but at two fixed scales Q^. The recent OPAL small-x data are 
taken from Ref. [98]. 
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11 A The Parton Content of Virtual Photons 



Next we turn to the somewhat more speculative concept and models of 'resolved' virtual 
photons (P^ 7^ 0). As motivated in Chap. 9 the dependence on the target virtuahty 
P^ will be entirely taken care of by the boundary conditions whereas the partons will 
be treated as if they were on-shell. Of course, this is quite familiar for the quarks and 
gluons, however, the latter rule also refers to the photon-parton V'^'^^'^\x, Q^) in the target 
photon, i.e., the 'direct' subprocesses, e.g., Y {Q^)^'^^^'^^ ~^ QQ have to be calculated with 
P'^ — (cf. rule (ii) in [6]). Thus we can implement the same DIS-y factorization scheme 
as for real photons in Eq. (9.9), and the structure function F2^^ \x, Q'^) of a transverse 
virtual target photon becomes formally very similar to Eq. (11.7): 



q=u,d,s 

+ + -Fif'M) (11.11) 

X ' X ' 

with Cq^'2,{x), Cg^'2,{x) given in (9.5). The 'direct' heavy (charm) quark contribution is given 
by Eq. (9.13) as for real photons and the 'resolved' charm contribution is as in 

Eq. (9.14) with the gluon distribution g'^{z,iip) g'^^^^\z, 

It is physically compelling to expect a smooth behavior of the observable structure 
function in Eq. (11.11) with respect to P^. Especially, in the real photon hmit Eq. (11.7) 

has to be recovered imposing the condition P^^'\x,Q^) P{x,Q^) as P^ ^ for the 
parton distributions. Note that this requirement is automatically satisfied as soon as it 
has been implemented at the input scale Qq since the PDFs are governed by the same 
evolution equations which is a consequence of our unified approach. 



11.4.1 Boundary Conditions 

The above smoothness requirements are fulfilled by the following boundary conditions for 
P(^'\ cf. Eq. (11.5), 

r^^'\x, = P^) = fS\x, P') = v{P')fU^, P') (11-12) 

in LO as well as in NLO. Here P^ = max (P^, //^) as dictated by continuity in P^ [HI] and 
?y(p2) — (l + p2/77j,2^-2 is a dipole suppression factor with = 0.59 GeV^. The scale P^ 
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is dictated not only by the above mentioned continuity requirement, but also by the fact 
that the hadronic component of p^^^\x, Q^) is probed at the scale — [97, 111, 154] 
where the pointlike component vanishes by definition. The boundary condition in Eq. 
(11.12) guarantees, as should be evident, a far better perturbative stability as compared 
to the situation in [111] where the NLO input differed drastically from its LO counterpart 
(cf. Eq. (8) in Rcf. [111]). 

The evolution to > P^ is now analogous to the case of real photons in the previous 
section and the general solution for the resulting parton distributions is similar to the one 
in Eq. (11.4) and Eq. (11.5), 

r^''\x,Q-) = f;pM')+fS'M') (11.13) 



- r/\x.Q')+r,{P')a 



G)rM') + 5f\{Gl-Gl)s-{x,Q') 



with 5f as in Eq. (11.5) and where f^^^ ■* refers again to the solution of the homogeneous 
RG evolution equations, being driven by the hadronic input in (11.12), which is explicitly 
given by Eq. (9.25). Its parametrization is fixed by the available parametrization for 
f^iXjQ"^) in App. D.l. The inhomogeneous 'pointlikc' solution in (11.13) is explicitly 
given by Eq. (9.24) where L = as{Q'^) / as{P'^)- A parametrization of /p/^ \x,Q'^) in LO 
is thus easily obtained from the one for the real photon /pj(x, Q^) in (11.5) in terms of 
In = In [as(/i^)/Q;s(Q^)], where now «s(yU^) has simply to be replaced by as{P^) as 
described in detail in App. D.2. Furthermore, since our NLO predictions for f"'^^'^\x, Q"^) 
turn out to be rather similar to the LO ones, as will be shown below, the simple analytic 
LO parametrizations for f^^^^\x,Q'^) can be used for NLO calculations as well. This is 
certainly sufficiently accurate and reliable in view of additional model ambiguities inherent 
in the parton distributions of virtual photons. 

It should be emphasized that the RG resummed results in (11.13) are relevant when- 
ever <^ (5^, typically [111,162] ~ ^ Q^, so as to suppress power-like (possibly 
higher twist) terms (P^/Q^)" which would spoil the dominance of the resummed loga- 
rithmic contributions. For P^ approaching Q^, the e'^e" — > e'^e'X reaction, for example, 
should be simply described by the full fixed order box 7*((5^)7*(-P^) ~^ QQ keeping all 
(P^/Q^)" terms (cf. Ghap. 8). A calculation of the full perturbative 0{as) corrections 
to this virtual box [108] will offer the possibihty to determine reliably at what values of 
P^ (and possibly x) this 0{as) corrected virtual box becomes the more appropriate and 
correct description. Similar remarks hold for a DIS process ep — > eX, i.e. ^{P'^)p — > X, 
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where 0{as) corrections to pointlike virtual 7(P^)-parton subprocesses have to be ana- 
lyzed in detail in order to decide at what P^ these pointlike expressions become the more 
appropriate description and the virtual photonic parton distributions (i.e., resummations) 
become irrelevant. 

Our unified approach implies that the 'direct' photon contribution to any process 
whatsoever should always be calculated as if this photon is real apart from the fact that 
its fiux should be evaluated according to Eq. (7.37) with 7^ [155,162,163]. This 
differs from the somewhat inconsistent procedure adopted by the HERA-Hl collabora- 
tion [157-161] where exact eq eqg and eg — > eqq matrix elements were used for 
the direct photon contribution to the dijet cross section. As long as ^ ^ Q'^, the 
exact treatment of matrix elements, however, should not differ too much [164] from the 
more appropriate treatment described above. To conclude let us stress that our unified 
approach, as illustrated by the foregoing examples, is not a free option but a necessary 
consistency condition for introducing the concept of the resolved parton content of the 
virtual photon as an alternative to a non-resummed fixed order perturbative analysis at 
P^ 7^ 0. This consistency requirement is related to the fact that all the resolved con- 
tributions due to p^^^\x,Q^) are calculated (evoluted) as if these partons are massless 
[7 — 11] (i.e. employing photon splitting functions for real photons, etc.) in spite of the 
fact that their actual virtuality is given by P^ 7^ 0. Thus the direct photon contribution 
should obviously be also treated accordingly. 

11.4.2 Quantitative Results 

Now we turn to our quantitative predictions and we first present in Fig. 11.7 detailed 
predictions for F2^^ \x, Q"^) for various virtuahties P^ and scales Q'^. Since the 'pointhke' 
component in (11.13) is uniquely calculable perturbatively, a detailed measurement of 
the X and P^ dependence at various fixed values of Q^, as shown in Fig. 11.7, would 
shed light on the theoretically more speculative and far less understood nonperturbative 
'hadronic' contribution in Eq. (11.13) and eventually estabhsh the absolute perturbative 
predictions. Our LO and NLO predictions in Fig. 11.7 show a remarkable perturbative 
stabihty throughout the whole x-region shown, except perhaps for P^ » 1 GeV^ where 
the perturbatively very stable [5,116] 'hadronic' component in (11.13) becomes strongly 
suppressed with respect to the 'pointlike' solution which is less stable in the small x region, 
X < 10~^, as is evident from the right hand side of Fig. 11.7. 
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Figure 11.7: LO and NLO predictions for the x dependence of the virtual photon structure 
function F^''^ ^ at = 10 and 100 GeV^ and various fixed values of ^ Q^, neglecting 
any heavy quark contribution. For comparison we also show the NLO Gluck-Reya-Stratmann 
(GRS) [111] results as well as the predictions for a real (P^ = 0) photon. Notice that the 
dotted curve for = referred to as GRS obviously coincides with the GRV-,, result [96, 112]. 
The results have been multiplied by the number indicated in brackets. 

The individual LO and NLO parton distributions of the virtual photon at — 
10 GeV^ are shown in Fig. 11.8 where they are compared with the ones of Gliick, Reya 
and Stratmann (GRS) [111]. The LO SaS expectations [97] are compared with our LO 
predictions in Fig. 11.9. 

In Fig. 11.10 we show our predictions for xu"''^^'^\x,Q'^) and xg'^^^^\x, 

Q2) with par- 
ticular emphasis on the very small x region. For comparison we also show the results for a 
real (P^ = 0) photon. Plotting the 'hadronic' component in (11.13) separately in the up- 
per two panels of Fig. 11.10 demonstrates that the perturbative 'pointlike' component in 
(11.13) dominates for x > 10^^. In the lower two panels of Fig. 11.10. at P^ = 5, 10 GeV^, 
the hadronic part is nearly vanishing even at small x [?7(P^(GeV^) = 5, 10) = 0.01, 0.003]. 
Furthermore the expected perturbative stability of our present LO and NLO predictions 
is fulfilled. This is in contrast to the GRS results which are unstable [111] throughout 
the whole x-region for P^ >, 1 GeV^, as illustrated in the lower half of Fig. 11.10 at 
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Figure 11.8: Comparison of our predicted LO and NLO(DIS^) distributions of the virtual 
photon at = 10 GeV^ and various fixed values of <^ Q'^ with the GRS densities [111]. 
The curves refer from top to bottom to = 0, 0.2 and 1 GeV^, respectively. The results 
for the real photon (P^ = 0) are very similar to the ones in Fig. 11.2 where the GRV^ curves 
practically coincide with GRS. 

— 100 GeV^, due to the very different perturbative (box) input in LO and NLO [111]. 
In general, however, as soon as the perturbatively very stable 'hadronic' component in 
(11.13) becomes suppressed for P^ » 1 GeV^, the remaining perturbatively less stable 
'pointlike' component destabilizes the total results for q'^^^^\x,Q'^) in the very small x 
region, x < 10-^ as can be seen in Fig. 11.10 for u'^^^^^ at = 100 GeV^ (cf. Fig. 11.7). 

Finally in Fig. 11.11 we confront our LO predictions for p^^^\x, Q"^) with the effective 
parton density 

f-^'\x, Q-) = {^''"^' + ^'''"'') + i 9''"^' (11-14) 

q=u,d,s 
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Figure 11.9: Our LO distributions as in Fig. 11.8 compared with the ones of SaS [97]. 

extracted in LO from DIS dijet data by the HERA-Hl collaboration [157-161] very 
recently. The predicted dependence on the photon's virtuality at the scale = 
{p't^ = 85 GeV^ agrees reasonably well with the measurements in the relevant kine- 
matic region <^ Q^. This is also the case at other scales Q"^ = ^Pr*) ^^'^ fixed values 
of X [157-161] not shown in Fig. 11.11. As discussed above, it should be kept in mind, 
however, that for larger values of approaching Q^, which refer to the dashed curves 
in Fig. 11.11, the whole concept of RG resummed parton distributions of virtual photons 
is not appropriate anymore. Since the resolved contributions of a virtual photon with 
virtuality as large as P^ = 10 — 15 GeV^ are by a factor of about 10 smaller than the ones 
of a real (P^ = 0) photon, it is reasonable to conclude from Fig. 11.11 that for P^ 10 
GeV^ the DIS ep — > eX process considered is dominated by the usual direct 7* = 7(P^) 
exchange cross sections and not 'contaminated' anymore by resolved contributions. 
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Figure 11.10: LO and NLO predictions for the up-quark and gluon distributions of a virtual 
photon 7(P^) at Q"^ = 10 and 100 GeV^. For comparison the results for the real photon 
(P^ = 0) are shown as well. In the upper half the NLO 'hadronic' contribution in (11.13) is 
also shown separately. The GRS expectations are taken from Ref. [111]. The DIS-y results for 
li'''^^^) can be easily converted to the MS scheme with the help of Eq. (9.9), whereas (y'^*^^^-' 
remains the same in both schemes. The results have been multiplied by the numbers indicated 
in brackets. 



11-4 The Parton Content of Virtual Photons 



130 




10 20 30 40 50 60 70 80 



P^/GeV^ 



10 



10 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

— 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

- 85 GeV^ , x = 0.6 


, , , , 1 , , , , 1 , , , , 1 , , , , 





10 20 30 40 50 60 70 80 

P^/GeV^ 



Figure 11.11: Predictions for the LO effective parton density xf'^'^^^\x,Q'^), defined in 
Eq. (11.14), at the scale = (^Pt*) ~ GeV^ and at two fixed values of x. The 
HI data [157-161] have been extracted from DIS ep dijet production. The solid curves re- 
fer to our predictions in the theoretically legitimate region P'^ ^ Q"^ = [t^t^ , whereas the 

dashed curves extend into the kinematic region of larger approaching Q'^ where the concept 
of parton distributions of virtual photons is not valid anymore (see text). 
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11.5 Summary and Conclusions 

The main purpose of the present chapter was to formulate a consistent set of boundary 
conditions which allow for a perturbatively stable LO and NLO calculation of the photonic 
parton distributions p^^'^\x, Q^) as well as for a smooth transition to the parton densities 
of a real {P^ — 0) photon. Employing the pionic distributions f'"{x,Q'^) of Chapter 10, 
required for describing, via VMD, the nonpointlike hadronic components of a photon, we 
arrive at essentially parameter-/ree predictions for f'^^^'^\x, Q^) which are furthermore in 
good agreement with all present measurements of the structure function Q"^) of real 

photons 7 = 7(P^ = 0). It should be noted that the experimentally almost unconstrained 
pionic gluon and sea distributions, g'^{x,Q'^) and q'^{x,Q'^), have been uniquely derived 
from the experimentally rather well known pionic valence density v'^{x,Q'^) and the (re- 
cently updated [26] dynamical) parton distributions of the proton. We have furthermore 
implemented these hadronic components by using a VMD ansatz for a coherent super- 
position of vector mesons which maximally enhances the contributions of the up-quarks 
to F2 as favored by all present data. Since these hadronic contributions are generated 
from the valence-like input parton distributions at the universal target-mass independent 
low resolution scale = /U^ — 0.3 GcV^, we arrive, at least for real (P^ = 0) photons, 
at unique small-x predictions for x <, 10^^ at > fi"^ which arc of purely dynamical 
origin, as in the case of hadrons. As has been demonstrated in Sec. 11.3 these small-x 
predictions are in perfect agreement with recent small-x measurements of the photon 
structure function Q^) [98] at all presently accessible values of Q^. Furthermore, 

since our universal input scale /J,^ fixes also uniquely the perturbative pointlike part of 
the photonic parton distributions, which dominates for a: ^ 0.1, the large-x behavior of 
photonic structure functions is unambiguously predicted as well. 

Our expectations for the parton content of virtual (P^ 7^ 0) photons are clearly more 
speculative, depending on how one models the hadronic component (input) of a vir- 
tual photon. The latter is usually assumed to be similar to the VMD input for a real 
photon, times a dipole suppression factor which derives from an effective vector-meson 
P^-propagator, cf. Eq. (11.12). Whenever a virtual photon is probed at a scale ^ P^, 
with its virtuality being entirely taken care of by the (transverse) equivalent photon flux 
factor and the boundary conditions for the photonic parton distributions, all cross sec- 
tions of partonic subprocesses (resolved and direct) should be calculated as if P^ = 0. In 
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other words, the treatment and expressions for p^^'^\x, Q^) as on-shell transverse partons 
obeying the usual RG Q^-evolution equations (with the usual splitting functions of real 
photons, etc.) dictate an identification of the relevant sub-cross-sections p^^^^X — > X' 
with that of the real photon, a{P'^^^^X X') — a{px — > X'). In particular, the 
calculation of F2^^ \x,Q'^) requires the same photonic Wilson coefficient C^{x) as for 

= 0. This allows to formulate similar boundary conditions in LO and NLO which 
give rise to perturbatively stable parton distributions, cross sections (i.e. also structure 
functions) of virtual photons 7(-P^) as long as they arc probed at scales ^ where 

= Ami, {p^T^ ' ^^^-i si'iid typically ^ TqQ'^- should be emphasized that only 
in this latter kinematic range P^ ^ is the whole concept of RG resummed parton 
distributions of (resolved) virtual photons appropriate and relevant. Parton distributions 
of virtual photons extracted recently from DIS ep dijet data are in good agreement with 
our (parameter-free) predictions. 

Finally, we present simple analytic parametrizations of our predicted LO and 
NLO(DIS^) parton distributions of real photons. Prom these LO parametrizations one 
can easily obtain also the ones for a virtual photon which, whithin sufficient accuracy, may 
also be used in NLO. Our NLO(DIS^) parametrizations of the parton densities of the real 
photon can be easily transformed to the MS scheme according to Eq. (9.9) which might 
be relevant for future NLO analyses of resolved photon contributions to hard processes 
where most NLO subprocesses have so far been calculated in the MS scheme. 

A FORTRAN package containing our most recent parametrizations can be obtained 
by electronic mail on request. 
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Has the QCD RG— Improved Parton 
Content of Virtual Photons been 
Observed? 

It is demonstrated that present e^e~ and DIS ep data on the structure of the virtual 
photon can be understood entirely in terms of the standard 'naive' quark-parton model 
box approach. Thus the QCD renormalization group (RG) improved parton distributions 
of virtual photons, in particular their gluonic component, have not yet been observed. 
The appropriate kinematical regions for their future observation are pointed out as well 
as suitable measurements which may demonstrate their relevance. The results presented 
in this chapter are taken from [8]. 

12.1 Introduction 

Recent measurements and experimental studies of dijet events in deep inelastic ep [161] 
and of double-tagged e+e~ [105] reactions have indicated a necessity for assigning a (QCD 
resummed) parton content of virtual photons 7(-P^) as suggested and predicted theoreti- 
cally [6, 97, 104, 111, 125-127, 153-155]. In particular the DIS dijet production data [161] 
appear to imply a sizeable gluon component g^'^^^\x,Q'^) in the derived effective parton 
density of the virtual photon, where refers to the hadronic scale of the process, Q ~ p^^ , 
or to the virtuality of the probe photon 7*((5^) which probes the virtual target photon 
7(P^) in e+e^ e+e^X. It is the main purpose of this chapter to demonstrate that this 
is not the case and that all present data on virtual photons can be explained entirely in 
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terms of the conventional QED doubly virtual box contribution 7*(Q^)7*(-P^) ~^ QQ in 
fixed order perturbation theory -sometimes also referred to as the quark-parton model 
(QPM). 

This is of course in contrast to the well known case of a real photon 7 = 7(P^ = 
0) whose (anti-)quark and gluon content has been already experimentally established 
(for recent reviews sec [89,90]) which result mainly from resummations (inhomogcneous 
evolutions) of the pointlike mass singularities proportional to In Q'^/rrig occurring in the 
box diagram of 7*(Q^)7 qq for the light q = u, d, s quarks. This is in contrast to 
a virtual photon target where 7*(Q^)7*(P^) qq does not give rise to coUinear (mass) 
singularities but instead just to finite contributions proportional to In /P^ which a 
priori need not be resummed to all orders in QCD. 

This chapter is organized as follows: 
In Sec. 12.2 we compare fixed order box and QCD resummed expectations for the effective 
structure function Feff with present e'^e" data, while Sec. 12.3 contains a comparison of an 
effective leading order parton density for the virtual photon with DIS ep data. Suggestions 
of experimental signatures which can probe the QCD parton content, in particular the 
gluon content of virtual photons are presented in Sec. 12.4 and our conclusions are finally 
drawn in Section 12.5. 

12.2 RG— Improved Parton Model Expectations for 
the Effective Structure Function F^^ 

As promised in Chap. 8 we shall now turn to a quantitative study of the various QED- 
box and QCD Q^-evoluted structure function expectations for a virtual photon target 
and confront them with all presently available e+e" data of PLUTO [106] and the recent 
one of LEP-L3 [105]. The 'asymptotic' and 'full' box results for Fes have already been 
presented in Chap. 8 and details of their calculation can be found there. The QCD 
resummed expectations for F^q according to our model in the previous chapter are given 

by 

F,^{x,Q\P^)^F^^^''"\x,Q^) (12.1) 
since we neglect any effects due to longitudinal target photons in our approach [6]. 
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Figure 12.1: Predictions for Fgg as defined in (8.14). The box results are as in Fig. 8.3. 
The QCD resummed NLO expectations of GRSc [6] for F2 in (12.1) turn out to be similar 
to the LO ones [6]. Also shown are the LO-resummed results of SaS ID [97] for F2 and 
Feff = F2 + |-Plt (see text). The total charm contribution to the latter two QCD results 

involves also a 'resolved' component -F^^^ ''''^^^{x.Q'^) according to Eq. (9.14) which turns out 
to be small as compared to the box contribution shown which dominates in the kinematic 
region considered. The PLUTO data are taken from [106]. 
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Figure 12.2: As in Fig. 12.1, but for ^ 120 GeV^ and = 3.7 GeV^ appropriate for the 
LEP-L3 data [105]. 



A different approach has been suggested by Schuler and Sjostrand [97]. Apart from 

using somewhat different input scales Qq and parton densities, the perturbatively exactly 

calculable box expressions for A^ ^ ^ (cf. App. C.2) are, together with their 

LO-QCD Q^-evolutions, extrapolated to the case of real photons = by employing 

some dispersion-integral-like relations. These link perturbative and non-perturbative 

contributions and allow a smooth limit P^ 0. (Note, however, that the LO Q^- 

evolutions are performed by using again the splitting functions of real photons and on- 

shell partons.) Since one works here explicitly with virtual (P^ 7^ 0) expressions, the 

longitudinal contributions of the virtual photon target should be also taken into account 

when calculating Pefr = P2 + fa^pLT with Fab = Q^/(47r^a)(a;/3)~Vafe (a = (L,T), h ~ 

(L, T)) and where is given in Eq. (8.11), as described for example in [111], which is in 

contrast to our approach in (12.1).^ 

"'^In an alternative approach [130, 132] one may consider the longitudinal component of the virtual 
photon target 7l(-P^) to possess, like the transverse component, a universal process independent hadronic 
content obtained radiatively via the standard homogeneous (Altarelli-Parisi) Q^-evolution equations with 
the boundary conditions for the pointlikc component at = given in Eq. (9.55). We have checked 
that the predictions for Fgg{x, Q^,P^) obtained in this approach differ only slightly (typically about 10% 
or less) from those of the standard fixed order perturbative approach at presently relevant kinematical 
regions (P^ < j^Q"^, x > 0.05) due to the smallness of F^.^ relative to F^^ deriving from a^^^ and a^rj,, 
respectively, in (8.17). 
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Despite the limited statistics of present data the box predictions for F^s in (8.14) shown 
in Figs. 12.1 and 12.2 appear to be in even better agreement with present measurements 
than the QCD resummed expectations of SaS [97] and GRSc [6]. Typical QCD effects 
like the increase in the small-x region in Fig. 12.1, being partly caused by the presence of 
a finite gluon content g'^^^'^\x, Q^), cannot be dehneated with the present poor statistics 
data. 

These results clearly demonstrate that the naive QPM predictions derived from the 
doubly virtual box 7*((5^)7*(-P^) qq fully reproduce all e+e~ data on the structure of 
virtual photons 7(P^). In other words, there is no sign of a QCD resummed parton content 
in virtual photons in present data, in particular of a finite gluon content 
which is absent in the 'naive' box (QPM) approach. 

Characteristic possible signatures for QCD effects which are caused by the presence 
of a finite and dominant gluon component g''^^^'' will be discussed in Sec. 12.4. 



12.3 Comparison of Theoretical Expectations with 
DIS ep Data and Effective Quark Distributions 
of Virtual Photons 

In oder to extract the parton densities of virtual photons from DIS ep dijet data, the HI 
collaboration [161] has adopted the 'single effective subprocess approximation' [183] which 
exploits the fact that the dominant contributions to the cross section in LO-QCD comes 
from the 2 — > 2 parton-parton hard scattering subprocesses that have similar shapes and 
thus differ mainly by their associated color factors. Therefore the sum over the partonic 
subprocesses can be replaced by a single effective subprocess cross section and effective 
parton densities for the virtual photon given by 

/>(^^)(x,g^) = [^'^''^\^^Q") ^Q'^'^Kx.Q')] +\9'^''\x.Q') (12.2) 

q=u, d, s 

with a similar relation for the proton /^(x, Q'^) which is assumed to be known. It should 
be emphasized that such an effective procedure does not hold in NLO where all additional 
(very different) 2 — > 3 subprocesses contribute [164] . This NLO analysis affords therefore 
a confrontation with more detailed data on the triple-differential dijet cross section as 
compared to presently available data [161] which are not yet sufficient for examining the 
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Figure 12.3: Predictions for the effective parton density defined in Eq. (12.2). The 'box' results 
refer to the universal g^^^ ■* in (12.3), and the 'effective' ones to q^lf as defined in (12.4) 
as derived from the full box expressions (8.11) including all 0{P'^ /Q"^) contributions. The 
LO-QCD predictions of GRSc [6] are shown by the solid curves which refer to the predictions 
in the theoretically legitinnate region <^ Q^, whereas the dashed curves extend into the 
kinematic region of larger approaching where the concept of QCD-resummed parton 
distributions of virtual photons is not valid anymore. (Note that the results for x = 0.6 
terminate at ^ 54 GeV^ due to the kinematic constraint W"^ > 0, with being defined 
in (7.1), i.e. x < (1 + P^/Q^)"^-) For illustration we also show the effective LO-QCD parton 
density /'"' of a real photon 7 = 7(P^ = 0) of GRSc [6] multiplied by the simple p-pole 
suppression factor rj{P'^) in (11.12) which clearly underestimates the Hl-data [161]. 
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relative contributions of q^^^^^XjQ^) and g'^^^'^\x,Q^). In Fig. 12.3 we compare our 
LO RG-resummed predictions for f^^^^^XjQ"^) with the naive non-resummed universal 
(process independent) box expressions [cf. Sec. 9.3.4] 

Q') - Q') - + (1 - ^f] In % . (12.3) 

Although the fully QCD-resummed results are sizeable and somewhat larger in the small 
region than the universal box expectations, present HI data [161] at = (Px*)^ ~ ^5 
GeV^ cannot definitely distinguish between these predictions. It should be furthermore 
noted that the QCD gluon contribution , Q'^) is suppressed at the large values of 

X shown in Fig. 12.3. Therefore present data [161] cannot discriminate between the finite 
QCD resummed component g'^^^'^\x,Q'^) and the non-resummed g^l^ \x,Q'^) = 0. 

It is obvious that these two results shown in Fig. 12.3 are only appropriate for vir- 
tualities < Q^^ typically ^ iq to 20 GeV^ at ^ 35 GeV^, since 0{PyQ^) 
contributions are neglected in RG-resummations as well as in the definition (12.3). In or- 
der to demonstrate the importance of 0{P'^/Q'^) power corrections in the large region 
let us define, generalizing the definition (9.34), some effective (anti)quark distributions as 
common via 

lp,V(a;,g^)^ Y: elU:P{x,Q-)+q:P{x,Q')] (12.4) 

X L J 

q=u, d, s 

where, of course, qjl^^ = q^jf and the full box expression for Plboxl^''^^) can be ob- 
tained from the general definition in Eq. (7.25) and the full box results in Eq. (8.11) 
utilizing m = ruq = 0, i.e. A = and summing over the light quarks. The full box expres- 
sions imply again gjjf \x, Q"^) = in contrast to the QCD resummed gluon distribution. 
The qllf "* introduced in (12.4) is, in contrast to (12.3), of course non-universal. The 
'effective' results shown in Fig. 12.3 clearly demonstrate the importance of the 0{P^/Q'^) 
terms at larger values of P^ ^ which are not taken into account by the QCD re- 
summations and by the universal box expressions in (12.3) also shown in Fig. 12.3. It is 
interesting that the non-universal qj^^ ^ defined via F2 in (12.4) describes the Hl-data 
at large values of P^ in Fig. 12.3 remarkably well. This may be accidental and it remains 
to be seen whether future LO analyses will indicate the general relevance of q^g' ^ {x, Q"^) 
in the large P^ region. 

As we have seen, present DIS dijet data cannot discriminate between the universal 
naive box and QCD-resummed expectations in the theoretically relevant region P'^ <ti Q^, 
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Figure 12.4: Expected charm contributions to F2. The naive 'direct (box)' result refers to 
(9.13) [or (C.12)] and the LO-QCD 'resolved' prediction is given in (9.14) with g^^^^\x, Ami) 
taken from GRSc [6]. This latter 'resolved' contribution is absent in the naive box (QPM) 
approach. 

mainly because these data are insensitive to the gluon content in 7(i-*^) generated by 
QCD-evolutions which is absent within the naive box approach. Therefore we finally 
turn to a brief discussion where such typical QCD effects may be observed and delineated 
by future experiments. 



12.4 Possible Signatures for the QCD Parton Con- 
tent of Virtual Photons 

Since e+e" and DIS ep dijet data cannot, at present, delineate the QCD-resummed parton 
content of a virtual photon, in particular not its gluon content, we shall now propose and 
discuss a few cases where such typical QCD effects may be observed and possibly confirmed 
by future experiments. 

Charm production in e'^e" e'^e" ccX would be a classical possibility to delineate 
such effects due to a nonvanishing g'^^^^\x, Q^). In Fig. 12.4 we compare the usual (fixed 
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Figure 12.5: Predictions for the total light quark E'''^^^^ = 2Sq=u_d,s and gluon contri- 

butions to the effective parton density in (12.2) at a fixed scale — 85 GeV^ and two fixed 
virtualities P'^ — 1 and 5 GeV^. The naive box results refer to the universal q^^^ ^ defined in 

(12.3) and to qj^^'^ defined in (12.4). The LO-QCD RG -resummed predictions are denoted 
by T."'^^ ) and g"'^^ ^ according to GRSc [6]. The latter gluon contribution is absent in the 
naive box (QPM) approach. 
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order) 'direct' box contribution to F| with the 'resolved' gluon-initiated one as given by 
(9.14). The 'direct' box contribution entirely dominates in the large-x region, x ^ 0.05, 
accessible by present experiments (cf. Figs. 12.1 and 12.2), whereas the typical QCD- 
resummed 'resolved' contribution becomes comparable to the 'direct' one and eventually 
dominates in the small-x region, x < 0.05. Thus a careful measurement of the charm 
contribution to F2 at x ^ 0.05 would shed some light on the QCD parton (gluon) content 
of virtual photons, since such a 'resolved' contribution in Fig. 12.4 would be absent within 
the naive box approach. 

The effective parton distribution f'^^^^\x,Q'^) in (12.2) at not too large values of x 
and P^, as may be extracted in LO from DIS ep dijet data, would be another possibility to 
observe QCD-resummation effects due to a nonvanishing gluon component g^'^^^^x , Q"^) . 
In Fig. 12.5 we show the quark and gluon contributions to f^^^^^i in (12.2) separately. The 
box (anti-)quark contributions, which are similar to the QCD-resummed ones, entirely 
dominate over the QCD-resummed gluon contribution in the large-x region, x ^ 0.4, 
accessible to present experiments (cf. Fig. 12.3). Only below x ~ 0.3 does the QCD gluon 
contribution become comparable to the (anti-) quark components and dominates, as usual, 
for X 0.1. It should be remembered that g^^^ \x, Q'^) — 0. Furthermore, the increase of 
the RG-resummed q^^^^\x^Q'^) at small x in Fig. 12.5 is induced by the vector-meson- 
dominance-like input for the Q^-evolution of the 'hadronic' component of the photon's 
parton distribution [6, 97] and is disregarded in our naive box (QPM) analysis. 

Thus a measurement of dijets produced in DIS ep reactions in the not too large x 
region, x ^ 0.3, would probe the QCD parton content of virtual photons, in particular 
their gluon content which is absent in the naive QPM box approach. In this region, 
and at not too large photon virtuahties 5 GeV^ shown in Fig. 12.5, the 'resolved' 

gluon-dominated contribution of the virtual photon to high jet production at scales 
Q = Et ~ 5 — 10 GeV exceeds by far the 'direct' box-like contribution of a pointlike 
virtual photon [162]. The production of prompt photons at HERA via a tagged DIS 
process ep — > ejX may offer additional probes of the gluonic content of virtual photons 
[184]. 
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12.5 Summary and Conclusions 

Virtual photons 7(-P^), probed at a large scale » P^, may be described either by 
fixed order perturbation theory, which in lowest order of QCD yield the quarks and 
anti-quarks generated by the universal part of the 'box' diagram, or alternatively by 
their renormalization group (RG) improved counterparts including particularly the gluon 
distribution g^^P''>{x,Q^). 

The results in Sections 12.2 and 12.3 demonstrate that all presently available e~^e~ 
and DIS ep dijct data can be fully accounted for by the standard doubly virtual QED 
box diagram and are not yet sensitive to RG resummation effects which are manifest only 
in the presently unexplored low-x region of the parton distributions in 7(P^). In fact, 
as shown in Section 12.4, these resummation effects start to dominate only at x < 0.3 
and may be observed by future measurements at — 0{1 GeV^) of a{ep —>■ ejjX) 
or <j(e~^e~ e+e~ ccX) at high energy collisions. These measurements could finally 
discriminate between the fixed order and RG-improved parton distributions of the virtual 
photon. 



Chapter 13 



Summary 



It was the main theme of Part II of this thesis to analyze the parton content of pions and 
real and virtual photons in leading order (LO) and next-to-leading order (NLO) of QCD 
within the framework of the recently updated [26] radiative parton model [26,41,57,92, 
93]. 

We started with a detailed consideration of the ee — > eeX scattering process being 
presently the main source of information on the structure of real and virtual photons. We 
defined structure functions for virtual photons in a model independent way and related 
them to the photon-photon scattering cross sections. These structure functions can be 
measured in deep inelastic electron-photon scattering (DISe7) ^® demonstrated the 
factorization of the cross section for the process ee — > eeX into a fiux of photons (produced 
by bremsstrahlung off the incident lepton) times the DlSe-y cross section for the general 

7^ case in the Bjorken hmit <^ Q^. However, we found that the neglected terms 
are of the order 0{-\/xP'^ /Q"^) only and not of the order 0{xP'^ /Q'^) as might have been 
naively expected. Clearly, a numerical comparison of the factorization formula with the 
exact ee — > eeX cross section would be interesting in order to clarify when xP'^/Q'^ is 
small enough such that the factorization formula is a good approximation. 

From these general kinematical considerations we turned to a calculation of the 
photon-photon cross sections a ah according to the doubly virtual box 7* ((5^)7* (P^) qq 
in lowest order perturbation theory. We utilized these results in order to compute the ef- 
fective structure function Pefr(a;, ■, P^) in fixed order perturbation theory (FOPT) which 
we compared with all present e+e^ virtual photon data [105, 106]. It came out that the 
LO-box predictions for Peff are in agreement with the present low statistics data even in 
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the case of the PLUTO data [106] with = 0.35 GeV^ which is just in the transition 
region from deeply virtual (P^ » A^) to real (P^ ~ 0) photons where non-perturbative 
effects are expected to become increasingly important, especially in the small-x region. 
More precise data for Peff in this 'transition region' with intermediate target virtualities 
would be highly desirable in order to observe the onset of such non-perturbative effects. 
On the other hand, NLO corrections to the doubly virtual box [108] offer the possibility to 
investigate the perturbative stability of the fixed order results and to see if the description 
of the data further improves in NLO accuracy. 

After presenting the complete theoretical framework we started our phenomenological 
studies of radiatively generated parton distributions with an analysis of the parton content 
of the pion in LO and NLO QCD. Since only the pionic valence density v'^{x, Q'^) is 
experimentally rather well known at present, we utilized a constituent quark model [94, 95] 
in order to unambiguously relate the pionic light sea q'^{x, Q"^) and gluon g'^ix, to the 
much better known (recently updated) parton distributions of the proton [26] and the 
pionic valence density v'^{x, Q"^). 

Next we formulated a consistent set of boundary conditions which allowed for a pertur- 
batively stable LO and NLO calculation of the photonic parton distributions p^^^\x, Q^) 
as well as for a smooth transition to the parton densities of a real (P^ = 0) photon. Em- 
ploying the above summarized pionic distributions f^iXjQ^), required for describing -via 
vector meson dominance (VMD)- the nonpointlike hadronic components of a photon, we 
arrived at essentially parameter-/ree predictions for p^^'^\x^ Q"^) which turned out to be 
in good agreement with all present measurements of the structure function F^(a;, Q^) of 
real photons 7 = 7(P^ = 0). More specifically, we implemented these hadronic com- 
ponents by using a VMD ansatz for a coherent superposition of vector mesons which 
maximally enhances the contributions of the up-quarks to P^(a;, Q^) as favored by all 
present data. Since these hadronic contributions arc generated from the valence-like in- 
put parton distributions at the universal target-mass independent low resolution scale 
Qq = /i^ ~ 0.3 GeV^, we arrived, at least for real (P^ = 0) photons, at unique small-x 
predictions for x <, 10"^ at > fi'^ which are of purely dynamical origin, as in the case of 
hadrons. A comparison with most recent small-x measurements of the photon structure 
function P^(a:, Q'^) [98] showed that these small-,x predictions arc in perfect agreement 
with the data at all presently accessible vahics of Q^. The agreement with the data on 
P^(a;, Q^) confirms our radiatively generated quark distributions q'^{x,Q^) = (p{x,Q^) 
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while the gluon is probed only indirectly at small-x due to the evolution. It would be also 
interesting and important to extend present measurements [89, 90, 180, 181] of the gluon 
distribution of the photon, g^i^XjQ"^), below the presently measured region 0.1 ^ x < 1 
where g'^{x < 0.1, Q^) is expected to be somewhat flatter in the small-x region than 
previously anticipated [96]. Our expectations for the parton content of virtual (P^ ^ 0) 
photons are clearly more speculative, depending on how one models the hadronic com- 
ponent (input) of a virtual photon. We assumed the latter to be similar to the VMD 
input for a real photon times a dipole suppression factor which derives from an effective 
vector-meson P^-propagator representing somehow the simpliest choice. We formulated 
similar boundary conditions in LO and NLO which gave rise to perturbatively stable 
parton distributions and cross sections (i.e. also structure functions of virtual photons) 
as long as they are probed at scales ^ P^ where is a hard scale of the process 
under consideration. We showed that effective parton distributions of virtual photons ex- 
tracted recently from DIS ep dijet data are in good agreement with our (parameter-free) 
predictions. 

Finally, we turned to more detailed tests of our model for the parton content of virtual 
photons and we compared these QCD-resummed expectations with the above summarized 
fixed order box expressions. Our results demonstrated that all presently available e"^e~ 
and DIS ep dijet data can be fully accounted for by the standard doubly virtual QED 
box diagram and are not yet sensitive to RG resummation effects which are manifest 
only in the presently unexplored \ow-x region of the parton distributions in 7(P^). It 
turned out that these resummation effects start to dominate only at x < 0.3 (the larger 
P^ the smaller x) and may be observed by future measurements at P^ = 0{1 GeV^) of 
a{ep — > ejjX) or a{e'^e'' — > e'^e~ ccX) at high energy collisions. These measurements 
could finally discriminate between the fixed order and RG-improved parton distributions 
of the virtual photon. 

FORTRAN packages containing simple analytic parametrizations of our most recent 
parton distributions of the pion and the (real and virtual) photon are available upon 
request by electronic mail from schien@zylon.physik.uni-dortmund.de. 
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Deep Inelastic Scattering on Massive 
Quarks at 0{al) 



A.l Real Gluon Emission 



We define a partonic tensor 

= J2 J2^Q2{P2),9{k) iQ^riv - A75)gi| gi(pi)) x (/. ^ uy (a.i) 

color spin 

which can be decomposed into its difi^erent tensor components as usual 



? g^'' + a;? p'ip\ + iu^ £ + o'? q^q" + u'^ (gX + ? V) ■ (A.2) 

uj^y can be easily calculated from the general Feynman rules for invariant matrix elements 
which are customarily expressed as functions of the Mandelstam variables s = {pi + qY 
and i = {pi — k)'^ to which we will refer in the following. Projection onto the individual 
'^1^1,2,3 ^1- (A- 2) is performed for nonzero masses and in n = 4 + 2£ dimensions with 
the following operators 

X 4A-^[ml,s,-Q^] } 
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such that Pi ■ ou — cuf. The normahzation in Eqs. (A.l), (A. 2) is such that real gluon 
emission contributes F/^ to the hadronic structure functions via 
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where [185] 



dPS 



1 s — ml 
Stt i 



r(i + £) 



(s — m 



2\2 
2) 



Atts 



[y{l - y)Y dy (A.5) 



is the partonic phase space. In Eq. (A.5) y is related to the partonic center of mass 
scattering angle 9* and the partonic Mandelstam variable i via 

1 



y = - (1 + cosr) 
1 



2AK,s,-g2 



2s{i - mj) 



s — rrin 



(A.6) 



We have chosen dimensional regularization for the soft gluon poles stemming from s — > m| 
which arise from propagators in the uji times phase space factors in dPS. In Eq. (A.4) we 
use [48] 



(s - m^) 



2\2£-l 



1-^ 



2e-l 



+ 0(e) (A.7) 



which separates hard gluon emission (~ ff) from soft gluon {Si) contributions in Eq. 
(3.10). Note that in Eqs. (3.9), (3.10) the integration variable ^, which is implicitly defined 
in Eq. (3.2), has been changed to ^' = x/^ for an easier handhng of the distributions. For 
the relation between s and ^' see Eq. (A. 39) below. 

Since all quark masses are kept nonzero, no poles in y (coUinear singularities) are 
contained in the integration volume. The which occur in Eq. (3.10) and which are 
given below in Eq. (A. 40) are therefore straightforward integrals of the 
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and the in Eqs. (3.10), (A. 37) pick up the pole in Eq. (A. 7) 



dy [y(i-y)]^ <^i(s-m2 



2\2 



(A.9) 



where the proportionahty is given by kinematical and phase space factors which must be 
kept up to 0{e). 

The normahzation of our hadronic structure functions in Eq. (A. 4) can be clearly 
inferred from the corresponding LO results in Eq. (3.3). Nevertheless, for definiteness we 
also give the hadronic differential cross section to which it corresponds 
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1 (Gf '^V (Cf ^? 2M^Ei 



dxdy ni 

X 



27r 



Si,+ (1 - y)F2 + Si,+ y'xF, ± Ri,+ 2y{l - ^)xFs 



(A.IO) 



where {Cfff = [^f//(Q^ + M^) g^g^/{Q^ + M|,)]V2 is the effective squared gauge 
coupling -including the gauge boson propagator- of the •y^iV — A75) lepton and quark 
current, respectively, and rii counts the spin degrees of freedom of the lepton, e.g. = 1, 2 
for I = v^e^ . The leptonic couplings 5*;,+, Ri^+ are defined analogous to the quark couplings 
in Eq. (3.6). As noted below Eq. (3.6), B = B' for non-interference (pure B scattering). 

A. 2 Vertex Correction 
A. 2.1 Results 



We have calculated the vertex correction in n = 4 + 2e dimensions at 0{al) for gen- 
eral masses and couplings using the Feynman gauge. The unrenormalized vertex Aq 
[Fig. 3.1(c.l)] has the structure 



a. 



A^ = CF-^r{l-e 



Q 



2 \ £ 



Co,- YL, + C+ YR^ 



+ Ci,_ m2 L5 + Ci,+ mi i?5 + Cg,_ mag'' L5 + Cg,+ m^q^" R5 > (A.ll) 



with L^ — {V — A ^^), R5— (V + A 75). The coefficients read 



Co,- = -(-l-E++/i) + 



A2 / , /g2 
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+ -In 



+ - In 



+ 



m% — m 



X 



In^ 



A-E, 



2Q2 



+ -ln2 
2 



2Q2 



iln 



mi 



- Lio 



2A 



— Lio 



2g2 

A-E_+ 
2A~~ 



In^ 



+ 



A + E, 



(A.12) 



+ Li2 



2g2 

A + E+. 
2A 



In^ 



A + E_ 



+ Li2 



2Q^ 

A + S.. 
2A 



Cn.4- = -TTT 



^9: 

with 



C+ = 2mim2/i 
-1 

-1 
1 
1 



E,_A + In ( 



E_ , /i - In 



m 



mi 



2\ 1 



(A^ - 2mlQ^) h - + E^_ In f 4 
(-A^ + 2m^Q2 - E_+g^) 7i + 2Q2 + (E_+ + Q^) In 



A = -ln 



- A, 



72 = /ilnA- — 
A 



(A.13) 

(A.14) 
(A.15) 



X < - In 
2 



Li2 



A - E, 



2Q2 

A - S+_ 
2A~~ 



In^ 



A-E_ 



2g2 



In^ 



A + E, 



2g2 



In^ 



A + E, 



2g2 



2A 



2A 



+ Li2 



A + E_ 
2A~ 



The renormahzed vertex [Fig. 3.1 (c.l)-(c.3)] is obtained by wave function renormal- 
ization: 



A^ = 7'^L5(Zi-l) + A(^ + 0(a^) 



(A.16) 

where Zi = a/^2(pi)^2(P2)- The fermion wave function renormahzation constants are 
defined on mass sheU: 

,2 



-[3-4£ + e»(£^)] 



such that 



a. 



2, = l + C.j|r(l-e, 



£ 2 \m{ 



(A.17) 



(A.18) 
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The final result for the renormalized vertex reads 



+ Ci,_ m2 Vi U + Ci,+ mi Vi Rb + - ^2 U + Cg,+ rm q>' R5 }{A.19) 



with C+, Ci^±, Cg^± as given above and 




4+ ^ 



(A.20) 









2Q2 


2 


2g2 



A + S, 
2A 



+ Li2 



A + S_ 
2A" 




A. 2. 2 Calculation 



Basic ingredient in the calculation of one-loop virtual corrections are one-loop integrals 
classified according to the number N of propagator factors in the denominator and the 
number P of integration momenta in the numerator. Integrals with N = 1,2,3,4 are 
usually called one-, two-, three- and four-point functions (IPF . . . 4PF). For P+4— 2A^ > 
these integrals are UV-divergent. In dimensional regularization these divergences are 
regulated by evaluating the integrals in general dimensions n 7^ 4 and the divergences 
become manifest as poles in the limit n — > 4. 

We define the following one-loop tensor integrals^ (see for example [186, 187]): 



Tai-apiPl^ ■ ■ ■ ,PN-i,mo,... ,mN-i) = 



(2^//) 



l-rt 



d"/ 



IT 



(A.21) 



ZTT- J DqDi ■ ■ ■ Dn-1 

where all the external momenta pi are defined to be incoming and the denominators 
stemming from the propagators in the Feynman diagram are given by 



A = {I + Pif - ml + ie , 



with 



Pi = Y,Pj^ (i = 0,...,A^-l). (A.22) 



^Noto that this definition is related to Eq. (4.1) in [186] as foUows: 

mo, . . . ,mjv-i). 



pi - MJ-f-Pi' • • • ,PN-i,mo, . . . ,mjv-ij = T^^'...^p (pi, . . . ,pn-i, 
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To achieve a cyclic symmetry we furthermore identify the 0-th with the A^-th propagator: 
Disr = Do, thn = trio and pN — po — 0. Furthermore it is convenient to define 

Pio = Pi , Pij ^Pi- Pj ■ (A.23) 

The +ie part in the denominators with an infinitesimal £ > is needed to regulate 
singularities of the integrand and its specific choice ensures causality. After integration it 
determines the correct imaginary parts of the logarithms and dilogarithms. The arbitrary 
mass scale /i has been introduced such that the integrals have an integer mass dimension. 
Conventionally is denoted by the A'"-th character of the alphabet, i.e. = A, 

= B, . . . , and the scalar integrals carry an index 0. 

The tensor integrals in (A. 21) can be related to the scalar integrals Aq, Bq, Cq and 
Dq by a Passarino-Veltman decomposition [188] which will be described in the following. 
Due to Lorentz covariance the tensor integrals (A.21) can be decomposed into tensors 
constructed from the external momenta pi and the metric tensor g^y. The choice of the 
tensor basis is not unique and we will stick to the following form [187, 189]: 

B^^ pIBi 

B''"^ p'ip\B2,+g^''B22 

C"'^ = ptp'iC^i + P^P^C22 + {prnVC^^ + g^-'^C^A 

C"^^" = pi'pXCai +P»^C32 + {P1PIP2Y'"'C^^ + {PlP2P2Y'"'C^A 
+ {pi^r^C35 + {p2yr''C36 

Dl'- = P>tp'iD2l + P^2P2D22 + P^P^^23 + {PlP2}'^"^24 + {PlP3}'^"^25 

D^^<^ = pW.d^i +p'2P2f2D32 +P^3PsfsDs3 + {pipmY''' + {pwrnY'^'D^, 

+ {pmP2Y'"'D^^ + {pmPzY^'D^-j + {p2P2P3}'^'^''^38 + {p2P3P3}^'^''i^39 

+ [pmp^Y^'D^.o + {PigY^'D^ii + {P2gY'''D,,2 + {p3^}'^'^''^313 

(A.24) 
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where we have used the shorthand notations 

{PiPjPk}nup = ^ Pa(i)fiPaij)uPa{k)p 

cr{i,j,k) 

with a{i,j,k) denoting all different peYmutaiions of {i,j,k) and 

{Pi9} jxi/p — Pip,9vp + Piv9p,p + Pip9p.v ■ 

For example {piPiP2}^^^ = P1P1P2 ^ P1P1P2 ~^ PiPiP'2 ^^'^ {PiO}^"'^ = Pig"'^ + Pig^^ + plg^" ■ 
The scalar coefficients -Bj(fc), Cjk and Dj^i-^ can depend on all possible invariants (built 
from the leg momenta) and the masses rrii. 

Using the Lorentz decomposition (A. 24) of the tensor integrals all scalar coefficients 
can be iteratively reduced to scalar integrals (with equal or less propagator factors) [187, 
189], see [186] for a general treatment. With other words, all one-loop tensor integrals can 
be expressed in terms of scalar integrals (with equal or less propagator factors) . In simple 
cases the reduction can be done by hand, however, in general the reduction algorithm has 
to be automized on a computer.^ 

We now present some details of the derivation of the unrenormalized vertex in 
Eq. (A. 11). In Feynman gauge the vertex correction diagram in [Fig. 3.1(c.l)] is given in 
n dimensions by 

A/^_r^ (27r/x)^-" f^nu.j. + H+^i + ^i ^ , . 

^° " ^^4;^ in' J "^"^ {k+p2Y-ml+ k' (k+p.r-mf" ^ ^ 

where L5 = (^^ — ^75), Cf = 4/3 and the scale /i has been introduced in order to maintain 
the coupling gg dimensionless (sec above). 

In the following we consider on-shell fermions, i.e. pf = ml and P2 — m\, and a 
spacelike gauge boson with 4-momcntum q = P2 — Pi, (f < 0. Furthermore, since Aq 
will be eventually sandwiched between external spinors u{p2)A.Qu{pi) we apply the Dirac 
equations u{p2)'^2 ~ f^2u{p2) and ]/i{u{pi) = miu{pi) whenever possible which will be 
indicated by a '=' instead of a '='. Utilizing an anti-commuting 75, {7;i,75} = [190], 
the numerator (under the integral) in Eq.(A.25) can be written as 

Num = |7^ [{n - 6)k^ + 2[{k + pif - m?] + 2[(A; + ps)' - m^j + 4pi • P2] 

+ 2m2YH - 2^ [^Q.^ + (n - 2)k^ + 4p^] jLg - 2miYHR5 + 4miFi?5 
(A.26) 

^We are grateful to Dr. I. Bojak for providing us his Mathematica package for the reduction of the 
tensor integrals. This package is described in [189]. 
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with i?5 - (^ + ^75). 

Inserting the numerator into Eq. (A. 25) and utihzing the Passarino-Vehman- 
decomposition (A. 24) of the encountered tensor integrals (and again exploiting the Dirac 
equation in order to eliminate and there remain only Dirac structures as given in 
Eq. (A. 11). Comparing the corresponding coefficients we obtain 

-1 

(n - 6)Bo{q, mi, 1712) + 2Bq{pi, 0, mi] 



Co,. 



Cn.4- 



r(^) 



(n-4)/2- 



+2So(pi + q, 0, ma) + 4(m? + pi ■ q)Co 
-2m\Cii + 1{m\ - ml)Ci2 - 2{n - 2)024^ 

2mim2Cii 

-4Ci2-2(n- 2)^23 

-4Cii + AC12 -2{n- 2)(C2i - C23) 

-2{n - 2)C22 

-4Cii + 4Ci2 + 2(n - 2)(C22 - C23) 



(A.27) 



with Cijfi = Cij,o(pi, g, 0, mi, m2). 

As described above the coefficients of the Passarino-Veltman-decomposition can be 
expressed in terms of the scalar integrals 



DCii 
DC12 
DC21 
D C22 



^ ^^-[So(p2, 0, m2) - Bo{q, mi, m2)] - ^^^[-Bo(pi, 0, mi) - Bo{q, mi, m2)] 



2 

Pi ■P2 



2 [^o(P2, 0, m2) - 5o(g, mi, ms)] + ^[5o(pi, 0, mi) - B^^q, mi, m2)] 



^o(^i)-r-r-^o(m2)- 



4mf 



4m^ +5o(g,mi,m2)g2(^_2^ 



= (Ao(mi) - ^0(^2)) + ^(^2) 



4g2 



Pi ■ g + mf 

4m2 



+So(g,mi,m2) 



mt 



2(n-2) 



Pi • q ^ ijnl - ml)pi ■ q 



-1 



A(^i)^ + ^o("^2) 
1 



i^C23 

C24 = Bo{q,mi,m2) 



Pi-q + mf 



4g2 

Bo{q,mi,m2)pi ■ q 



n — 3 
2{n - 2) 



2(n-2) 



23 



with D — m\q — {pi- q) ■ It is easy to convince oneself that the coefficients Cn, ... ,Q 
are finite, whereas C24 is UV-divergent. Furthermore, Co(pi, q, 0, mi, m2) is IR-divergent. 
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At this place the problem has reduced to the calculation of the scalar integrals Aq, Bq 
and Cq. Introducing Feynman parameters and performing a Wick rotation the iV-point 
scalar integral can be brought into the following form: 



^ T^{p,, . . . mo, ... , m^r-i) = ^^^^ / d'^^ E ^ 

= (-l)^r(iV-n/2) (4V)^ [\daU {M'y^'-'' (A.28) 

[(iajAT = dai ■ . . . ■ da^ 5 11 — j 

N 

= - ^ Oiiajplj + ^ ajm- - is (A.29) 

N>i>j>l i=l 

where pij has been introduced in Eq. (A. 23). 

The calculation of the scalar one-point function AQ{m) is trivial: 

^(m) = -m^ j ' r (l - ^) = (^Aw - In ^ + 1^ + 0(n - 4) (A.30) 

with the UV-divergence contained in 

_2 

Auv = 7 - 7e + In 47r (A.31) 

n — 4 

and 7e is Euler's constant. Note that Eq. (A.30) imphes Aq{0) — 0. The terms of order 
0{n — 4) are only relevant for two- or higher-loop calculations. 

The evaluation of the needed two-point functions is still easy 



with 



, . I mim2 rril — rni m2 
Bo{q,mi,m2)\q2<o = Auv + 2 - In \ In — 



(A.32) 



with 

^2 _ (mi +m2f - q^ 
(mi — m2)^ — q^ 

and 

,2 



m 

Bo{p, 0, m) |p2=^2 = Auv - In — + 2 . (A.33) 

A* 
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Finally the required three-point function is given in n = 4 + 2s dimensions by (M^ = 
m\x + m|(l — x) — q^x{l — x) — is) 

Coipi, q, 0, mi, m2)\p2=rnl{pi+q)^=mW<o = (47r//^)~^ r(l (A.34) 

with 

7= /" dx[mlx + ml{l-x)-q^x{l-x)\-^+^ . (A.35) 

The Feynman parameter integral / is finite at £ = since X = = mfx + m|(l — x) — 
q'^x{l — x) > ioT X e [0, 1] as long as > and m| > 0. Therefore we can expand I 
in e up to neghgible terms of the order 0{e'^) 

7 = f dx = f dx + = /i + eh + ■ (A.36) 

Jo Jo ^ 

The integral Ji = Q dx ^ is elementary and can be found in Eq. (A. 14). I2 = Jq dx 
can be calculated by writing X in linear factors and by partial fractioning the denominator. 
After some algebra one finds the result given in Eq. (A. 15). 



A. 3 Real and Virtual Contributions to Structure 
Functions 



The soft real contributions Si to the coefficient functions in Eq. (3.10) are given by 

2A \ ^ / 2A 



S, = i(-2 + E++7i) + 2 + % 
e A 



A 7i + Lis 



A + E++ 



-•2,3 



+ In 

= Si 



(-2 + E++7i) 



(A.37) 



with 7i given in Appendix A. 2 and where x is given in Eq. (3.7). The virtual contributions 
are derived from the renormalized vertex in Eq. (A. 19) by using the projectors in Eq. (A. 3): 

- 2S+mim2 

Vl - L.R + — — — 0+ 

o+2^_,_+ — 2S^mim2 



1 S'_ 
V2 = Cr,^ + - {ml Ci,+ + ml + — 



C, + ^(C,, + C,.) 



(A.38) 
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where the Cs are given in Appendix A. 2. Note that the soft poles {l/e) of Si, Vi cancel 
in the sum Si + Vi in Eq. (3.10) as must be. 

The massive matrix elements f^{C) ^^^'^ most conveniently given as functions of the 
Mandelstam variable si(C') = (pi + qY — ml, i.e. f^{C') = f^[si{i')] with 



(A.39) 



From the real graphs of Fig. 3.1 (b) one obtains 



E++ + §1 ^ si + ml " 



+ ^4 



2 A'si 

2-++ (A2 + 4m^E+_)- 



-m^E 

{Si+ml)Si''" ' 4(si + 7712 



A" + 2E+_E++ + {ml + 
1 



3E;^E_+ + 4777^(10E++E^ 



- E+_E_+ - 777fE++) + Si[-7E++E_+ + 18A^ - 4777f(7g^ - 4777^ + 7777^ 



+ 3s^[E+_ - 2777^] - s\ 



+ 



Si + 7772 

2A' 



Sl 



■E++ (A^ + 2E+_E++) 



+ (4777^777^ - 7E+_E++) - 4E+_Si - S^] L^/ 



16 



f^ih) = - 2 AX^^,+ 2777177725. 1^ 



h + ml / 2 



6777 



+ 
+ 



^'ol^'l + (2^" - + ^-)) I + f - - 6m?Q^)(si + 777^) 

2(si + 7772) V / \ 

2 (777? + 777^) §1 - 9777^E^_ + A^ (2E++ - 777^) + 2si {2A'^ + {ml - 5ml) ^+-) 

(A'2 - 6Q2 {rnl + h)) S++ {§1 + 2A' 



2{si + ml) 



Sl 



(A^ + 2(2777^ + si)E+_) 



(si + 7771) ^^^2(^2 ^ 2E^_E^^) - 2si(A^ - 6777?g^) 



Sl 



- (A'^ - 18mjQ^)J:++ - 2A^ (E++ + 2E+_)] L^, 



16 



fsi-^l) = -2AX/^,+ 2777i7772i?_ ( 1 



+ 



(si + 777I) (si + E_, 



+ E_ 



- 3E+_ - ^ (A^ + 2777lE+_) - 



Sl A'si 

5i-E_)(,si + E, 



2(si + 777i) 



+ 



Sl + 7772 
A'Si 



[-§1 + 4 (777?E_+ - A^) - 3siE+_] L^, 



(A.40) 
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with 

E++ + si - A' 



L^i = In 
and 



S++ + si + A' 



- i2 + A/S2 ^++ 



A is given below Eq. (3.4) and A' = A[m?, s, -Q^]. 
Finally, the normalization factors in Eq. (3.10) are 

1 2A (A')2 ' A' ^ ^ 

A. 4 Comparison with E. Hoffmann and R. Moore, 
Z. Phys. C20, 71 (1983): 
Detailed List 

HM: Refers to E. Hoffmann and R. Moore, Z. Phys. C20, 71 (1983) 

KS: Refers to S. Kretzer and 1. Schienbein, Phys. Rev. D58, 094035 (1998), as 

documented in this thesis in Chapter 3 and the above Appendices A.1-A.3. 

Replacements: 
mi,2(KS) —>■ m 
V,A{K^) ^1,0 

As in HM: A = mVQ^. 

The bracketed numbers refer to the equations in the Hoffmann & Moore article and in 
this Thesis, respectively. 

Vertex Correction 



A^ = r'^ 

(A.19) = (31) 

^ KS ^ ^ H 



up to a Gordon decomposition 
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Virtual Contributions 



2-K 



Cf [A(A.38)]ks 



IV 



-2 v^rriA a^^^ 



IV 



= -2 V1 + 4A (39) 



VI + 4A 6{l-z) 



Soft Real Contributions 



V1 + 4A 
V1 + 4A 



cr 



+ /(A) 



5(1 -Z) 



(48), 



+ /(A) 



5(1-^) 



where the discrepancy /(A) 



/(A) = ^ 2 HI + 4A) (-1 - iu 

^ ^ 27r ^ ^ ^ V yrriA vmA + 1 



0; A^O 



arises from the different choice of the convolution variables and is exactly cancelled if the 
identity 



is used in Eq. (48) or in Eq. (3.10) 



Hard Real Contributions 



Pis 

2^ 



-"2 



a 



2. P-IDksI^.,, 



1 + 4Az^ 

V1 + 4A ^7^1 

1 + 



^ = v^rriA (-2af)+af)) 

27r 



a 



2. 



Vl +4A [-2 (52) + (51) II ^, 

L V / HM ^ HM J 10^1 

where (52) should be multiplied by an obvious factor z. 
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Summary 

We (KS) agree with HM up to: 

• The sign of the virtual contribution to the longitudinal structure function (which is 
irrelevant for F|) 

• The overall normalization of F^. According to the above equations our (KS) nor- 
malization differs from HM by a factor Vl + 4A. This does not really mean a 
discrepancy but stems from the choice of normalizing the 0{a],) coefficient function 
relative to a Born term contribution of 5(1 — z) (HM) or 5(1 — ^') (KS), correspond- 
ing to (37) or (3.8) . However, we regard our normalization as preferable since 
it corresponds according to (3.2) directly to the notion of c(^) carrying a light 
cone fraction ^ of the nucleon's momentum and is therefore directly related to 

— the original derivation of the conventionally used intrinsic charm component 
in Phys. Rev. D23, 2745 (1981) {Intrinsic heavy-quark states) 

— the all-order proof of Hard-scattering factorization with heavy quarks in Phys. 
Rev. D58, 094002 (1998) 

• The normalization of which is -besides the y/1 + 4A normalization ambiguity- 
wrong by a factor (1 -|- 4:Xz^)/{l + 4 A) in HM. The HM expression corresponds 
to the partonic structure function multiplying the ~ p'^p'^2z'^ / Q'^ partonic tensor. 
This normalization is only correct in the massless case and neglects finite mass 
terms which arise when contracting with the leptonic tensor and extracting the 
structure function as the properly normalized coefficient of the ~ {1 — y) term. 
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Matrix Elements for Real Gluon 
Emission off Massive Quarks 



The projections of the partonic Matrix Element onto the structure functions are most 
conveniently given in the Mandelstam variables si and ti which are defined below Eq. 
(4.6): 



A'2 

+ 2mim2S. 



SI 



i\ Siii 



'mlsiih + 2E+_) ^ A'^ + (mi + Q'^)si + 2E+_E 



++ 



^1 



A'^ - si(mf + + si) + 2miS+_ . {ml + Si] 
+ ^ ^1" 



Sl 



Sl 



mi2siS+_(si + 2E++) -sl - 4s2S+_ + si{Amlml - 7E+_E++) 
+ '-'+ ^ 1 



n 



^ 2E^^(-A^-2E^_E^J ^ 



2ti 

(A2 + 2miE+_)E 



++ 



16 
A^ 



- 2A^5, 



2ti 

4mi + 2mjsi - ^+-{ml + E+_) 
A'^-2(mi + Si)E^ 



-'++ 



+ ^ + 



2si 

+ 2mim2<S'_ 



(A^^ - 6mfg^)5i 



+ 



2{A"-3Q■'{s^ + E^^)) 



Sl 
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+ 
+ 



-2m2si[(A2 - 6m^g' )gi + 2A%_] ^ -2A^( A^ + 2E+_E++) 



^1 

-si[2(A2 - QmlQ^Ysi + (A'^ - l^mlQ^)T.++ + 2A2(3S++ - Ami) 



ti 

.9 9x 9 9 9 2A2(A2 + 2m^E,_) 
- 2(m^ + ml)sl - 9mlEl_ ^ + 2si[2A2 



Sl 



+ (m^-5m2)E+_]+A2(2E++-m^) 



++ 



16 



A'2 
+ R 



ml mf E^ 



[A"-6Q\ml + h)]i: 

Sl 

Sl + E+_ ^ Sl — E_ 



+ 2mim2R- 



ti 



+ 



-2mpiE+_ ^ -s^ - 4(A^ - m?E_+) - 3siE^ 



''1 



+ 



2{m\ — ml) 



2(A2 + miE+_) 



+ ti 



Sl — E_ 

Sl 



(B.l) 



where we conveniently use the shorthands A = A[m^, m|, —Q'^] and A' = A[m^, s, — Q^]- 
In order to obtain the inclusive results f^ih) in Appendix A the f^{si,ii) have to be 
integrated over < y < 1, i.e. 



Jo 



dyf?(si,ii) , 



(B.2) 



where y is defined via the partonic center of mass scattering angle 9* and related to ti 
through 



y = - (1 + cosr) 



Sl 



Sl + 7712 



A'{y-yo) , 



(B.3) 



with yo — [1 + (E+_,. + si)/A']/2 being the would-be coUinear pole of the ^-channel 
propagator. 



Appendix C 

Limits of the Doubly Virtual Box 



Due to the particular choice of variables, the results in Eq. (8.11) are especially well suited 
for deriving various important limits wide-spread over the literature which are relevant 
for some of the chapters in this thesis: 

• Most important for our purposes is the Bjorken limit {Q^ oo, u oo, x = const) 
in which we can study structure functions of the real and virtual photon as has 
already been discussed in Chapter 7. 

• Also needed is the real photon P'^ — {6 — 0) case in which the general virtual 
box results in (8.11) reduce to the standard box-diagram 7*((5^)7 — QQ expressions 
for a real photon 7 = 7(P^ = 0). Keeping the full mass dependence in (8.11) we 
obtain expressions relevant for the heavy quark contribution to the photon structure 
functions. 

• Finally, the general light quark mass limit (for arbitrary P^, Q^) can be easily 
obtained from (8.11) by setting m — {X — 0, (5 — 1) and needs no separate 
discussion. 

In the following, our main concern will lie on the Bjorken limit. Practically this limit 
means that is much larger than the other scales and P^. Beside the general case 
m^, P^ ^ which will be studied in Section C.l additional orderings = 0, -C 
(Section C.2) and P^ = 0, -C (Section C.3) are of interest and the expressions 
further simplify under these circumstances. One must be careful in handling these limits 
because for some of the box cross sections the result depends on which of the two limits 
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m — > and — > is taken first. For example, below we will find the following asymptotic 
expressions for (Jtl derived from (8.11): 

= lim lim cttl 7^ lim lim cttl = ^c^n 4x^(1 — x) 
Mathematically, the origin of such a behavior is easily identified. Terms like 



Ax6 + \i3^ 

occurring in (8.11) (not only as the argument of the logarithm) require a careful treatment. 
In general they are not negligible even for small and and, viewed as a function of 
and P^, they are discontinuous at (m^,P^) = (0,0). Finally, in Section C.3 we also 
derive expressions for the heavy quark contributions to the photon structure functions in 
the real photon limit by setting = (5 = 0) but keeping the full mass dependence in 
(8.11). 

All results for the photon-photon cross sections will be given for a single quark with 
charge Cg and mass m. The photon structure functions F^, (i = 1, 2, L) can be obtained 
from these expressions with help of the relations in Section 7.2 which simplify in the 
Bjorken hmit (P^ < ^ ;9 ~ 1): 

^2^' = r^^2L , = -^aLL (C. 1) 

with (72T = ctt + clt, = ctl + cll- The commonly utilized expressions for a spin- 
averaged target photon are given by 

^<7>^^7T_1^7L^(i^l^2,L). (C.2) 

Finally, the structure function can be deduced from F^ = F^—^xF^. Since the structure 
functions arc (apart from the normalization factor Q"^ /A-k'^o) simple linear combinations 
of the photon-photon cross sections they will only be written out in those cases where 
they are needed explicitly in the thesis. 
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C.l General Bjorken Limit: m?,P'^<^Q'' 



In the general Bjorken limit the normalization factor N given in Eq. (8.10) and the 
logarithm L from Eq. (8.9) are given by 



(C.3) 



L = In 



4:XS + X 

Keeping this in mind and using P — 1 + 0{X), P — 1 + 0{S) the following results can be 
easily deduced from Eq. (8.11): 





= N4 










ttt 




'TT 





Atto^ 



-X 



x^ + {1- xf L + Ax{l -x)-l 



Ax5 



AxS + X 



+ 0{S, A) 



AxS 



-X 



^ ^Ax5 + X 
{-2x^ + 0{5, A)} 



+ 0{5, A) 



X < {2x - 1)L + 3 - Ax 



AxS 



+ 0{d,X)} ■ 



AxS + X 

The remaining expressions are suppressed by powers oi 5 — xP'^/Q'^ (see Eq. (8.11)) 
At the price of a slightly worse approximation (at larger x) one can further use 



(C.4) 



AxS 



P^x(l-x) 



AxS + A P^x{l - x) + 
Q^l-x) 



L = ln 



x[P'^x{l — x) + w?] 



+ 0{^,X), 



(C.5) 



For example we can write: 



-X 



Q 

Ph-il ,r) 



l-2x(l-x] 



In 



Q\l-x) 



x[P'^x(l — x) + m^] 



+ Ax{l - x) -1 



P'^x{l — x) + m? 
.2 _ n d2 ^ /n2 



+ 0{^,,X)\. 



(C.6) 



C.2 = 0,P' <^Q' 



In this section we consider the asymptotic virtual (P^ ^ 0) box expressions for light quarks 
in the Bjorken limit for which the expressions in Eq. (C.4) further reduce. Noticing that 
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the logarithm L is given by L = In we can write in this case (neglecting terms of the 
order 0{5)): 



4 ^ 

(7TT ^ NA^^x { [x" + (1 - xf\ In + 4x(l - x) - 2 



(7tl ^ (7lt ^ Nee ——X [4:x{l - x)] 



Ttt ^ce;-^:r [-2a:2] 



TTT^N^e^.^x <j(2a;-l)ln-p|-^ + 2-4a: 



(C.7) 



Summing over q = u,d,s and utilizing Eqs. (C.l) and (C.2) we recover the well known 
asymptotic results for the virtual (P^ ^ 0) box structure functions for the light q = u,d,s 
quarks in the Bjorken limit P'^/Q^ <^ 1: 

F^j>ii^, Q\ p') ^ E 4 ^4 ^ " ^^'^ & ^ ^^^^ " " 4 ^^-^^ 

g2,p2)~Ar^EeJ^x|[x2+(l-x)^] ln-^ + 6a;(l-a;)-2| . (C.IO) 



C.3 Real Photon Limit: = q 



For p2 = the virtual box results in (8.11) reduce to the standard box-diagram 
7*(Q^)7 — > qq expressions for a real photon 7 = 7(P^ = 0). 

4m? X 



The heavy quark contribution becomes, utilizing 5 = 0, /3 = 1 and A 



, 4-7rrv^ 



+ (1 — xY + x{l — X 



Ami 28ml 



X 



In 



" Q^(i-xy 

1-/3 



Ami 



Ana' 



A'KOp' 



Ax{\ — x) — 1 — x{\ — x) 

Q 

xQ{0^)i-x^^-^\n 



ml, 1 + /5 



TTT = N^el-^x e(P') 



—X 



l-(3 
28ml ^Sml 



+ /3 Ax{l - x] 



— X 



1 ^ + ^ R 



2x^ + x(l-x)^ 



r^T = AT^et^x e(/3^) {(2x - 1) In + /3(3 - 4x)} 



(C.ll) 
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i.e., according to (C.2) (or (C.l)) 



' TT 



a 



Ami 



+ (1 - xf + x{l - 2>x) 



8mi 



2 "'"-/i 



X 



In 



1-P 



8x{l-x)-l-x{l-x)^ 

Q 



F2U^, Q') - 3 ei %Q{/3') In 1±| + /? 4x(l 



X 



(C.12) 



which are the famihar massive Bethe-Heitler expressions [117] relevant for the heavy 
quark contributions to the structure functions of real photons (cf. [6], for example). 

In the hght quark sector where A <^ 1, i.e. = <^ Q^, the logarithm can be 
written as 



l-l3 



rrigX 



and we obtain from (C.4) or (C.ll) the following results (neglecting terms of the order 

o{x)y. 



axT^iVeeJ^o; |[a;^ + (l-a;f]ln 
(7lt ^ N^e ——X [Ax{l - x)\ 



Q \l-x) 
rrigX 



+ 4a;(l - x) - 1 



TTT ^ A^ceJ^x [-2x'] 



47ra2 f_ g^^l-x) 



T^T - Ncet^x \ {2x - 1) In 



niqX 



+ 3-4:X} , 



(C.13) 



i.e., according to (C.2) (or (C.l)) 

F^l^ix, Q') ^N^J24 -4 + (1 - ^)'] 1^ ^^^^ + 8x(l - x) - 1 1 . (C.14) 
Note also that cttl vanishes like cttl oc P'^/m^. 



Appendix D 
Parametrizations 

D.l Pion Distributions 

D.1.1 Parametrization of LO Part on Distributions 

Defining [26] 



to be evaluated for /Xlo ~ 0-26 GeV^, all our resulting pionic parton distributions can 
be expressed by the following simple parametrizations, valid for 0.5 ^ ^ 10^ GeV^ 
(i.e. 0.31 < s ^ 2.2) and 10~^ ^ x < 1. For the valence distribution we take 



s = 



in 



ln[gV(0.204GeV)^] 
ln[/.2Q/(0.204GeV)2] 



(D.l) 



xv''{x,Q^) = N x^'il + A^/^ + Bx){l - x) 



D 



(D.2) 



with 



N 



1.212 + 0.498 s + 0.009 



a 



0.517 -0.020 s 



A 



-0.037 - 0.578 s 



B 



0.241 + 0.251s 



D 



0.383 + 0.624 s. 



(D.3) 
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The gluon and light sea-quark distributions are parametrized as 

b 



For w — g 



c" {A + By/^ + Cx) ^Ini^ 



+ exp -E 



E'snn 

X 



a = 0.504, 

a = 2.251 - 1.339^5, 

A = 2.668- 1.265 s + 0.156 



/3 = 0.226, 
6=0, 

B = -1.839 + 0.386 s, 



C 

E 



-1.014 + 0.920s- 0.101 s^, D 



E' 



1.245 + 1.833 s, 

and for the hght sea w — q 

1.147, (3 

0.309- 0.134^, b 

0.219 - 0.054 s, B 

1.100 - 0.452 s, D 

4.521 + 1.583 s, E' 

The strange sea distribution s'^ = s is parametrized as 



a 

a 

A 

C 

E 



-0.077 + 1.466 s, 
0.510 + 3.844 s, 

1.241, 

0.893 - 0.264 ^s, 
-0.593 + 0.240 s, 
3.526 + 0.491s, 
3.102. 



xs 



"(^' ^ ITTVa (1 + + Bx) (1 
(In-)" 



with 



a = 0.823, (5 

a = 1.036 -0.709 s, A 

B = 5.580- 1.281s, D 

E = 5.101 + 1.294 s, E' 



x)^ex.p I -E+ \l E'snn^ 



0.650, 

-1.245 + 0.713 s, 
2.746 -0.191s, 
4.854 - 0.437 s. 



(1- 



(D.4) 



(D.5) 



(D.6) 



(D.7) 



(D.8) 
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D.1.2 Parametrization of NLO(MS) Parton Distributions 

Defining [26] 

ln|QV(().299GeV)T 
lBtt.S„o/(0.299GeV)^] ^ 

to be evaluated for /^nlo — O-^O GeV^, our NLO predictions can be parametrized as the 
LO ones and are similarly valid for 0.5 ^ ^ 10^ GeV^ (i.e. 0.14 s <, 2.38) and 
10^^ ^ X < 1. The valence distribution is given by (D.2) with 

N = 1.500 + 0.525 s- 0.050 
a = 0.560 - 0.034 s 
A = -0.357 - 0.458 s 
B = 0.427 + 0.220 s 

D = 0.475 + 0.550 s. (D.IO) 
The gluon and light sea distributions are parametrized as in (D.4) where ior w — g 
a = 0.793, P = 1.722, 

a = 1.418 - 0.215^/s, 6=0, 

A = 5.392 + 0.553 s - 0.385^2, B = -11.928 + 1.844 s, (D.U) 
C = 11.548 - 4.316 s + 0.382 s2, D = 1.347 + 1.135 s, 
E = 0.104 + 1.980 s, E' = 2.375 - 0.188 s, 

and for the light sea w — q 

a = 1.118, P = 0.457, 

a = 0.111 - 0.326^5, & = -0.978 - 0.488^5, 
A = 1.035 - 0.295 s, B = -3.008 + 1.165 s, (D.12) 
C = 4.111 - 1.575 s, D = 6.192 + 0.705 s, 
E = 5.035 + 0.997 s, = 1.486 + 1.288 s . 

The strange sea distribution is parametrized as in (D.7) with 
a = 0.908, p = 0.812, 

a = -0.567 - 0.466 s, A = -2.348 + 1.433 s, 

(D.13) 

B = 4.403, D = 2.061, 

E = 3.796 + 1.618 s, = 0.309 + 0.355 s . 
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Let us recall that in the hght quark sector — — — dl , — — 
u""' = d''' and = (/"^ + /"")/2. 

D.2 Photon Distributions 

Simple analytic parametrizations in LO and NLO of the ' hadronic' piece of the real and 
virtual photonic parton distributions, being proportional to f'^{x,Q'^) in Eqs. (11.5) and 
(11.13), respectively, are already known according to the parametrizations for f'^{x,Q'^) 
in Appendix D.l. Therefore we only need to parametrize the remaining 'pointlike' com- 
ponents in Eqs. (11.5) and (11.13). 

D.2.1 Parametrization of LO 'pointlike' Photonic Parton Dis- 



In LO the dependence of the 'pointlike' /Jj(a;, Q^) term in (11.5) enters, apart from an 
overall l/agiQ"^) factor, merely via the combination L = as{Q^) / Oisil^to) evident, for 
example, from Eq. (2.12) in Ref. [112]. Therefore, we prefer to parametrize the quantity 
fZ{x, Q^) in terms of 



where [26] hIq — 0.26 GeV^, which will later provide us a parametrization also for the 
virtual 'pointlike' component in (11.13). Our resulting 'pointlike' distributions in Eq. 
(11.5) can be expressed by the following simple parametrizations, valid for 0.5 ^ ^ 



tributions 



s = In 



In [gV(0.204GeV)2] 
ln[/x2o/(0.204GeV)2] 



(D.14) 



10^ GeV^ (i.e. 0.31 <s< 2.2) and 10"^ < 



X < 1 : 




9 , Q' 



2 s°x"(A + Sv^ + Cx'') 



47r (0.204 GeV) 




(D.15) 
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where for /Jj = = 

a = 0.897, a' = 2.626, 

(3 = 0.413, 

a = 2.137- 0.310 ^/i, 6 = -1.049 + 0.113 s, 

A = -0.785 + 0.270^5, 5 = 0.650 - 0.146 s, 

C = 0.252 - 0.065 Vs, D = -0.116 + 0.403 s - 0.117 

£; = 6.749 + 2.452 s- 0.226 s^, E' = 1.994 s - 0.216 s^ , 

for /Ji = = (i^j = s^i = s^i 





= 1.084, 


a' 


= 2.811, 




= 0.960, 






a 


= 0.914, 


b 


= 3.723 - 0.968 s, 


A 


= 0.081 - 0.028^5, 


B 


= -0.048, 


C 


= 0.094- 0.043^5, 


D 


= 0.059 + 0.263 s- 0.085 s^ 


E 


= 6.808 + 2.239 s- 0.108 s^. 


E' 


= 1.225 + 0.594 s- 0.073 s^ 



(D.16) 



(D.17) 



(D.18) 



and for f^^ = g^^ 

a = 1.262, a' = 2.024, 

(3 = 0.770, 

a = 0.081, b = 0.848, 

A = 0.012 + 0.039^5, B = -0.056 - 0.044 s, 

C = 0.043 + 0.031s, D = 0.925 + 0.316 s, 

E = 3.129 + 2.434 s- 0.115 s^, E' = 1.364 + 1.227 s - 0.128 s^ . 

With these parametrizations at hand in terms of s in (D.14), the appropriate ones 
for the 'pointhke' distributions /j/^ \x, Q^) of a virtual photon appearing in Eq. (11.13) 

arc simply given by the same expressions above where in (D.14) /i^o has to be replaced 
by = max(P^, /iLo)- As discussed in Sec. 11.4, these parametrizations can, within 
sufficient accuracy, also be used for the parton distribution of virtual photons in NLO. 
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D.2. 2 Parametrization of NLO 'pointlike' Photonic Parton Dis- 
tributions 

In NLO the dependence of the 'pointhke' distributions of real photons in (11.5), 
/pi(x, Q"^), can be easily described in terms of the following 'effective' logarithmic ratio 

In [QV(0.299GeV)2] 



In 



(D.19) 



ln[^2^Lo/(0-299GeV)2] 

0.40 GeV^. Our NLO(DIS^) predictions can now be 
parametrized as the LO ones and are similarly valid for 0.5 ^ ~ 10^ GeV^ (i.e. 0.14 
^ s ^ 2.38) and 10~^ ^ a; < 1. For convenience we include now the NLO a, in the r.h.s. 



to be evaluated for /i^ 



NLO 



of (D.15), i.e. 



a 



+s°'exp I-E+ ^ E'sP In ^ 



[1-x)^ 



where for /J; = = li^j 

a = 1.051, 
(3 = 0.970, 



a 
A 
C 
E 



0.412-0.115^5, 
-0.028 + 0.019 s^, 

6.726 - 3.264 v^- 0.166 
4.122 + 3.170 s -0.598 s^. 



b 
B 
D 
E' 



= 2.107, 



for r 



pi 



"pi 



"pi 



'pl 



"5 pi 



a 

a 
A 



1.043, 
0.457, 

0.416-0.173 v^, 
-0.010^/^ + 0.006 



C = 1.577- 0.916 Vs, 



a' 

b 

B 

D 



E 



5.240 + 1.666 s- 0.234 E' 



4.544 -0.563 s, 
0.263 + 0.137 s, 
1.145-0.131^2, 
1.615 s- 0.321 s2. 



1.812, 

4.489 - 0.827 s, 

0.064 + 0.020 s, 

1.122 - 0.093 s -0.100 s^, 

1.284 s, 



(D.20) 



(D.21) 



(D.22) 
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and for f^^ = 



a = 0.901, 

(5 = 1.666, 

a = 0.844- 0.820^/8, 

A = 0.194, 

C = 0.330 - 0.177 s, 



a' 

h 

B 

D 



E = 2.895 + 1.823 s - 0.441 E' 



1.773, 

2.302 -0.474 s, 
-0.324 + 0.143 s, 
0.778 + 0.502 s -0.154 s^, 
2.344 -0.584 s. 



(D.23) 
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